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Summary 


The known results for the propagation of sound in plates with Poisson’s constant o=0.5 
(liquids) and o=0.35 (aluminium) are completed by dispersion curves with o=0.47 (rub- 
ber). The dispersion curves are easily understandable because of the large difference 
between the velocities of transverse and longitudinal waves occurring at this value of o. 

i Measurements of the propagation of elastic waves in rods of rubber with square 
( cross-section in the frequency range 0.1 to 300 kc/s are reported. The amplitude is meas- 

ured as a function of distance by partly submerging the rod into a waterfilled vessel within 
which a hydrophone is arranged. 

The results with respect to phase velocity and attenuation can be explained in terms 
of the propagation of extensional waves at low frequencies and compressional waves at 
high frequencies separated by a range with high attenuation. 


Zusammenfassung 


Die bekannten Ergebnisse fiir die Schallausbreitung in Platten mit einer Poissonschen 
Konstante von 6=0,5 (Fliissigkeiten) und o=0,35 (Aluminium) werden durch die Disper- 
sionskurven fiir o=0,47 erganzt. Durch den grofen Unterschied zwischen Transversal- und 
Longitudinalwellengeschwindigkeit, der bei diesem Wert der Poissonschen Konstante auf- 
tritt, werden die Dispersionsverhiltnisse tibersichtlicher. 

Weiterhin wird iiber Messungen der Schallausbreitung in Gummistében mit quadrati- 
schem Querschnitt im Frequenzbereich von 0,1 bis 300 kHz berichtet. Die Schwingungs- 
amplitude langs des Stabes wird dadurch ermittelt, dai dieser allmahlich in ein wasser- 
gefulltes Gefa8 eingetaucht wird, in dem sich ein feststehendes Wasserschallmikrophon 
befindet. 

Die MefSSergebnisse fiir Phasengeschwindigkeit und Dimpfung lassen sich durch das 
Auftreten von Dehnwellen bei tiefen und von Kompressionswellen bei hohen Frequenzen 
erklaren. Die Ubertragungsbereiche sind durch einen Frequenzbereich hoher Dampfung 
yoneinander getrennt. 


Sommaire 

Les résultats connus concernant la propagation du son dans des plaques ayant un 
coefficient de Poisson o=0,5 (liquides) et o=0,35 (aluminium) ont été complétés par les 
courbes de dispersion relatives au cas =0,47. Les caractéristiques de dispersion sont bien 
mises en évidence du fait de la grande différence qui existe, pour cette valeur du coefficient 
de Poisson, entre les vitesses des ondes transversales et longitudinales. 

On a également mesuré la propagation du son dans des barreaux de caoutchouc de sec- 
tion carrée, dans la gamme de fréquences 0,1 4 300kHz. On a déterminé l’amplitude 
vibratoire le long du barreau en le plongeant peu a peu dans un récipient rempli d’eau, 
et ou se trouvait un microphone fixe. 

Les résultats des mesures de la vitesse de phase et de l’amortissement peuvent s’expli- 

* quer par l’apparition d’ondes d’allongement aux basses fréquences et d’ondes de compres- 
sion aux fréquences élevées. Les domaines de transmission sont séparés par une gamme 
de fréquences ot l’amortissement est plus grand. 


1. Introduction plate or rod respectively on the sound propagation 

within it. In the first case care must be taken to 

The propagation of elastic waves in plates and obtain a simple relation between the propagation 
rods has very often been investigated theoretically constants and the elastic constants of the material. 
and experimentally. In these investigations the main For instance from velocity and damping of extensio- 
interest was either in the determination of the elastic nal or flexural waves in thin plates and rods (wave 
constants of the material of the sample, or in finding _ length at least 6 times as large as the thickness of 
the influence of the geometrical dimensions of the _ the plate or the linear dimensions of the cross-section 
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of the rod resp.), Young’s modulus and its loss tan- 
gent can be calculated by means of simple formulas. 

Measurements and calculations of the propagation 
of elastic waves in plates and rods, aimed at in- 
vestigating the influence of geometry, have been 
made up to now only with metal plates and metal 
rods with rectangular, square or circular cross- 
sections [1], [2], [5], [6], [7]. The use of metal 
for the investigations has two substantial advantages: 
lst their loss tangent is very small (mostly smaller 
than 10-3) and therewith allows a very exact deter- 
mination of resonant frequencies. 2nd their elastic 
properties are mostly independent of frequency and 
therewith the velocities of the pure transverse and the 
pure longitudinal wave also. In this case the 
phase velocity of a certain type of wave in samples of 
equal shape of cross-section is a function of the ratio 
d/4 alone, where d means a linear dimension of the 
cross-section, e. g. the diameter in the case of cylin- 
drical rods or the plate thickness respectively, and 4 
the wave length. 

These advantages are absent when measurements 
are made with materials having large loss tangents 
and elastic properties which alter with frequency. 
Such materials are, for example, rubber-elastic ma- 
terials such as soft rubber and numerous plastics. 
Measurements of the sound propagation in such 
materials are reported in many publications too, but 
there the interest was nearly exclusively directed 
towards the determination of the elastic constants. 
The present paper mainly deals with the influence 
of the geometry of the samples on the propagation 
of elastic waves. When compared with the con- 
ditions in metals substantial differences are given 
— in addition to the frequency dependence of the 
elastic constants and the large loss tangents already 
mentioned — also by the fact that the numerical 
values of Lamé’s constants /* and wu differ by several 
orders of magnitude from each other (for soft rub- 
ber 4*/u=4000, while for aluminium 4*/w = 2.3). 
The corresponding loss tangents moreover largely 
differ from each other, and always 7i<7,. 

The damping of elastic waves in general is de- 
termined by both these loss tangents 7, and 7, . For 
some simple cases (extensional wave, plane shear 
and longitudinal waves) their influence may be ex- 
pressed by a loss tangent of the experimentally de- 
termined elastic constants as e. g. Young’s modulus 
and “sound stiffness” as discussed in details by 


Sxuprzyk ([8], p. 771). 
2. Theoretical considerations on wave propaga- 
tion in plates 


For the determination of the propagation of elas- 
tic waves in solid plates of infinite extension, the 
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assumption of obliquely propagating wave fronts 
has proved very suitable. This model yields exact 
results for propagation in isotropic media with 
negligible losses. The latter fact, of course, cannot 
be assumed in our experiments but in spite of that, 
this method may be mentioned here. A detailed 
representation has been given e. g. by Scuocu [6]. 
The conditions are simplest when the considera- 
tions are first made for a medium without shear 
tensions («=0) e. g. for a liquid [9]. The boundary 
conditions, demanding vanishing sound pressure at 
the two plane surface — having a distance d from 


each other which is-equal to the thickness of the — 


plate -- lead to the following equation for the angle ? 
between the plate-normal and the wave-normals of 
the two symmetrically incident waves: 


sin E d cos a cos (" d sin 0) =0 
io Z 
with the solutions 


for 2=0, 15 233 


kd 
From the angle of incidence of the obliquely pro- 
pagating waves we immediately get the trace velo- 
city of the waves which means the phase velocity 
of the “plate-wave”’ along the surfaces 


v=c/sin a-e/ V3 = ( He 
kd 


where c is the velocity of the longitudinal waves, 
and é their wave number. For the phase velocities 
of the different types of waves we thus get the 
hyperbolas represented in Fig. 1. 


0 1 2 3 4 
fd/o——= 


Fig. 1. Phase velocity v of free waves in a “plate” 
without shear tension (“~=0) (according to 


Tamm [9]); ¢ sound velocity, f/ frequency, 
d thickness of the “plate”. 


The result can be easily understood since the 
crossing waves form a standing wave pattern in the 
direction normal to the plane of symmetry, so that 
the boundary condition of vanishing pressure is 


d 


| 
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” fulfilled periodically in parallel planes, the distances 
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between which are given by the angle of incidence. 
Vice versa the distance between these planes which 
is given by the thickness of the plate, determines 
the discrete angles of incidence and therewith also 
the trace velocity (phase velocity of the “plate- 
wave’). For parallel incidence (corresponding to 
infinite plate thickness) the latter is equal to that of 


_ the longitudinal wave, while it is infinite for the 
_ case of normal incidence (corresponding to a thick- 


\ 


ness equal to half a wave length) this means dis- 
placements with equal phase along the direction of 
propagation. The course of the phase velocity be- 
tween these limiting cases is given by the hyperbola 
already mentioned above. The limit of infinite phase 
velocity determines a minimum thickness of the pla- 
tes or cut-off frequency respectively below which no 

_ progressive wave can occur but only displacements 
with equal phase decreasing exponentially with dis- 
tance. This “suppression”’-attenuation given thereby 
increases with decreasing thickness or decreasing 
frequency respectively. The wave form can be de- 
scribed by a superposition of two plane waves with 
imaginary values of #,,. Waves with phase velocities 
smaller than that of the longitudinal waves are not 
present. 

The equation for the phase velocity can be written 
in a form which is well suited for checking the ex- 
perimental results. If the limiting frequency f, and 
the corresponding wave number k, of the wave 
type n is introduced by the condition 


In 


nce 


Od 
v=c/V1— (fr/f)? 


or after some transformation 
Ae lic -(f? = f,3). 


This is the equation of a straight line with the 
slope 1/c?, if the square of the measuring frequency 
is plotted along the abscissa and the reciprocal of 
the square of the experimentally determined wave 
lengths along the ordinate. The relation also holds 
for propagation in rods with arbitrary cross-sections. 
The measurements made with rods of rubber which 
are described later, have been evaluated according 
to this method. 

In infinitely extended media capable of bearing 
shear forces two types of waves can occur which 
propagate independently from each other in the case 
of small amplitudes: the transverse wave (shear 
wave) with the phase velocity 


S|) 7H, €: 


we get 


cy=Vu/o (@ density of the medium) 


and the longitudinal wave with the phase velocity 


a= VOr+21)/0. 
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While in infinite media both types of waves occur 
independently, the boundary conditions at plane 
surfaces (vanishing stress- and shear-tensions) can 
be fulfilled only by superposition of both types: If 
a plane wave of one type hits the surface, in ad- 
dition to the reflected wave of the same type (angle 
of incidence equal to angle of reflection) a reflected 
wave of the other type is observed, the angle of 
reflection of which is given by the condition of 
equal phase velocity along the boundary surface. 
Besides the angle of reflection, amplitude and phase 
of the additional wave are also determined by the 
boundary conditions. In addition further hetero- 
geneous waves (Rayleigh’s surface waves) occur 
near the boundaries. 

If now the propagation in plates (of thickness d) 
is considered, the boundary conditions have to be 
fulfilled at two opposite boundary surfaces as in the 
case of the medium without shear tensions. Using 
the abbreviations 
q=(cr/e,)? and s= (cy/v)? =sin? B, 
this yields [6] the two equations 


for symmetrical modes: 


tan Ones pas ae 
2c _ _4sV1l—s Vq-s 
pe od Vee (1-25)? 
2 cy 
for antisymmetrical modes: 
od .; 
tan —— — S$ —— : 
Dey V4 _ _4sV1—sVq-s 
ces aod V1— (1—2:s)? 
2, CT 


where f is the angle between plate normal and 
direction of propagation (i.e. the angle of inci- 
dence) of the transverse waves. 

The roots of the equations thus determine the 
phase velocities 


v=cry/Vs=cry/sin B 


of the different kinds of plate waves, if the proper- 
ties of the material and the thickness of the plate 
are given. 

An evaluation of the equations can only be made 
numerically which has been done for example by 
Frrestone [1] for aluminium plates with o=0.35 
(see Fig. 2). The results have been stated by Scnocu 
[7] with measurements of the transmission of sound 
through plates. Since the nodal planes of the stress 
(or pressure) distribution of the two fundamental 
types forming the plate wave are no longer (as those 
of the “compressional” wave in the case of a liquid) 
situated on the boundary surfaces (for this purpose 
each type of wave ought to fulfil the boundary condi- 
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tions separately) , the thickness of the plate no longer 
corresponds to half a period of a sinusoidal pressure 
distribution or a multiple of it. 


b= C, 
V=Cy 
TN=p 
0 1 2 3 4 
fd/ey ——— 


Fig. 2. Phase velocity v of free waves in solid plates 
with a Poisson’s number o=0.35 (azcording to 
Firestone [1]); cr velocity of transverse waves, 
{ frequency, d thickness of the plate. 


The course of the phase velocity therewith de- 
viates from the hyperbola. Only for the case of nor- 
mal incidence can the boundary conditions be ful- 
filled for each type of wave separately, so that the 
cut-off frequencies of the longitudinal waves are 
unaltered and harmonic. They are completed by 
corresponding cut-off frequencies of the shear waves. 
As a whole we get the wellknown rather complicated 
picture of the dispersion curves which is shown in 
Fig. 2. It may be mentioned that here also waves 
with phase velocities can occur which are smaller 
than those of the longitudinal waves. These waves 
are superpositions of plane shear waves and inhomo- 
geneous longitudinal waves. The lower limit for 
the phase velocity of superpositions of this type is 
the velocity of the shear waves. A difference from 
the medium without shear tensions is given by the 
occurrence of extensional and flexural waves i.e. 
of waves without cut-off frequency, the phase velo- 
city of which depends only to a small extent on the 
dilatational stiffness 2*. It runs from a value some- 
what higher than the velocity of the shear waves or 
zero resp. to that of the Rayleigh waves; the latter 
being smaller than the velocity of the shear waves 
means that in this case the shear wave component, 
too, is inhomogeneous. 


3. Waves in rubber-elastic materials 


In the present paper the equations have been 
evaluated for the value o=0.47 for Poisson’s num- 
ber!. The resulting picture (see Fig. 3) is much 


1 The numerical calculations have been made by 
cand. phys. Gruser. 
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Fig.3. Phase velocity v of free waves in solid plates 
with a Poisson’s number of o=0.47; cr velocity 
of transverse waves, cy, velocity of longitudinal 
waves, cp velocity of Rayleigh’s surface waves. 


clearer than that for metals with a Poisson’s number 
of o=0.35. This is due to the fact that the lowest 
limiting frequency calculated from the shear-wave- 
velocity is much lower than that from the longitudi- 
nal-wave-velocity. The course of the phase velocity 
of the different types of waves may be understood 
as follows. Below the lowest “longitudinal-wave- 
limiting-frequency”’ the compressibility of the mate- 
rial is of low influence, so that the phase velocity 
is determined almost exclusively by the shear mo- 
dulus. The shear modulus on the other hand is 
so small that above this limiting frequency the 
waves excited scarcely deviate from the longitu- 
dinal waves in plates without shear tensions though, 
of course, the waves are superpositions of longitu- 
dinal and shear waves. Only for values of d// for 
which large amplitudes of the shear waves are 
necessary to fulfil the boundary conditions (thick- 
ness-resonance) do the waves behave with respect 
to the phase velocity like pure shear waves in an 
incompressible medium. This alternating influence 
of longitudinal and transversal components pro- 
duces the steplike course of the phase velocity. 

Above the limiting frequency of the longitudinal 
waves, a phase velocity largely fluctuating with 
frequency should be expected. Experience, however, 
shows, that the large number of dispersions is not 
found experimentally. This is obviously due to the 
attenuation which increases so much near the dis- 
persion frequencies that the corresponding wave 
types are damped out and cannot be measured. The 
increase of attenuation is caused by the higher 
energy losses of the shear component which here 
overcomes the longitudinal one. Apart from these 
small regions the attenuation is determined by the 
loss factor of the longitudinal waves. 
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While in plane waves a constant loss tangent 
means constant energy loss per wave length, the 
attenuation of the plate waves increases with phase 
velocity. The reason, of course, is that, correspond- 
ing to the reduced group velocity, the energy is pro- 
pagated more slowly, so that the energy loss per 
period, given by the loss tangent, is effective within 
a smaller distance so that the energy decreases more 
rapidly with distance. This process can be under- 
stood also from the picture of the superimposed 
obliquely propagating waves, if we take into con- 
sideration that the boundary conditions cannot be 
- fulfilled with plane damped waves, but that ad- 
ditional inhomogeneous waves of the kind of the 
surface waves are necessary which are described 
mathematically by complex angles of incidence. 
Physically this process also can be taken as the in- 
creased attenuation in the vicinity of a transverse 
resonance. In the case of the plate without shear 
_ tensions we therewith are able to determine the loss 
tangent in a simple way also from the values of the 
attenuation in regions of strong dispersion near the 
| limiting frequency. The influence of the loss tangent 
can be easily seen in those cases when the propaga- 
tion process is determined by one of the two wave 
types only. Phase velocity and attenuation of the 
extensional wave are given by Young’s modulus and 
its loss tangent. These two values are composed 
(s. Skuprzyk [8], p.771) in a simple way by 
Lamé’s constants (2* and sc) and their loss tangents. 
For soft plastics (60.5) Young’s modulus is 
mainly determined by the shear modulus («) and 
4 its loss tangent is practically equal to that of the 
shear modulus. 

If on the other hand both types of waves parti- 
cipate to more or less the same extent in the for- 
mation of the wave mode, the attenuation is de- 
termined in a rather complicated way by both the 
loss tangents. For the calculation of the 

attenuation in the most general case it is 
~ necessary to calculate the distribution of 
the potential energy of the compressions 
and the distortions separately and to de- 
termine the energy losses from the cor- 


} 


Worm drive 
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plates at one of their edges with uniform ampli- 
tude and phase along this edge. An excitation 
of this kind was realized by Kuni and Meyer ([3]. 
chap. III) by exciting a U-profile iron in its lowest 
cross-sectional resonance. In this resonance both 
sides of the cross section vibrate with opposite 
phase while the central part is exposed to bending 
moments, For the measurements, a strip of rubber 
which was so wide that it could be looked upon as 
an infinite plate was fixed with one of its edges at 
one side of the profile iron, the latter being excited 
to the mentioned resonance vibration. This method, 
however, only allowed measurements at single fre- 
quencies, while in the present work a very large fre- 
quency range (0.1 to 300 kc/s) should be investi- 
gated. 

Fortunately the experimental results with rods 
may be compared with the theoretical results for 
plates at least qualitatively, because it is known from 
calculations for rods of circular cross-section which 
can be and have been made exactly ([2], [5]) that 
here the dispersional behaviour of the different wave 
types qualitatively agrees with that in a plate. 

The apparatus which was used in the frequency 
range | to 300 ke/s is very similar to that used by 
Kuni and Meyer. The scheme is given in Fig. 4. 
The excitation of the rubber rods to vibrations in 
longitudinal direction was made from one end. For 
this purpose the rods were connected at their front 
surfaces with an aluminium block, having the same 
cross-section and about 2 cm length. The aluminium 
block was fixed on a transmitter for structure borne 
sound. In the range 3 to 300 ke/s a piezo-electric 
system was used as a transmitter, which mainly 
consisted of a block of Rochelle salt plates. In the 
frequency range below 3 kc/s the amplitudes of this 
system were too small to get a sufficient excitation, 
so that here a moving coil transmitter was used, the 


Synchronous motor 


_ responding loss tangents. This calcula- 
tion has not been made in detail up to 
now. 
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4. Measuring apparatus 


Contrary to the calculations which 
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could be made exactly only for the 


case of plates, the measurements have 
been made with rods of square cross- 
section, which is much easier than 
with plates. It is difficult to excite 


- 


Sound absorbing lining 


Fig. 4. Scheme of the apparatus for phase velocity and damping 
measurements, 
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coil of which bore a conical coupling mass of alu- 
minium. 

To reduce the radiation of air-borne sound, both 
systems had to be surrounded by thick-walled cap- 


sules of aluminium. The space between capsules and 


systems was filled with sound absorbing materials. 
The systems were supplied by an RC-generator with 
power amplifier with an output power of 80 watts 
for frequencies below 20kc/s and 15 watts above 
20 ke/s available. 

To pick up the vibrations of the rod it was plung- 
ed slowly into (or drawn out of) a water-filled 
vessel, 40 cm in depth, within which a hydrophone 
was arranged. If now the rod is excited to elastic 
vibrations from its upper end an elastic wave travels 
a varying distance along the rod where it is at- 
tenuated by the losses of the material, so that its 
vibration amplitude at the surface of the water is 
altered correspondingly. The part of the rod under 
water now radiates sound energy into the water by 
its contraction and at the place of the hydrophone 
a sound pressure is generated which is linearly com- 
posed of the sound waves coming directly from the 
rod and those reflected from the surface and the 
walls of the vessel. This sound pressure is therefore 
proportional to the vibration amplitude at any point 
of sound field e. g. the cross section of the rod at 
the surface of the water, if the sound distribution 
in the rod and in the water volume does not alter 
during the motion of the rod. This is fulfilled if the 
lower front side of the rod does not radiate sound. 
The attenuation in the rod therefore must be so 
large that no measurable sound energy reaches the 
lower end of it. Moreover the rod must not move 
in a transverse direction and the conditions of pro- 
pagation within the vessel must not be altered by 
the differences in the depth of immersion of the 
rod. For this purpose all walls of the vessel were 
lined with a sound absorbing layer of the “Fafnir’’- 
type (described by Meyer and Tamm [4]) so that 
the alteration of the sound absorption within the 
vessel caused by the varying length of the rod was 
negligible against the constant absorption of the 
lining. 

With the mentioned conditions fulfilled, this 
measuring method has the advantage when com- 
pared with other methods using vibration pick-ups 
brought into direct contact with the specimen at 
different places, that a completely uniform and re- 
producible contact between bar and microphone is 
secured. With vibration pick-ups in direct contact 
with the specimen, it is difficult, especially in the 
case of measurements at high frequencies and phase 
measurements, to get reproducible conditions. 

The measuring vessel filled with water was sur- 
rounded by another vessel filled with a liquid of 


=e 
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controlled temperature so that measurements could 
be made under constant conditions. For all meas- 
urements the water température: was kept constant 
at 20° C. Thus the specimen, too, had a temperature _ 
of 20° C, since it was taken off the water only for 
a short time and since its temperature altered only 
very slowly if the room temperature differed from 
20° C. Fig. 5 shows a photograph of the measuring 
vessel with the rod of rubber, the exciting system 
and the driving mechanism for the movements of 
the rod. 


Fig. 5. View of the measuring vessel. 


The output voltage of the hydrophone was ampli- 
fied by means of a resonance amplifier and a hetero-. 
dyne amplifier and recorded by means of a level 
recorder. A high selectivity of the receiving set is 
necessary to suppress the harmonics produced espe- 
cially by the piezo-electric transmitter, which are 
disturbing when the attenuation within the rubber 
rod is large for the basic frequency and small for 
its harmonics. Moreover the vertical amplifier of an 
oscilloscope is connected with the resonance ampli- 
fier, while the voltage of the RC-generator is given 
to the horizontal plates, thus allowing measurements 
of the phase shift between.transmitter and receiver 
voltage and measurements of the wave length on the 
bar when the latter is moved. The rubber bar to- 
gether with the transmitter is moved vertically (velo- 


a ey eel elie & 
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city 0.2 or 0.4mm/s) by means of a worm drive 
coupled to a synchronous motor. The decrease of 
the amplitude on the bar occurs exponentially so 
; that the level recorder records — on a logarithmic 
_ scale — a straight line from the slope of which and 


_ the attenuation (dB/cm) can be determined. 

As already mentioned, this procedure can no 
_ longer be applied when the attenuation within the 
~_ rods becomes so small that a considerable part of 
‘ the sound energy reaches its lower end and is re- 
_ flected there or radiated into the liquid. With the 
_ special sorts of rubber investigated we thus got a 
lower frequency limit of about 1 ke/s. To extend the 
| measuring range to frequencies below 1 kc/s some 
resonance measurements have been made in addition. 


_ two stretched steel wires was made without mechani- 
_ cal contact by means of an electromagnetic system 
ee acting on a small iron plate fixed at the front of the 
_ bar so that the latter was excited to vibrations in 
longitudinal direction. The vibrations were detected 
_ electrostatically by using the other end of the bar, 
_ which had been made conductive, as the moving 
electrode of a condenser microphone. The measure- 
ments were made by continuously altering the fre- 
quency of excitation and recording the amplitude 
_ with a level recorder, the resulting resonance curves 
being evaluated with respect to resonant frequency 
and half energy width. The lowest resonant frequency 
for extensional waves of the rod, 30 cm in length, 
- was 70 c/s in the case of the softest rubber type. 
Because of the relatively high loss tangent of the 
extensional waves used in this experiment (7~0.1) 
_ the additional losses caused by the support of the 
rod and the radiation of airborne sound do not 
need consideration. 


*. 


5. Experimental results 


The final experiments were made with five diffe- 
rent types of rubber”. Of these types, two were 
__ based on synthetic and three on natural rubber each 
type having a different degree of hardness which 
has been given in terms of DVM-hardness*. The 


2 The ‘samples were placed at our disposal by the 
firm Continental-Werke A.G., Hannover, which we grate- 
_ fully acknowledge. 


__ 8 DVM-hardness is a measure which is used in ad- 
_ dition to the Shore-hardness to characterize the elasticity 
_ of rubber and plastics. It is determined by means of a 
metal sphere, 10 mm in diameter, which is pressed on 
to a plate of 6 mm thickness of the material investigated 
on a rigid foundation, first with a weight of 50 ¢ and 
then with an additional weight of 1000 g. The difference 
between the depths of penetration measured in hund- 
- redths of a millimeter gives the DVM-number. Thus 


a the softer material has the larger DVM-number. 
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types of synthetic rubber (Neoprene) investigated 
had the DVM-numbers 25 and 65 and the natura! 
rubber the numbers 25, 65 and 95. Of each type 
one rod of 5 x5 and another of 10 x 10 mm? cross- 
section were measured, their length being about 
30 cm. Typical examples for the experimental results 
with respect to attenuation and phase velocity are 
given in Fig. 6. Figs.6a and 6b show the results 
found with a hard type of rubber, Figs. 6c and 6d 
with a soft type. 

We found out that for all types of rubber in- 
vestigated, the whole frequency range of the experi- 
ments from 0.1 to 300 ke/s can be subdivided into 
four parts which differ from each other with respect 
to the courses of phase velocity and attenuation. The 
limits of these ranges depend on the type of rubber 
and on the cross-section of the rods. 


1. The lowest part of the frequency range is 
characterized by attenuation and phase velocity be- 
ing equal for both rods of different cross-section of 
each type of rubber within the limits of measuring 
accuracy. The attenuation (measured in dB/cm) in- 
creases linearly of somewhat more than linearly 
with frequency. The phase velocity remains constant 
or shows a slight increase with frequency. We here 
have ,,genuine” extensional waves, the wave length 
being at least 6 times as large as the diameter of 
the rod. Within this frequency range we can de- 
termine from the experimental values of attenuation 
and wave length a loss tangent which is characte- 
ristic for the material. We also can use the experi- 
mental values of the phase velocity to determine 
Young’s modulus. 

Since the loss tangent is comparatively large, the 
usual approximation formulas cannot be used for 
calculating it from the attenuation, but an exact. 
relation must be applied derived from the solution 
of the wave equation for a damped wave propa- 
gating in the x-direction in a medium with a com- 
plex modulus (called “Sound Stiffness”) D(1 +i) 


A(x, t) = Ay e 1@-192 eiot — 
= Ay e7ie VolD(1+in) z eiot 


For arbitrary values of 7 we then get the attenuation 
per wave length 


piea gn Oa) a eee oath at 

U7] dB Yt 4 
This relation holds for all types of waves the phase 
velocity of which is proportional to the square root 
of an elastic modulus, e. g. for extensional waves 
and for longitudinal waves, but not for the calcula- 
tion of the loss tangent from the attenuation of 
flexural waves. The relation given by this formula 
is represented in Fig. 7. For small loss tangents 
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Fig. 6. Phase velocity and damping of elastic waves in bars of two types of natural rubber for two different cross- 
sections; @ @ @ phase velocity //, 000 damping a, 
(a) 55 mm?, hardness 25 DVM, (b) 10x 10 mm?, hardness 25 DVM, 
(c) 5X5 mm?, hardness 95 DVM, (d) 10X10 mm’, hardness 95 DVM. 


(7? <1) we get the approximation 10 ] 
aa 
BI le 08 | 
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With large loss tangents also the simple relation cae 
between frequency, wave length and phase velocity 
c=f2/ is no longer valid. From the equation for a nag 
damped wave propagating in the x-direction we get 26 
with 2 defined by e'” =e?" : 
pea Pi, a 10 
1+?) 1 
This relation too is a ini Big suit ihe 
2. For higher frequencies, i.e. smaller wave & Le x SOdBkmom 


CA, 
lengths, we may still use the term extensional waves, 


Fig. 7. Loss factor 7 and ratio c/f A as function of the 
but here the geometry of the rod has a large effect 


attenuation per waye length. 
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on the measured attenuation of propagation and the 
latter therefore no longer depends exclusively on 
the material of the rod. Here the calculation of a 
loss tangent is fruitless and therefore the product 
a4 was plotted as a function of frequency. Some 
experimental results for different types of rubber 
and for two different cross-sections are given in 
Fig. 8. The plotting extends also over the frequency 
range (1) showing there a good agreement of the 


~aA-values for both cross-sections. At a ratio of 


A/d~5 (for the larger cross-section), the experi- 
mental results for rods of different cross-sections 
deviate from each other. Both of them show linear 
increases but with different slopes, the rod with the 
larger cross-section having the steeper one. The 
experimental results have been plotted as far as a 
finite wave length could be determined, the know- 
ledge of which is necessary for the calculation of « /. 
The corresponding frequency gives the upper limit 
of the frequency range (2). | 
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Fig. 8. Damping in bars of soft rubber of two different 
cross-sections. 
(a) Natural rubber, hardness 65 DVM, 
(b) Natural rubber, hardness 95 DVM, 
(c) Neoprene, hardness 65 DVM. 
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As can be seen from a comparison between the 
magnitudes of the pure material damping in the 
frequency range (1) and the damping additionally 
influenced by geometry in the frequency range (2), 
the latter is considerably higher. The linear increase 
of a/ with frequency seems worth mentioning and 
it would be of interest to find the reasons for this 
by means of theoretical considerations on the me- 
chanism of the propagation of extensional waves in 
that frequency range within which the wave length 
is comparable with the diameter of the rod. Up to 
now, however, such considerations have not been 
made with success. 

As well as the lower limit, the upper limit of this 
frequency range (2) not only depends on the cross- 
section but also on the hardness of the rubber rod. 
While for the rod 10 x 10 mm? in cross-section with 
hardness 95 DVM finite phase velocity only can be 
measured up to 7kc/s, finite phase velocities are 
observed in a rod with the hardness 25 DVM and 
the same cross-section up to 30 kc/s. This, of course, 
is a consequence of its larger phase velocity, i.e. 
longer wave length at the same frequency. The re- 
sulting ratio of the limiting frequencies (30:7) is 
not the same as the ratio of the phase velocities 
(300 : 25) because of the fact that the loss tangent 
of the harder rubber type is higher than that of the 
softer type. Therewith the constant “geometrically 
induced” attenuation of the “compressional wave” 
of the next frequency range (3) is only reached for 
higher values of d/A. This type of wave then pre- 
dominates because of its lower attenuation. 

3. The frequency range (3) is characterized by 
the infinite phase velocity of the “compressional 
wave” which means that the rod vibrates in equal 
phase at any point. This is obvious from a com- 
parison between Fig.9a and Fig.9b, which show 
the phase of the vibration along a distance of nearly 
lcm. Fig.9a shows the conditions for the exten- 
sional wave in the frequency range (2). In the 
immediate neighbourhood of the upper end of the 
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Phase shift between transmitter and receiver 
voltage dependent on the transmitting length of 
the rod of rubber; 

(a) frequency range (2), 

(b) frequency range (3). 


Fig. 9. 
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rod the transverse contraction is hindered because 
there the rod is connected to an aluminium block. 
Thus the phase velocity is increased here, but be- 


ginning at a distance of some millimetres, a con- ° 


stant value is observed which has been given as the 
result. Fig. 9b holds for the “compressional wave” 
in the frequency range (3). There is no phase al- 
teration along the small distance within which 
— because of the high attenuation — the vibration 
can be observed at all. 

In the frequency range (3) the attenuation (meas- 
ured in dB/cm) is nearly constant. Near the upper 
limit of this range, however, it decreases first slowly 
and then very steeply, the steep decrease being cha- 
racteristic for the upper frequency limit of this fre- 
quency range. The same course is observed for the 
wattless attenuation when a liquid column with soft 
boundaries is excited to vibrations below its limiting 
frequency. This limiting frequency of a liquid co- 
lumn with square cross-section and without losses 
in the liquid and at the boundaries is given ([3], 


chap. VIII) by 


her a y2 5: 


with d being the edge length of the cross-section 
and c the sound velocity in the free medium. Below 
the limiting frequency the value of attenuation is 
given by the formula 


= 
(ofa) : 


Co 


and approximates to the constant value a at very 
low frequencies: 


This value does not depend on the properties of the 
liquid but only on the cross-section of the liquid 
column. If this limiting value ap is calculated for 
the rods of rubber as if they were columns of liquid 
we get 

at 10 x 10 mm? cross-section: 38.5 dB/em and 

at 5 5mm? cross-section: 77 dB/cm. 
A comparison with the experimental results for 
those soft types of rubber for which the constant 
attenuation could be measured over a larger fre- 
quency range, shows that within this range the 
elastic properties of the rubber can very well be put 
equal to those of a liquid. The somewhat higher 
attenuation observed is at least partly due to the 
losses of the material. The experimental values found 
by Kuni and Tamm ([3], chap. VIII) for water 


columns with soft boundaries also are somewhat 
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higher than the calculated values in the range of the 
constant attenuation, which may be due to energy 
losses at the “sound-soft” walls. 

As was stated in the last section, the lower limit 
of the frequency range (3) depends on the softness 
of the rubber for rods with a given cross-section. 
The upper limit on the other hand is given by the 
velocity of the longitudinal waves which is nearly 
independent of the softness of the material. We 
therefore get a frequency range (3) with infinite 
phase velocity and nearly constant attenuation which 
is the larger the softer the type of rubber. 

4, In the frequency range (4) above the limiting 
frequency of the compressional wave once more a 
finite phase velocity can be measured and we find 
a dispersion curve which corresponds to that in a 
liquid and which begins with very large values of 
the phase velocity at the limiting frequency. This 
limiting frequency is about 210 ke/s for rods with 
5 x5 mm? cross-section and about 105 ke/s for the 
rods with 10 x 10 mm? cross-section, i. e. inversely 
proportional to the linear dimensions of the cross- 
section as is to be expected for a liquid. With in- 
creasing frequency the phase velocity approaches 
asymptotically the phase velocity of longitudinal 
waves in an infinitely extended medium, which 
equals about 1500 m/s for all types of rubber in- 
vestigated. 

Systematic deviations of the dispersion curves 
from their hyperbolic forms have not been found. 
A plotting of the experimental results in the above 
described way (1/2? as a function of /?) shows, at 
least for the softer types of rubber, also quantita- 
tively a good agreement of the phase velocity with 
that in water (see Figs 10a and 10b). The inter- 
sections of the straight lines with ordinate and ab- 
scissa give — 1/d? and f,?, their slope 1/c? (e velo- 
city of the longitudinal wave). 

The change of the attenuation connected with the 
dispersion is very steep and equals a factor of 100 
for change of frequency of only some ke/s — de- 
pending on the special type of rubber. It does not 
seem impossible to use this large alteration in at- 
tenuation for the construction of acoustic filters. 

In spite of this large alteration of the value of 
the attenuation, a nearly constant value for the loss 
tangent is obtained if the attenuation — measured 
in dB/em — is multiplied by the quotient c.,/f and 
the loss tangent is calculated from this product (as 
from the product a 4 in section 3.1). By this method, 
the fact that the vibration energy travels with the 
group velocity is considered. 

The results of the measurements have been plot- 
ted in Figs. 11 and 12 resp. Fig. 11 shows the loss 
tangent of the two types of synthetic rubber. The 
values derived from the product «4 are plotted for 


x 


_ Fig. 10. Dispersion of 
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Fig.12. Loss tangent measured with “vompressional 


waves” in bars of natural rubber of different 
hardness. The values are calculated from a cgr/f; 
© 10x10 mm?, 

® 5 x5 mm? cross-section. 


comparison. Fig. 12 shows the loss tangent of the 
three types of natural rubber. It is larger the har- 
der the rubber and equals between 0.01 and 0.05, 
i.e. it is considerably smaller than the loss tangent 
of the extensional waves. 

Fig. 13 shows a comparison between the theore- 


tical course of the attenuation in the whole fre- 


quency range and the experimental results got with 
a rubber rod with the hardness 65 DVM and a cross- 
section of 10 x 10 mm?. The full drawn straight line 
at low frequencies has been calculated for pure exten- 
sional waves with the constant loss tangent 7p = 0.44 
and the constant phase velocity cp=85 m/s. The 
calculated course of the blocking attenuation in a 
liquid column of the same dimensions is given by 
the cross-section of the rod and the measured velo- 


0? kk 10° 


Fig. 13. Comparison between the theoretical course of 
the attenuation and the experimental results 
in bars of rubber (neoprene 10 X 10 mm”, hard- 
ness 65 DVM). 
ap calculated for extensional waves 
(cp =85 m/s, np =0.44), 

ay calculated for longitudinal waves 
(c.=1590 m/s, mL=2...3% 10-2), 

ap calculated for a liquid column below limit- 
ing frequency (cL =1590 m/s, 7=0). 
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city for longitudinal waves in the infinitely extended 
medium, being 1590 m/s. For the attenuation of 
longitudinal waves calculated from the group velocity, 


a linearly increasing loss tangent 7, =2to3 x10? ° 


has been assumed. 

The experimental results can be summarized by 
stating that for an excitation in the longitudinal 
direction the sound propagation in rods of rubber 
or plastics is determined by either extensional or 
compressional waves. At low frequencies only the 
extensional waves can be excited, so that phase 
velocity and attenuation are given by Young’s mo- 
dulus (in the case of rubber nearly exclusively de- 
termined by the shear modulus) and its loss tan- 
gent. 

For high frequencies the attenuation of the 
pressional waves” is much smaller than that of the 


ce 


com- 


extensional waves so that only the first ones are 
observed, the sound propagation being determined 
by the bulk modulus and its loss tangent. In this 
frequency range the rods behave like columns of 
liquid. 

The ranges of transmission of these two types of 
waves are separated by a more or less wide {fre- 
quency range of low transmission of both the waves. 
In this frequency range the extensional wave shows 
an attenuation considerably increased by the geo- 
metric influence of the cross-section while the pro- 
pagation of the compressional wave is blocked, be- 
cause it is excited below its limiting frequency. The 
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limits of this intermediate range are given by the 
cross-sectional dimensions. of the rods, which are 
comparable with the wave length of the shear waves 
at the lower and that of the longitudinal waves at 
the upper limit. 


The work described in this report was carried 
out under arrangements with the Department of 
Scientific and Industrial Research, London, whose as- 
sistance and interest are gratefully acknowledged. It 
is published with the permission of the Department. 
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A NOTE. ON RECIFROGIVY 
IN LINEAR PASSIVE ACOUSTICAL SYSTEMS 
by J. H. Janssen 


Technical Physics Department T. N. O. and T. H., Delft, Netherlands 


Summary 


It is shown theoretically that linear acoustical passive systems can behave as so called 
reciprocity-violating systems; an example of such a system is a sound absorbing material 


that is porous as well as flexible. 


Sommaire 


On montre par la théorie que les systémes acoustiques passifs et linéaires peuvent ne pas 
satisfaire au théoréme de réciprocité. On en donne un exemple constitué par des matériaux 


absorbants a la fois poreux et flexibles. 


Zusammenfassung 


Das Reziprozitatsgesetz braucht nicht immer erfiillt zu sein fiir passive lineare akustische 
Systeme, wie zum Beaispiel porése Materialien mit elastischem Skelett. 


1. Introduction 


It is common engineering practice to calibrate 
electro-acoustic transducers by the so called reci- 
procity method (e.g. [1]). 


In doing so the transducers are supposed not to 
violate the principle of reciprocity. If need be the 
same assumption is made as to the surrounding 
medium. Reciprocity is also — be it tacitly — as- 


3 De eet, taal? 
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sumed in the case of the measurement of the airborne 
sound insulation in buildings. 

With this in view it seems worthwhile to draw 
the attention to two cases, that one meets even 
rather frequently in practice, in which a set-up is 
used that does violate the reciprocity principle. One 
case is that of flexural waves in plates [2]. Another 
case that will be dealt with in this paper is that of 


_ flexible porous materials. 


taciaieeeeeedlieelll 


in | 


d 


As yet we cannot overlook the practical conse- 
quences of the possibly non-reciprocal behaviour of 
such materials, but it seems worthwhile to study 
this case, since flexible porous materials are used 
a great deal in acoustics, e. g. in anechoic rooms. 


2. Reciprocity 


Let us consider a closed vessel, the walls of which 
ere rigid, containing air (Fig. 1). In the vessel there 
are two sound sources A and B and a cylindrical 
rigid, inert tube with open ends. The sound pres- 
sures Py and p; as well as the particle velocities vg 
and v; are for simplicity assumed to be constant all 
over the open ends. They may vary with time how- 
ever. This simplification does not contain a loss of 
generality. 


Fig. 1. An idealized linear acoustical fourpole, in the 
form of an inert tube, inside a closed vessel con- 
taining air. 


The sound field ~4(yg) that occurs when the 
sound source A(B) is in action satisfies (see list 
of symbols at the end of the paper) : 


V? path’ pa=—44; 
V? va t+? op = — 4B. 
Outside the region A (B) the specific source strength 


44 (qp) vanishes. 
According to GREEN: 


/ (~4 Oyp/dn — pp Oe-4/On) do = 


= if (g4PB— 4B Pa) At, 
where o and t stand for surface and volume respec- 
tively. The surface integral vanishes when Oy/On 
is zero everywhere on the surface (rigid surface). 


It is already sufficient, however, when the surface 
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is locally reacting, i.e. when the ratio of y and 
S/n only depends upon the properties of the 
surface at the point where these quantities are con- 
sidered. 


We then have 


pa=f(o) dv4/On, 
pp=f(o) Opz/On, 


which obviously results in a zero value of the sur- 
face integral. 

If the regions A and B are very small in com- 
parison to the wave length we may write 


px [ qadt—ga | qpde=0, 


which yields the well known reciprocity theorem in 
its simplest form 


paVa=ppVa 
or since @ is proportional to the sound pressure 


paVe=psV4. 

Reciprocity is apparently due to the vanishing 
of the surface integral which condition may be re- 
written in p and v, since these quantities are pro- 
portional to m and O@/On respectively: 


| (P42my —PB em ,) d=. (1) 


We are allowed to exclude the volume occupied by 
the cylinder in Fig. 1 from the total volume t in 
question. Reciprocity between A and B then how- 
ever requires that the surface integral (1) vanishes 
over the outer surface of the cylinder, its open ends 
included, The outer cylindrical surface of the cylin- 
der does not contribute to the integral since the 
cylinder is supposed not to move. It follows that 


(2) 


is the condition for reciprocity, equivalent to (1). 
This condition is satisfied for a tube filled with air. 
If, however, the tube is filled with an absorbing 
material conditions might be changed. 

Suppose py=api+bv, and vy=cp,+duy, and 
ad—bc=A., Then one has for this fourpole 


(P04 VoB — Pos VoA) — (P14 Vip — PiB ViA) =9 


(po4 Vos — Pos Vos) — A(pra Vip — Pip Via) =O. 


One sees at a glance that only for J=1 reciprocity 
will hold, because only then (2) is satisfied. 


3. Elastic porous materials 


The theory of the absorption of sound in porous 
and elastic materials is fully treated in [3], [4], 
[5]. The equations of motion for the material of 
the skeleton and for the air inside the pores are 
given by 
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01 Ov,/Ot+ 5 (vy — vg) — (1—h) (Ko/j w) 0? v/d2? — 
—[K,+ (1—h)?(Ky— Py) /h] (1/j @) 0?0,/da? = 


(Cra)on 


Or v9/Ot +s (vg—04) — (h Ks/j w) 9?v9/Bx? — 


— (1—A) (Ky — Pp) (1/j @) 0v,/dx? =0 (3b) 


These equations may be treated as the Lagrangian 
form of the description for the motion of continuous 
mechanical systems (see for example [6]). RayLeicH 
proved the existence of reciprocity [7], unless so 
called gyrostatic terms occur. These terms are found 
in (3a) from the coefficients for 02v,/Oz? and in 
(3b) from the coefficients for 0?v,/d2z respectively 
(taking 0/dx = — 7): 


[—s/y?— (1—h) K,/j] 
[ —s/y? — (1—A) (Kg —Py)/j ]. 


It will be seen by comparison with [7] that Ray- 
tEIGH’s f,s is formed for this case by (1 —h) Py/2 j w. 

Already in 1951 it was found theoretically by the 
present author that in deducing the fourpole coeffi- 
cients for layers of elastic and porous materials 
results were obtained which violated the reciprocity 
principle, unless h—> 1. This fact may now be ex- 
plained. The differential equations describing the 
physical phenomena in such materials, as derived 
in [3], are of such a form as to result in fourpole 
coefficients a, b, c and d, which give 41 unless 
the layer or blanket is symmetrically built or h— 1. 
It may be expected therefore that for low frequen- 
cies and small values of h (the porosity) non-reci- 
procal fourpoles in the form of flexible blankets 
experimentally can be demonstrated; for example 
when one side of the blanket is open to the air and 
the other side is closed by a very thin (e. g. mass- 
less) membrane or film which is fixed on the flexible 
skeleton of the material. 


and 
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List of symbols used in this paper 


~A,9VB_ Velocity potential for the sound field originat- 
ing from source A, resp. B, 

qdA>QB Specific source strength of source A, resp. B, 

n normal to boundary of vessel, 

k wave number (=/c), 

Oo boundary surface of vessel, 

G volume of vessel, 

V4,Vz_ source strength of source A, resp. B, 

On mean velocity of the solid material at a, 

01 density of the solid material as given by the 
mass of solid material per unit volume of 
porous material, 

s coupling coefficient, 

Ve mean velocity of the air at 2; this quantity 
is 1/h times greater than the velocity of spe- 
cific volume displacement, 

h cavity or porosity factor, i.e. ratio of the 
volume of the accessible holes to the total 
volume of the porous material, 

Ky normal modulus of air, if necessary to be 
taken complex in the case of losses being 
present (thermal or viscous), 

K, specific stiffness of the frame, to be taken 
complex if necessary, 

1 total equilibrium pressure of the air, 

Oo density of the air as given by the mass of the 
air content per unit volume of porous mate- 
rial; it is h times the density of free air, 

? =j oe, 

c complex sound yelocity, 

a) pulsatance. 
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ON IaE THEORY OF THE RAYLEIGH-DISK 
AND THE SOUND PRESSURE RADIOMETER 


by H. Hoseaarp Jensen and K. Sarrmark 


Physics Department, Technical University of Denmark, Copenhagen 


Summary 


Expressions are derived for the mean torque H and the mean force K on a thin, rigid 
obstacle, placed in a plane, progressive sound wave of arbitrary wavelength. Numerical cal- 
culations are carried out for a long strip and a circular disk, in the latter case based upon 
earlier numerical results by Koran1, Bouwxame and Kawat. It is found that, when the wave- 
length is varied, both H and K go through maxima and minima owing to resonance effects 
in the diffraction by the strip or disk. 


Sommaire 

On établit des expressions donnant le moment de torsion moyen H et la force moyenne K 
sur un obstacle mince. et rigide qui se trouve dans une onde sonore plane et progressive, 
de longueur d’onde arbitraire. On fait le calcul numérique de H et K pour une bande longue 
et un disque circulaire; ce dernier cas est traité a partir des résultats numériques antérieurs 
de Korant, Bouwxamp, et Kawai. On trouve que lorsque la longueur d’onde varie, H et K 
passent tous deux par des maxima et des minima, a cause d’effets de résonance intervenant 
dans la diffraction par la bande ou le disque. 


Zusammenfassung 

Es werden Ausdriicke gegeben fiir die mittlere Kraft K und das mittlere Drehmoment H 
auf einer diinnen Platte, die von einer fortschreitenden, ebenen akustischen Welle will- 
kirlicher Wellenlange getroffen wird. Numerische Berechnungen werden fiir das lang- 
gestreckte Band und die runde Scheibe durchgefihrt. Im letzten Falle sind die Berechnun- 
gen auf friihere numerische Resultate yon Koran1, Bouwkamp und Kawai aufgebaut. Wenn 
die Wellenlange variiert wird, durchlaufen sowohl H als K Maxima und Minima, ver- 


ursacht durch Resonanzeffekte bei der Beugung an der Platte. 


1. Introduction 


A sound wave impinging on a rigid body will 
exert a mean pressure, the radiation pressure, on the 
surface of the body. Consequently the resultant tor- 
que and force on the body usually will have non- 


zero mean values H and K. This is utilized for 
measurement of sound intensities in two kinds of 
apparatus. 

In the first, the Rayleigh-disk, the mean torque H 
on a plane disk or strip is measured. When the dia- 
meter of the disk or the width of the strip is so much 
less than the wavelength of the sound that the flow 
around the body is unaffected by compressibility, 
the torque will be independent of wavelength and is 
known in terms of the sound intensity from the 
theory of Kénic [1] and Cisorr: [2]. When the 
wavelength decreases and becomes comparable to 
the dimensions of the disk, the torque will vary and 
ultimately vanish, The Rayleigh-disk, therefore, is 
primarily an apparatus for absolute measurement 
of low-frequency sound. In the other apparatus, the 


_sound-radiometer, the mean force K on a plane disk 


is measured. The theory of this instrument is only 


simple when the wavelength is so much smaller than 
the disk that the shadow of the disk is sharp. In this 
case K is independent of wavelength and related to 
the sound intensity in a simple and well known way. 
It follows that the sound-radiometer is especially 
suited for measuring high frequency sound. When 
the wavelength increases, the force will vary and 
ultimately vanish in the limit of very long wave- 
lengths. 

However it may sometimes be necessary to use 
the instruments in regions of wavelength where the 
above mentioned conditions for the validity of the 
simple theories are not fulfilled. It is the aim of this 
note to discuss the way in which H and K depend on 
wavelength in such cases. We treat the two-dimensio- 
nal case of the strip and —in less detail — the three- 
dimensional case of the circular disk. 

Parts of the problem have to some extent been 
treated in the literature. Koran1 [3] has given an 
analytical expression for the torque on a circular 
disk and a series expansion in powers of ka, where 
k is the wave number and a the radius of the disk. 
This expansion is valid only for large wavelength. 
Awaranr [3] and Kawar [3] have treated the circu- 
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lar disk and calculated H and K numerically for 
several values of the wavelength and the angle of 
incidence. Furthermore, Levine [3] has obtained 


an asymptotic formula, valid for very short wave-— 


lengths, for the torque on a strip, while Ketter (Re- 

search report from the institute of Mathematical 

Sciences of New York University) has shown, how 

improved asymptotic formulae for H and K may be 

obtained. His methods are applicable down to wave- 
lengths as large as the scattering obstacle, and may 
thus be used for all ka-values not dealt with in this 
paper. 

We make the following assumptions: 

1. The disk or strip is freely suspended in a practi- 
cally unlimited radiation field. 

2. It can be treated as rigid. This requires that it 
has a mass much larger than the equivalent acou- 
stic mass of a radiating body of the same form 
and size. This condition can usually be fulfilled 
in gases, but not in liquids. The effect of the 
forced vibrations and rotations of the circular 
Rayleigh-disk has been discussed theoretically by 
Kine [4], Woop [5] and Kawar [3]. 

3. Sound absorption, vortex movement and turbu- 
lence are assumed negligible. We therefore treat 
the radiation field as a potential field, and the 
suspended body as non-absorbing. 

4. Though some of the results in the following are 
valid for a wider class of radiation fields we 


usually consider only the case of the incoming » 


wave being plane, progressive and monochro- 
matic. 


2. General expression for the torque 


We consider a periodic undisturbed radiation 
field with velocity potential ~). The plane disk or 
strip P (Fig. 1) is placed in this field. P will pro- 


Fig. 1. Geometry and nomenclature for the scattering 
problem. 


duce a scattered wave y. The resulting velocity po- 
tential is p=@)+y.- If my and w are known (see 
section 4) the torque on P around the vertical z-axis 
can be calculated in the following way: 
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Consider a closed surface F surrounding P. Angu- 
lar momentum is transported into the interior of F 
through action of the pressure from outside and by 
flow of matter carrying angular momentum. It is 
removed through P by the torque acting on P. 
Owing to the periodicity of the radiation field, the 
mean value of the angular momentum of the medium 
inside F is constant. Consequently the mean rate of 
transport from outside F must equal the mean rate 
of removal of angular momentum through P, but the 
latter is equal to the mean torque H. Furthermore, 
if we let F be a surface of revolution with axis co- 
inciding with the z-axis} the torque due to the pres- 
sure from outside F vanishes, and H is then equal 
to the mean rate of transport of angular momentum 
into F through material flow. Denoting the particle 
velocity by wu, the density of the medium by @ and 
using the notation of Fig. 1, this can be expressed 
by the equation 


H= — [r Quy uz df~ — Qo [run u df, 


< 


(1) 
EF 
where a bar denotes mean value over a whole num- 
ber of periods. In eq. (1) we have put @ equal to 
0 , the undisturbed density, as we shall only be 
interested in calculating H to second order in the 
particle velocity uw. 

We now express u in terms of the velocity poten- 
tial 


u=grad¢p, (2) 


and at thesame time we introduce the usual complex 
notation with the time-factor e '®’. We then get 


2 1 
Hl eae 0 | rgrad, p gradyy* df +conj. (3) 
F 


This expression is valid for arbitrary bodies and 
periodic radiation fields. Further simplifications are 
possible, when the scattering body is plane and of a 
thickness in the y-direction, small compared to the 
wavelength. We first consider the twodimensional 
case of P being a strip. F is then most conveniently 
taken to be a cylinder, and eq. (3) leads to the fol- 
lowing expression 


ai 1 ) x : 
H=— 400 | (P+) ao (Pot) df+conj , 
F 


4 
where @ is defined in Fig. 1. on 
The integral in eq. (4) is evaluated by means of 
the asymptotic value of the field for large values of 
r. In this case we have, because y represents an 
outgoing progressive wave, that 


Oy/Or=iky, (5) 


where k& is the wave number. 
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We also note that % depends on @ only through 


the combination 0 —v, where s sl is the angle 


of incidence of the incoming wave. Consequently 


Op 9/90 = — Op4/Sv . (6) 

Finally we use the fact, that 
y(0) = —y(—8). (7) 
This follows, because the scatterer is plane, thin and 


rigid. 

Now eqs. (5), (6) and (7) are introduced into 
eq. (4). Owing to eq. (7) only terms containing 
both the incoming and scattered waves contribute 
to the integral. The terms with 0/30 and dw*/d0 
are simplified by partial integration over 0. We get 


= 1 " ; : 
H-- 2 a| |v = fa +ikgy}| dj +con : 
re 3 (8) 


This expression is very convenient, because the ex- 
pression between curly brackets when resolved in 
partial waves will have the same asymptotic be- 
haviour for roo as y. 

In the threedimensional case F is taken to be a 
sphere. By the same methods as above it is then seen 
that eq. (8) still holds. Only r now means the radius 
of the sphere. 


3. General expression for the force 


In this section we shall show, that there is a close 
connection between the mean force K on a thin disk 
or strip, and the scattering cross-section Q; in fact 

K =o) u,2sinv'Q, (9) 
where wy? is the mean square of the particle velocity 
in the incoming plane progressive wave. WESTERVELT 
[6] has derived a general expression for the radia- 
tion pressure force on an obstacle of arbitrary shape. 
In the special case of a plane scatterer it can be re- 
duced to eq. (9). For the sake of completeness we 
give a direct proof of eq. (9). We start by investigat- 
ing the mean value Ap of the sound pressure at the 
surface of the scatterer, 

In general Ap is determined by the following ex- 
pression, which is correct to second order in wu: 


Ap== alee), (10) 
where c is the velocity of sound. 
In complex notation eq. (10) is written 
Ap = 4 oo 2 py" — grad p grad mp"). (11) 


From Ap we can calculate K by integrating over the 
entire surface of the scatterer 
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K= | Apndf, (12) 

3 ib 
where n is a unit vector in the direction of the inner 
normal to the surface of the scatterer. 

The expression to be used for @ in eq. (11) is 
(compare Fig. 1) 

P=%+y= a CASS Ia EEN IEE (13) 
where u,, is the velocity amplitude of the incoming 
wave. 

Combining eqs. (13), (11) and (12) and remem- 
bering, that when we pass from the front to the rear 
of the plane scatterer, ~ is unaltered while w 
changes sign, we get 


ian 1 2 * <) <) ) 
<< 00 (Peov" ae) ae “Jn df-+conj (14) 


The last term in eq. (14) is integrated by parts; y 
is zero at the edges of the disk, so the integrated 
parts vanish. The second derivative of Gp is obtain- 
ed from eq. (13), whence 


. ia BO tin . 
K=- ov | (Reuy + SEP wt) nal +-conj = 


(15) 


By means of eq. (13) and the boundary condition, 
expressing that the normal component of the velocity 
is zero at the surface of the disk, we can now elimi- 
nate P,. In fact we get 


eo is OP Oy 
Boni: Gyeon 
whence by means of eq. (15) 
ayiksin { (y" = —y oy we) af (17) 


It 


y= (16) 


+iksinv pp 


ae 
4 


where Ow/On is equal to + Oy/dy for y> ¥0. 
ey The: expression (10) is well known from the lite- 


rature (see e.g. [4]). 
It is easily derived from Bernoulli's equation for a 
compressible medium 


Pp 
: 1 K 
P+ s+ | 0, 
Po 


where py ist the pressure far from the radiation field. 
We expand 1/o in powers of dp=p—po, breaking off 
the series at the second member. Remembering that 
Oo0/dp =1/c? we get 

1 2 

52 (4p)? =0 


PAY Olt oes 


ay tee 
p+ —u?+ —Ap— 


2 Qo 


From this eq. (10) 
because p=0 and Ap? = 


is obtained by taking mean values, 
(09 g)® to secoad order. 


82 H. HOJGAARD JENSEN and K.SAERMARK: THEORY OF THE RAYLEIGH-DISK 


The expression eq. (17) is now compared with 
the power scattered by the disk. If the disk were 
vibrating in such a way, that the radiated wavefield 


was y, the mean power transmitted to the medium ° 


would be 


— [Apundi ron [Pay (8) 
P P n 
where the expression Ap = — 0 p for the first order 
pressure has been used and only second order terms 
in W are retained. 
In complex notation eq. (18) is 


W= 0 Bae] y: SY af 4 conj . (19) 


Using the relation w=ck we get from eqs. (17) 
and (19) 


==. SIMU 
K= ——F, 
Cc 


(20) 


but the cross-section Q is defined as the ratio bet- 
ween the scattered power W and the incoming sound 


intensity Q9 C ug”. Therefore eq. (20) proves eq. (9). 


The actual calculation of K may again be made 
by means of the field far from the scatterer. Here 
the scattered wave is very nearly plane, and conse- 
quently the intensity is by eq. (5): 


1 * 
ee G hey (21) 


W is then obtained by integrating eq. (21) over a 
closed surface surrounding the scatterer and every- 


where lying in the distant field. From W and 
eq. (20) we then get 


1 
5% pene hae (r+). 


Ke (22) 
However, in case of the circular disk and the strip, 
Q is known from the literature and K is then obtain- 
ed directly from eq. (9). 


4. Calculation of the scattered wave wy 


In order to derive explicit expressions for the tor- 
que and force on a rectangular strip, exerted by an 
incoming plane wave with the wave vector lying in 


2 
Fig. 1) it is necessary to know the wave w scattered 
by the strip. wy is obtained from that solution of the 
wave equation 


the xy-plane and the angle of incidence (F — v) (see 


Ap+hp=0 (p=Hot+y¥) (23) 
which satisfies the boundary conditions 
Oy/dy=0 for y=0 
p|oy or y (24) 


0): /Yr for roo. 


pr tl 
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The solution of these equations is well known (see 
[7]). Accordingly, we shall only —in order to fix 
the notation? — sketch the procedure. 

Eq. (23) may be separated in elliptic cylinder co- 
ordinates 


x=acosh “cost, y =asinh wsin®?, 


where a is half the width of the strip, the surface of 
which is given by «=O. The separated equations 
are second order differential equations of the 
Mathieu-type. The equation for ® defines, for all 
values of (ka)? =s, a complete set of eigenfunc- — 
tions 


even functions: 


Sem (#, s) = De De, (s,m) cos(n?), m=0,1,2,..., 
n=0 


(25) 


odd functions: 


Son (, 8) = oD Do, (s,m) sin(n?), m=1, 2,. 


n=1 
In the following we shall refer to these functions 
simply as Se(s,%) and So(s,%). The functions are 


6) 
59 So(s,#) =1 


for ?=0. For the “-equation, which may be derived 
from the equation for # by the substitution 0 =i, 
one gets the corresponding solutions Je,,(s, “) and 
Jo,,(s, 4). For large r these functions asymptotically 
approach 


normalized so that Se(s, 7) =1 and 


Eaengaik re (2m+1) 
4 


In order to satisfy the boundary conditions eq. (24) 
it is necessary to introduce a second set of indepen- 
dent solutions to the equation for w. Here we choose 
the functions denoted by Ne,,(s, 1) and No» (s, “) 
which asymptotically behave like 


ve im kr=— 7 (2m+l). 


The linear combinations 

Ho» (s, “) =Jom(s, “) +iNom(s,“), (26) 
analogous to the Hankel-Functions, will satisfy the 
boundary condition on wy for rc. 


In this notation the required solution to eqs. (23) 


and (24) is, [7]: 


YP =PotY; 
go< Van my im [rane Sent 9) Jem (s, 4) + 
hae Mi, 
Som(s,v) v) 
ys o Som(s, 9) Jom (s, 19), (27) 


2 Our notation follows that of “Tables relating to 
Mathieu-Functions” as closely as possible [8]. 
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1 3 hiee! 
(eo 


gers REDO (S52) 
y= —V8a "Sy i” te 
3 ba. Ms, 


X Som(s, 7) Hon (s, “), 


2% 4 
where Me p [Sem (s, 7) ]? dd 

0 

27 

and M®, = [ [Som(s, 9) ]? dd 

6 
and fo, m= Nom (s, 0) with “ meaning d/dw. It is seen 

Jom(s, 0) 


that also fo,» and M-° are functions of s= (ka)?. 


5. Explicit expressions for the torque and force 
on a strip 


With the expression for the scattered wave given 
by eq. (27) we easily get an explicit expression for 
the torque on a Rayleigh-strip. Introducing eq. (27) 
in eq. (8), using the asymptotic values of the radial 
functions Je, Jo and Ho together with the fact that 
the Mathieu-functions form an orthogonal set of 
eigenfunctions we get for the numerical value of H: 


Saris. e) “ Sree) 


po Vv 


ae M°, 


yy. 8 Qo Uo? 
eet ep) 
™m 


(28) 


where we have introduced wy” instead of u,,. 

__ Similarly we could get an expression for the force 
K by means of eq. (27) and eq. (22). However, 
using the known scattering cross-section for the strip 
(see [7]) together with the general expression eq. 
(9) we immediately have 


1 [Som(s, 0) ]*- (29) 


mi+fom My, 
In the next section we shall give the results of nume- 
rical calculations based on the expressions eqs. (28) 
and (29). 

We have further found an expansion of the tor- 
que H in powers of s= (ka)?. However, as the ex- 
pansion converges very slowly, we shall only give 
the first terms for the special case v = 45°. We find 


=e = 1 Teg k 
H=xoyugety|t+(5 pink) (ka)?+... 
(30) 


with y = 1.7811. As will be seen by comparison with 
the numerical calculations in the next section, this 
expression is only valid when ka © 0.5. 


6. Numerical results 


By means of eqs. (28) and (29) we have per- 
formed the following calculations: 
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a) The torque H as a function of ka for the con- 
- _ »} =45°, 


b) The torque H as a function of the angle of inci- 


stant angle of incidence ( 


dence (F — ») for constant value of (ka)? =1.4. 


c) The force K as a function of ka for v= = : 


The functions So(s, v) and = So(s,») were cal- 
v 


culated by means of the coefficients Do, (s,m) (see 
eq. (25)) tabulated in [8]. However, for some 
values of s it was necessary to make a graphical 
interpolation of these coefficients. This was done by 
means of the coefficients F”’(=Do,(s,m)) tabulat- 
ed in [9]. The quantities /,,,, were taken from [8]. 


Table I. The torque H on a rectangular strip (v=45°). 


(ka)? | m=1 | m=2 |) m=3 | H 

0 1.0000 
0.25 | 1.120 
0.36) 1.1875 fe. PST 
0.60} 1.2373 1.2373 
0.80] 1.2278 | |. 1.2278 
1.00} 1.1739: | 0.0053 | —0.0010 | 1.1782 
1.2 | 1.0845 | 0.0076 | — 0.0014 1. 0907 
1.4 | 0.9730 0.0104 | — 0.0019 0.9816 
1.6 | 0.8541 | 0.0137 | — 0.0025 0.8653 
2.0 | 0.6315 | 0.0215 | — 0.0039 0.6491 
2.5 | 0.4203 | 0.0341 | — 0.0061 0.4483 
3.0 | 0.2814 0.0496 | — 0.0088 0.3222 
4.0 | 0.1346 | 0.0885 | — 0.0154 0.208 
5.0 | 0.0705 -| 0.1337 | — 0.0236 0.181 

6.0 | 0.0397 | 0.1744 | — 0.0332 0.181 

7.0 | 0.0236 | 0.1982 | — 0.0439 0.178 
8.0 | 0.0146 0.1999 | — 0.0555 0.158 
9.0 | 0.0094 | 0.1815 | —0.0676 0.123 
14.0 | (including m = 5) — 0.0610 


The torque H as calculated by eq. (28). The terms in eq. (28) with 
m =1,2,3 are shown, together with the final value. H is measured 
in the unit given by eq. (31). 


The results of the calculations for H are shown in 
Table I. To give an impression of the convergence 
of eq. (28) we have shown the contributions from 
the terms m=1, 2 and 3. It will be seen that the 
convergence is rather good. Further, for (ka)?=14. 
we have calculated the terms up to and including 
m=5. We find H = —0.0610 in the unit given by 
eq. (31). The surprising fact, that H may become 
negative, was first demonstrated by Awarani [3] 
for the circular disk. For this value of ka, Levine’s 
asymptotic expression gives a vanishing torque; 
using Keier’s methods, however, correction terms 
may be deduced securing a better agreement. For 
still larger values of ka the torque H oscillates bet- 
ween positive and negative values. In Fig. 2 we have 
plotted H vs. ka. The unit used for H is the well 
known low frequency limit: 


(31) 


Hy =X Gy Up” a” sin v cos v 


2S) ——— 
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with v=45°. For ka<0.5 the curve is based on 
the expansion eq. (30). 


14 J 


O04 


0.2 


0 
-01 


Fig. 2. H vs. ka for a) the circular disk and b) the rec- 
tangular strip, both in units of the low fre- 
quency limit. The angle of incidence is 45°. 


In Fig. 3 we have plotted H vs. v for (ka)? =1.4. 
H is here measured in units of x Q) up? a2. In Fig. 4 
we show H vs. v with H measured in the unit given 
by eq. (31). In this case H should be equal to 1 for 
all values of v if eq. (31) were valid for all fre- 
quencies. t 

Finally in Fig. 5 we have plotted the force K vs. 
ka for v=2x/2. K is here measured in units of the 
high frequency limit 4 09 U2 a. Here, too, we have 
neglected terms with m = 4 in eq. (29). However, 
this is not admissible for large values of ka. Accor- 
dingly the curve for ka®&3 is shown broken. 


: 


| se 
0.1 | a 
0 08 pp 30° 60° 90° 
0° 30° 60° 90° 
vy 
ae Fig. 4. 
Fig. 3. H vs. v for the rectangular 


strip, (ka)? =1.4. A is 
here measured with the 
unit J Q9 Up? a” sin v Cos v. 


H vs. v for the rectangular 
strip, (k a)? =1.4, The unit 
for H is 1 Qp Up? a”. 
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Fig. 5. The force K vs. ka for a) the circular disk and 
b) the rectangular strip, both in units of the 
~ high frequency limit. 


7. The circular disk 


In Fig. 2 we have also drawn a curve showing the 
torque on a circular disk? as a function of ka for 
the case of v=45°. It is based on the series expan- 
sion (valid to ka~1) of Koranr and the numerical 
calculations of Kawat, supplemented by our own cal- 
culations, made by means of “Tables of spheroidal 


functions” [10], [11]. The unit for H is the value 


given by Kénies formula, ine. 5 Qo te? a>. It will be 


seen that the circular disk shows qualitatively the 
same behaviour as the strip does, the main differen- 
ces being a somewhat larger value of the maximum 
torque and a displacement of the curve towards lar- 
ger values of (ka). 

The problem of the scattering cross section of the 
circular disk has already been treated numerically 
by Bouwxame [12] and by Levine and Scuwincer 
[13]. From their results K is determined by eq. (9). 

In Fig. 5 we have, for comparison with the strip, 
shown the curve thus obtained from Bouwxamp’s 
data. It will again be noted that the two curves are 
qualitatively similar. 


8. Discussion 


8.1. The torque H. 


It will be seen from Fig. 2 that the torque H for 
ka # 0.1 shows considerable deviation from the 
simple low frequency formula eq. (31). After pas- 
sing through a maximum around ka~0.8 it decrea- 
ses rapidly towards zero. For ka>1.2 the exact 
value of H is less than the value given by eq. (31). 
It is evident that we have to do with a resonance 
phenomenon of the type well known from related 
diffraction problems (see [12], [13], [14], [15]). 
The ka value for the first maximum corresponds to 
a ratio between the with of the strip and the wave- 
length of 2a/A~1/4 whereas the next maximum — 
or point of inflection—corresponds to 2a/A~3/4. 


3 For the circular disk a is the radius, 
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The deviation of H from the value of eq. (31) for 
ka=—0.8 amounts to ~25%, a rather large and pre- 
sumably measurable effect. In spite of this the effect 
does not seem to have been found experimentally. 


The torque on a Rayleigh-strip has been. measured 
for varying values of 2a/d by Harrmann and Mor- 
TENSEN [14], [15]. For a constant sound field they 
have measured the ratio Ep/E,, where Ep is the 
sound intensity as measured by a Rayleigh-strip, and 
Ex is the sound intensity as measured by a suitable 
sound-radiometer. Their curve, which, however, is 
based on only four points, does not show any trace 
of a resonance. The decrease in numerical value of 
H occurs already for ka 0.6, whereas below this 


value H appears to be constant, In comparing this 
with our calculations it is however important to 
have in mind the following points. 


a) We have neglected a possible vortex and tur- 
bulence formation, b) the sound field used by Harr- 
MANN and Mortensen is not entirely monochromatic 
and c) their measurements consist of only four 
points rather widely spaced with regard to the value 


of ka. 


Concerning the first point, vortex and turbulence 
formation, it is very difficult to form a definite opi- 
nion as to the possible consequences of these effects. 
Experimentally Merrineron and Oarttry [16] and 
Scorr [17] at low frequencies (ka®0.2) find for 
the circular Rayleigh-disk a definite, essentially con- 


stant deviation of H from K6n1c’s value, amounting 
to ~10%. They attribute this deviation to the for- 
mation of vortices. This formation is demonstrated 
by photographs. The experiments are made at very 
low frequencies. The effect of vorticity presumably 
will diminish at higher frequencies, because the 
thickness of the boundary layer will be of the order 
of magnitude //w , where » is the kinematic visco- 
sity and w the angular frequency. At frequencies, 
where it is necessary to correct for the diffraction 
effects, the boundary layer thickness will be very 
small indeed, and one is inclined to assume that vor- 
ticity and other effects of viscosity will be of minor 
importance. With regard to b), Harrmann and Mor- 
TENSEN state that their sound field, produced by an 
acoustic air jet generator, contains higher harmo- 
nics. This will have the effect that the resonance 
maximum is lowered and moved to lower values of 
ka and thus diminish the discrepancy between our 
curve and that of Hartmann and Mortensen. How- 
ever, the experimental material does not seem to be 
conclusive. Finally the consequences of point c) are 
obvious. 


For the torque on the circular Rayleigh-disk, Bar- 
nes and Wesr [18] have found experimental results, 


. 
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which qualitatively agree with the curve shown in 
Fig. 2. They find for the deviations from Kénic’s 
value, the latter realized by means of a standard disk 
with ka very small, the following results: 
0.1<ka<0.6: uniform distribution of deviations 
around a mean value zero, within 
+ 12%, 
: all deviations positive, maximum 
deviation for ka~1.2 amounting 
to 20 — 30%, 
all deviations negative and of order 
10-20%. 
This is qualitatively, and also with regard to 
order of magnitude, in agreement with the theoreti- 


0.9<ka<2 


kax2A: 


cal values for H. 

Also with regard to the angular dependence of H 
we find rather large deviations from eq. (31). This 
is most clearly seen from Fig. 4. Here H should be 
constant and equal to 1 if eq. (31) were valid for all 
frequencies. Instead of this, H varies between 0.87 
and 1.1. It may be possible to investigate the theory 
through this angular dependence of H. 


8.2. The force K at perpendicular incidence. 


The force K on a strip, exerted by the perpendicu- 
lar incoming wave, shows a series of maxima and 
minima, the first maximum occurring at ka~1.7 
corresponding to a ratio of 2a/A~1/2 between the 
width of the strip and the wavelength. This is ap- 
proximately twice the corresponding value for the 
torque. For the circular disk Harrmann and Mor- 
TENSEN have investigated the dependence of the force 
on ka [14], [15]. They measured the ratio Fp/E, . 
Here Fy is the intensity of the sound field measured 
by the radiometer, whereas Ex is the intensity of the 
same sound field mesured by the torque on a Ray- 
leigh-strip, ka being small enough for eq. (31) to 
be valid. Their results are in excellent agreement 
with the curve in Fig. 5, although some of their cur- 
ves for Ep/Ep vs. ka for large ka tends to a value 
<1. However, this has hardly any direct relevance 
in this connection. Harrmann and Morrensen find 
maxima and minima of Ep/Ep — i.e. of K — for 
the following ka values [14] 

Maxima: ka~5.7; 8.8; 

Minima: ka~4.4; 7.5; 10.7. 
Evidently, a first maximum has escaped observation. 
This is due to the fact that the measurements are 
restricted to kaw 3. 

This is now to be compared with the theoretical 
curves based on the results of Bouwkamp [12] 
Fig. (5). Here the maxima and minima are given by 

Maxima: ka~2.4; 5.5; 9.1; 
Minima: ka~4.3; 7.3. 
The agreement is excellent. In view of this, it ap- 
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pears reasonable to expect a similar agreement in 
the case of the force on a strip. With regard to the 
force K at perpendicular incidence it thus appears 
safe to neglect vortices and turbulence. Whether this 
holds in general may be dubious (compare [6]). 
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ABAQUES POUR LE CALCUL DE LABSORE TION 


ACOUSTIQUE DES PAROIS PERFOREES 


par P. Lizenarp 


Chef de Division Acoustique, Office National d’Etudes et de Recherches Aéronautiques, 
ChAtillon-sous-Bagneux, France 


Sommaire 


Les calculs des structures perforées tant continues que cloisonnées (employées comme 
absorbants aux basses fréquences) sont assez complexes et fastidieux pour des applications 
nombreuses et diverses. C’est pourquoi l’auteur a établi des abaques permettant un calcul 
rapide et d’un mode d’emploi plus général. Des vérifications expérimentales (en cours) 
donnent déja de bons résultats. 


Zusammenfassung 


Die Berechnung von gelochten Strukturen mit sowohl durchgehendem als auch unterteil- 
tem Hohlraum, wie sie fiir Tiefenabsorber angewandt werden, ist fiir die zahlreichen An- 
wendungen ziemlich kompliziert und langwierig. Deshalb hat der Verfasser Rechentafeln 
aufgestellt, die eine schnelle und allgemein anwendbare Berechnung gestatten. Gegenwartig 
durchgefiihrte experimentelle Nachpriifungen liefern gute Ergebnisse. 


Summary 

Calculations of the low-frequency absorptions of the widely-used perforated structures 
(whether continuous or in the form of bulkheads) are complex and difficult. For this reason 
the author has set up abaci suitable for rapid calculation and general use. Tests now in 
hand of these are already showing good results. 


1. Introduction 


On sait que les absorbants usuels, fibreux ou 
poreux, sont assez peu efficaces en basse fréquence, 
et qu'on emploie de plus en plus dans ce domaine 


des panneaux perforés a faible taux de perforation, 
formant autant de résonateurs d’Hetmuotrtz, dont les 
cavités sont les volumes d’air compris entre la paroi 


perforée et le mur rigide. 


4 


. 


| 
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Ces structures ont l’avantage d’étre assez peu coi- 
teuses, peu fragiles, de supporter les gaz chauds et 
humides et de pouvoir étre accordées sur des fré- 
quences méme infra-sonores, en donnant aux volu- 
mes d’air des dimensions suffisantes. 

Ces qualités nous ont fait utiliser ce type d’absor- 
bant pour les sorties de gaz de laboratoires de com- 
bustion construits par |’O.N.E.R.A., qui sont sources 
de bruits trés intenses et doivent recevoir une in- 
sonorisation trés soignée en raison de leur proximité 
de zones résidentielles. 

On sait qu’il existe deux types de structures perfo- 
rées: les structures continues, ot les divers volumes 
placés derriére chaque orifice ne sont pas séparés, et 
les structures cloisonnées, analogues 4 des résona- 
teurs juxtaposés. 

Les caractéristiques acoustiques de ces structures 
sont calculables, mais les calculs fastidieux que cela 
entraine nous ont conduits a\rechercher des métho- 
des graphiques. 


2. Calcul des structures perforées 


L’étude de ces parois perforées a été faite par 
plusieurs auteurs, expérimentalement et théorique- 
ment. Pour l’étude d’insonorisation signalée plus 
haut, nous avons utilisé les formules établies par 
J. Brittourn [1], que nous rappelleron ci-dessous, 
ainsi que les notations: 

A coefficient d’absorption cherché (rapport d’éner- 
gies), 

c vitesse du son, 

fy fréquence propre d’un résonateur a la tempéra- 

ture ordinaire, 

fréquence propre a la température ¢° C, 

température centésimale, 

rapport de la fréquence d’excitation a la fré- 

quence propre, 

résistivité relative de la paroi perforée, ou rap- 

port de la résistivité des trous a la résistance 

spécifique de l’air, 

coefficient d’épaisseur de la structure, 

taux de perforation de la paroi, 

épaisseur de la lame d’air entre paroi et mur, 

longueur de chaque col de résonateur, c’est-a-dire 

généralement, épaisseur de la paroi perforée, 

Uv longueur corrigée du col, 

0 angle d’incidence du son. 

La longueur corrigée I’ se calcule a partir de / 
par des formules assez empiriques: 

pega > pour un’ col unique circulaire 
2 de rayon r 


ou [=140,96 Vs 


Spa i ee fre pour plusieurs cols circulaires 
t voisins. 


bes 


Se a oe 


pour un col de section s 


= 
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La fréquence propre est donnée par 
= 


ease ai Ga? 
lo= 20 Ver 
avec c=34000cm/s a 20°C et le coefficient d’ 


épaisseur par 
K- Ve (2) 


Ces deux formules seront utilisables pour les 
structures continues comme pour les structures cloi- 


(1) 


sonnées; dans ce dernier cas on tiendra simplement 
compte de l’épaisseur des parois transversales, dont 
la tranche doit étre déduite de la surface totale de la 
paroi perforée quand on calcule o. 

La vitesse du son c variant avec la température, 
il y aura lieu d’en tenir compte dans le calcul de f 
si la température est notablement différente de la 
température ordinaire, comme cela se produit a la 
sortie de chambres de combustion. 

La fréquence de résonance sera calculée par 


t 
i= do yj + 973° 


Les calculs faits par J. Brrttourn [1] aboutissent 
aux expressions suivantes: 
Pour les structures continues 


4H cos 9 


(3) 


ce (7 cos9+1)?+N,? (4) 
avec N,= 7-0 cos 8— cot K 2 cos 0. 
Pour les structures cloisonnées 
4 H cos 0 (5) 


igo Sa es oe 
(H cos9 +1)? +N,? cos? 6 
avec 


Np= 52 — cot K.Q. 


La validité de ces formules a déja été vérifiée par 
leur auteur. 

Un travail analogue de Incarp et Botr [2] abou- 
tit 4 des expressions presque identiques, dans les- 
quelles le terme en 1/K des expressions NV, et Vy, est 
remplacé par cot K. Les deux formules sont prati- 
quement semblables pour K petit; méme si K prend 
des valeurs notables, l’écart obtenu pour A n’est pas 
trés grand (par exemple pour K=0,8, H=1, 
cos9=0,5, on obtient respectivement par les deux 
formules A= 0,38 et 0,36 pour 2=1; 0,84 et 0,9 
pour 2=2; 0,24 et 0,35 pour Q=4). 

Ces formules, sans étre d’une grande complexité, 
entrainent cependant des calculs numériques fasti- 
dieux, surtout quand les études d’un projet obligent 
a les appliquer trés souvent en faisant varier l'un 
ou l'autre paramétres. 
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Il faut remarquer de plus que le rapport H ne 
peut étre connu qu’expérimentalement et qu’alors 
les expressions (4) et (5) doivent étre résolues en H. 

Les graphiques donnés dans la référence [2] ne 
correspondent qu’a quatre valeurs du rapport H 
(désigné alors par 9 ) et ne permettent guére l’inter- 
polation. 

Nous avons donc été amenés a établir des abaques 
plus généraux, permettant le calcul rapide, avec une 
précision analogue a celle des mesures acoustiques, 
des expressions (1) a (5). 


3. Caleul graphique de / et K 


Les équations (1) a (3), sous forme de produits, 
se prétent bien a la réalisation d’un abaque a points 
alignés, les échelles étant paralléles et a graduations 
logarithmiques. 

La construction de cet abaque ne présente pas de 
difficulté et son emploi est immédiat: les valeurs de 
o et l’ ayant été déja calculées pour un systéme de 
résonateurs donné ou supposé, on joint par une droite 
les points correspondants de ces deux échelles, et le 
point d’intersection de cette droite avec l’axe o/I’ est 
joint 4 la valeur donnée ou choisie de e (l’axe o/l’ 
n’a pas été gradué). 


o, Kk t a/v" fo 
001-0) 
Hz 
10 
002 20 
30 
003 
température t 
0,04 
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On obtient ainsi les valeurs de K et de f a 20°C 
sur les axes correspondants, Un axe gradué en tem- 
pérature permet d’obtenir /; a partir de /g90. 


4. Calcul graphique de l’absorption des 
structures continues 


Le calcul des expressions (4) et (5) est plus 
compliqué et a été résolu par des réseaux de cour- 
bes. 

Les produits 2cos@ et Hcos@ sont calculables 
par points alignés: il a semblé plus commode de 
marquer cette correspondance par des réseaux de 
courbes, pour obtenir un graphique homogeéne. Avec 
des graduations logarithmiques de module con- 
venable, ces courbes sont des droites 4 45°. 

La partie essentielle de l’abaque porte les courbes: 


2 
pe eet Ke OICOS 0) : 
K 


K prenant les valeurs successives d’un parametre 
variable, et les courbes: 
_ 4H cos 8 

A 


N? (H cos0+1)?; 


A étant le paramétre. 


f ie Les abscisses sont graduées en 


% valeurs de {2cos@ ou H cos? 

pera ac et N?, 
be Les premiéres courbes présen- 
a as tent une variation trés rapide au 
F 40 voisinage de V=0. L’échelle de 


N? ne peut étre logarithmique, 
en raison de la forme trop 
pointue que prennent alors les 
courbes K=const, ni linéaire, 
en raison de |’étalement des va- 
leurs que prend ce nombre. 
oes On a pris une échelle arbitraire 
pour NV, telle que le courbe 
K =0,5 soit représentée par deux 
v droites. A partir de ces droites, 
on a construit une courbe d’éta- 
lonnage des longueurs de l’axe 


3 des ordonnées en fonction des 
valeurs de N? (non reproduites) et 
2 la, les autres courbes K = const. 


L’usage de l’abaque, pour les 

structures continues, s’en. déduit 
immédiatement. 

La partie (1) du graphique 

permet de passer des valeurs de 

H et 2 a celles de H cos 0 et 

Q cos 0, reportées sur le méme 

axe des abscisses de la partie (2). 


Tableau I. Abaque pour le ealcul de la fréquence propre et du coefficient 


d’épaisseur des parois perforées. 


oe 
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Partant de Hcos@ et 2cos@, K étant connu, 
Vhorizontale de l’intersection de K avec cos @ coupe 
H cos@ sur la courbe de la valeur cherchée de A. 

Inversement, on peut a partir de A donné, déter- 
miner K ou H. ; 

On remarquera que la valeur maximum de l’ab- 
sorption A est obtenue quand on se place le plus 
bas possible sur le diagramme, soit pour V=0: le 
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Nous préciserons plus loin comment on peut cal- 
culer l’intervalle quadrantal. 


5. Calcul graphique de l’absorption des 
structures cloisonnées 


Les mémes abaques sont utilisables pour le cal- 
cul des structures cloisonnées moyennant une adjonc- 


Nn2 
500, 
pr | | @ 
4 
| 
300 | ~ 
e 
200) 7 
2 K- 
150 ° + 
SE <S . 
100 B Ry 
OxXe wW 
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LI 50 Z 
A io > 
40 t SS i 06 
Y 30 ow 0,7 
[| 08 | 
a ot 
VW) tA 15 ow 
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ALA 8 ro mn : S 
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EE SRE g 
iS DEG 08 0 
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ee! 03 96) Qo 
= s OOS) 
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Tableau II. Abaque pour le calcul de l’absorption des parois perforées. 


maximum de A est donné en fonction de H cos @ 
par les deux graduations en regard de l’axe des 
abscisses. 
La forme des courbes K = const permet aussi de se 
rendre compte de l’acuité des courbes d’absorption 
en fonction de la fréquence quand K est petit. 


tion: en effet les équations (5) permettent de cal- 
culer No? a partir de 2 et K comme cn a calculé 
N,2 a partir de 2cos@ et K, et les graduations 2 
et Qcos@ sont les mémes. On calcule également A, 
a partir de H cos@ et N.* cos? @ comme 4A, a par- 


tir de H cos @ et N,?. 
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La partie (3) de l’abaque permet le calcul de 
N? cos? @ par des courbes continues (dont la forme 
dépend de l’échelle arbitraire de NV”). 

L’emploi de l’abaque pour les structures cloison- 
nées s’en déduit ainsi: pour chaque valeur de Q, la 
courbe K étant choisie sur la partie (2), on déter- 
mine l’horizontale d’ordonnée NV?; en suivant la 
courbe correspondant a cette valeur jusqu’a |’abs- 
cisse de @ dans la partie (3), on détermine lordon- 
née V? cos? 0, la valeur de A correspond a celle des 
courbes (interpolée au besoin) passant par l’inter- 
section de cette ordonnée avec l’abscisse H cos 0. 


6. Etude de l’intervalle quadrantal 


On sait que, sur une courbe de résonance, |’inter- 
valle quadrantal est l’écart des deux fréquences pour 
lesquelles ’amplitude (ou l’absorption dans le cas 
présent) est égale a la moitié du maximum. 

Dans le cas des structures continues, si H , 0 et K 
sont fixes, la valeur de Hcos@ détermine sur le 
méme axe la valeur du maximum de 4, comme on 
Ya indiqué plus haut, et l’intersection de la courbe 
ay Ae avec la droite d’abscisse H cos @ définit l’or- 
donnée des deux points qui, sur les deux branches 
de la courbe K , définissent les fréquences cherchées, 
dont l’écart en savarts est directement indiqué sur 
laxe des abscisses. 

Pour les structures cloisonnées, le maximum de 4 
est indépendant de 0; V intersection de la droite 
d’abscisse H avec Ja courbe = Amax donne la valeur 
de N? cos? 9, d’ou on remonte par la partie (3) de 
Pabaque, a NV, qui définit ’horizontale coupant la 
courbe K aux deux points, dont l’écart en fréquen- 
ces est la valeur cherchée. 

On voit facilement par ces constructions que 
Pintervalle quadrantal diminue quand 0@ croit, pour 
les structures continues et augmente au contraire 
pour les structures cloisonnées. 


7. Emploi des abaques a partir de données 
expérimentales 


Il est parfois nécessaire de déterminer expérimen- 
talement les caractéristiques d’une structure et d’en 
déduire en particulier le coefficient H de résistivité 
relative, ou plutdt le produit H cos 0. 

Si le maximum d’absorption (sommet de la 
courbe de résonance) est obtenu avec une précision 
suffisante, les graduations en regard de l’axe des 
abscisses du diagramme (3) donnent immédiate- 
ment cette correspondance; deux valeurs de H 
toutefois restant possibles: le choix entre elles sera 
fait en se rappelant qu’un orifice assez grand pré- 
sente une résistivité trés faible, et qu’il faut un 
amortissement assez important par grillages ou pro- 
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duits fibreux dans les cols de résonateurs ou des 
cols extrémement fins, pour obtenir une résistivité 
supérieure a 1. La construction de la courbe d’ab- 


‘sorption permet ensuite une vérification. 


Les mesures expérimentales seront faites géné- 
ralement sous incidence normale, ou au moins de 
valeur connue: on obtiendra ainsi le coefficient H. 
Si le son est diffus, on aura une premiére approxi- 
mation en admettant un angle d’incidence équiva- 
lent de 55°. 

Il est parfois plus précis, ayant tracé une courbe 
de résonance, de déterminer H a partir de J’inter- 
valle quadrantal mésuré, car une erreur sur le 
maximum An ax de la courbe se traduit par une 
erreur beaucoup plus faible sur la longueur de la 
courbe a |’amplitude ai I] suffit alors de pla- 


cer entre les deux branches de la courbe K (K. 
étant connu) une longueur égale a cet écart en fré- 
quences, au moyen de la graduation en savarts, tra- 
cée au bas du diagramme, et reportée sur rhodoid 
par exemple. 

Si K n’est pas connu, H étant calculé a partir du 
maximum de la courbe d’absorption, on cherchera 
quelle est celle des courbes K qui admet entre ses 
branches l’intervalle quadrantal mesuré, correspon- 
dant aux valeurs H (ou Hcos@) et = Amax déja 


connues. 


8. Verifications expérimentales 


Une étude de résonateurs est en cours, pour le 
projet d’insonorisation de laboratoires de com- 
bustion. Cette étude se base sur l’emploi d’une part 
d’éléments de petite dimension (20cm 20cm) 
étudiés sous incidence normale au tube a ondes sta- | 
tionnaires, d’autre part de parois de briques per- 
forées disposées dans une galerie, puis dans une 
salie réverbérante, sous incidence aléatoire. 

Les premiers résultats montrent quelques légéres 
divergences entre les résultats théoriques des formu- 
les ci-dessus et les résultats expérimentaux, mais 
sans que l’écart constaté puisse étre génant pour 
Pétablissement d’un projet. L’emploi des abaques a 
permis tous les calculs trés rapidement, et les résul- 
tats de ces mesures d’absorption feront l’objet 
d’une communication ultérieure. 


(Recu le 14 Fevrier 1957.) 
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DAS AUSGLEICHSVERHALTEN VON GEIGEN 


UND SEINE BEZIEHUNG ZU DER RESONANZKURVE 


Von W. Lorrermoser 
Mitteilung aus der Physikalisch-Technischen Bundesanstalt Braunschweig 


Zusammenfassung 


Zwei neue mittelmaBige Geigen und zwei hervorragende Stradivari-Geigen wurden am 
Steg impulsartig bei verschiedenen Frequenzen erregt. Die Einschwingvorginge verlaufen 
je nach Lage der Tragerfrequenz zu den Resonanzstellen verschiedenartig. Die Auswertung 
nach E. Meyer-Hetmsotp ergibt, daB die Stradivari-Geigen schneller einschwingen als die 
neuen. Dies wird dadurch erklart, da8 trotz geringerer Dampfung und daher groferen 
Wirkungsgrades der altitalienischen Geigen diese eine gré8ere Zahl von Resonanzen haben, 
die bandfilterartig dicht nebeneinander liegen. 


Summary 


Two new violins of medium quality and two excellent violins of Stradivarius were excited 
on the bridge by impulses at several frequencies. The transient phenomena differ accord- 
ing to the relation of the carrier-frequency to the resonances. The interpretation according 
to E. Meyver-Hetmsorp shows that the Stradivarius violins have shorter transients than the 
modern ones. This can be explained by the fact that the old Italian master-pieces have, in 
spite of a less damping and therefore greater efficiency, resonances lying close together 
like band-filters. 


Sommaire 


On a excité par impulsions sur le chevalet, 4 des fréquences différentes, deux violons 
modernes de qualité médiocre et deux violons supérieurs de Stradivarius. Les phénoménes 
transitoires de mise en vibration different suivant la position de la fréquence porteuse par 
rapport aux résonances. L’étude faite suivant le procédé de E. Mryer-Hetmsorp montre 
que les Stradivarius entrent plus rapidement en vibration que les violons modernes. Cela 
s’explique du fait que les anciens violons italiens, bien qu’ayant un amortissement plus 
faible et par conséquent un rendement plus grand, possédent un nombre plus grand de 


résonances qui sont situées trés pres l’une de |’autre 4 la maniére d’un filtre de bande. 


1. Problemstellung 


Aus Gesprachen mit Geigenbauern und Geigen- 
spielern tiber die Gite von Geigen wird deutlich, 
eine wie groe Bedeutung der Ansprache einer Geige 
beizumessen ist. Eine Violine, welche leicht anspricht, 
bei welcher sich also die Teilténe beim Anstrich 
schnell zur vollen Amplitude entwickeln, wird in 
jedem Fall besser beurteilt als eine solche, die 
»schwer kommt“. Es ist natiirlich vorteilhaft, wenn 
die Geige in méglichst vielen Tonlagen leicht an- 
spricht, dagegen ungiinstig, wenn es eine grofere 
Zahl von Toénen gibt, bei denen die Ansprache mih- 
sam ist. Ganz frei von derartigen Fehlerstellen ist 
keine Geige, doch kennt der geiibte Spieler die 
Schwichen seines Instruments und sucht sie durch 
besonderen Bogenstrich auszugleichen. 


H. v. Hermnorrz [1] teilt in seiner Lehre von 
den Tonempfindungen Beobachtungen mit, welche 
in diesem Zusammenhang von groftem Interesse 
sind. Er studierte die Schwingungsform angestriche- 
ner Saiten und fand, daB es schwierig war, stabile 
_ Schwingungen iiber langere Zeit auf einem Mono- 


chord zu erzeugen. Die Grundform der Sagezahn- 
schwingung wurde immer wieder unterbrochen und 
plotzlich auftretende Verschiebungen und Verande- 
rungen der Schwingungsform traten auf. Auf einer 
Guadagnini-Geige dagegen konnten mit Leichtigkeit 
ruhig verlaufende, periodische Schwingungen her- 
vorgebracht werden. 

In seiner grundlegenden Arbeit tiber die Aus- 
gleichsvorginge der Akustik hat H. Bacxnaus [2] 
u.a. die Entwicklung der Teilténe beim Anstrich 
einer Landolfi-Geige auf der Grundfrequenz a’ 
= 435 Hz verfolgt. Wahrend sich die Teiltone der 
1. und 2. Ordnung stetig bis zur vollen Amplitude 
in etwa 100 ms entwickelten, vergroSerten sich die 
der héheren Ordnungen in mehreren Schwebungen, 
erreichten jedoch nach etwa 10 bis 20 ms gréfere 
Amplituden als die erstgenannten Teilténe. Betracht- 
liche Schwankungen der Obertonintensitéaten wurden 
auch wahrend des sogenannten stationaren Teils des 
Bogenanstrichs festgestellt. 

Die Resonanzkurve einer Geige zeigt, wie be- 
kannt [3], eine groBe Zahl von Maxima mehr oder 
weniger groBber Amplitude, Frequenzbreite und -lage. 
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Je nachdem, ob die Harmonischen der Saitenschwin- 
gung auf Maxima oder Minima fallen, werden sie 
starker oder schwacher vom Resonanzkorper abge- 
strahlt. 

Die Dampfungen von Geigenhélzern hat zuerst 
E. Routorr, spater E. Skuprzyx [4] gemessen, wo- 
bei letzterer auf die damit im Zusammenhang stehen- 
den Ausgleichsvorgange aufmerksam gemacht hat. 
F. A. Saunpers [5] hat das Ausklingen beobachtet 
und versucht, in Beziehung zur klanglichen Giite zu 
setzen. J.G.Hetmporp [6] hat die Einschwing- 
vorgange von Lautsprechern gemessen und brachte 
deren Dauer in eine enge Beziehung zu ihrer Uber- 
tragungskurve. Oszillographisch wurden dabei Ein- 
schwingvorgange beim Einschalten von Sinusspan- 
nungen verschiedener Frequenz registriert. Zur Aus- 
wertung wurden die Hiillkurven zugrunde gelegt und 
der Einschwingvorgang als praktisch beendet an- 
gesehen, sobald die Amplitude in den Bereich von 
+ 20% des stationaren Wertes gelangt war. Die Ein- 
schwingzeit wurde als positiv bewertet, wenn sich die 
Hillkurve asymptotisch diesem Wert naherte. Wenn 
aber die Amplitude vortibergehend grofere Werte 
erreichte und sich darauf erst dem stationadren 
naherte, wurde die Einschwingzeit als negativ an- 
gesehen. Das Ergebnis zeigte, wie zu erwarten war, 
dal} die Einschwingzeit bei den Resonanzspitzen am 
langsten, und an den Stellen, wo kein ausgepragtes 
Maximum vorhanden ist, am kiirzesten ist. 

In seiner Einfthrung in die Akustik gibt F. Tren- 
DELENBURG [7] einige typische Falle wieder, aus 
denen ersichtlich ist, in welcher Weise sich der Ein- 
schwingungsvorgang und seine Hiillkurve verandert, 
wenn die rechteckig geschaltete Anregungsfrequenz f 
in verschiedenem Verhaltnis zu der Resonanz- 
frequenz fy des angeregten Systems steht. Ist f=fy. 
so wachst die Amplitude mit (1 —e °*) an. Ist f+ fo, 
so entwickelt sich die Amplitude nicht kontinuierlich, 
sondern erreicht voriibergehend grofere Werte als 
im eingeschwungenen Zustand. 

Bezeichnet k(t) die auBere, auf das System wir- 
kende sinusformige Kraft K (t) = Ky sin @ ¢, so lautet 
bekanntlich die Schwingungsgleichung 


+Dx=K(t) ? 


wobei M die Masse, R die Reibungskraft, D die 
Riickstellkraft des angeregten Systems ist. Die L6- 
sung der Gleichung fihrt zu 


A=Apsin(wi+) +A’e~*sin(w@yt+y) , 


d.h. es treten zwei Summanden auf, von denen der 
erste die Schwingungen im stationaren Zustand, der 
zweite die frei abklingenden Schwingungen mit der 
Dampfungskonstante 6 und der Eigenfrequenz , 
des Systems darstellt. Man erhalt den Einschwing- 
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vorgang durch Superposition beider Schwingungen 
nach Amplitude, Frequenz und Phase. F. Vitsre [8] 
gibt dazu ein Verfahren an, mit dessen Hilfe der- 


' artige Vorgange graphisch ermittelt werden kénnen. 


Allgemein lassen sich Probleme dieser Art bekannt- 
lich durch die Laplace-Transformation behandeln, 
wortiber ein ausftihrliches Schrifttum [9] vorliegt. 

In dieser Arbeit werden einige Aufnahmen ge- 
bracht, an Hand derer das Ausgleichsverhalten von 
Geigen diskutiert wird. 


2. Die Resonanzeigenschaften der Geige 


Der nach Form und Aufbau komplizierte Kérper 
der Geige lat genauere Berechnungen der Frequenz- 
lagen und Dampfungen der Resonanzen nicht zu. 
Auf experimentellem Weg hat man ermittelt, daS 
solche teils vom Hohlraum, teils von den Platten- 
eigenschwingungen herriihren, wobei komplizierend 
hinzukommt, dafi die Decke zur Verstaérkung mit 
dem Bafibalken versehen ist und dafi sie mit dem 
Boden durch den unsymmetrisch eingespannten 
Stimmstock gekoppelt ist. AuSerdem wird durch die 
Saiten tiber den Steg ein statischer Druck auf die 
Decke ausgeiibt, wodurch sich deren Resonanzen ver- 
lagern. Decke und Boden sind am Rand durch die 
Zargen und Eckklotzchen eingespannt, wobei aber 
die Zargen selbst keine Knoten der Schwingungs- 
zustande bilden. Durch frithere Versuche [10] weib 
man beispielsweise, das die sogenannte Hohlraum- 
resonanz durch das Volumen des Hohlraums und 
seiner Offnungen in Gestalt der beiden f-Lécher her- 
vorgerufen wird, und daS mit dem Hohlraum die 
Platten, besonders der Boden eng gekoppelt ist. Labt 
man namlich den Stimmstock weg, so fehlt das 
Maximum der Hohlraumresonanz nahezu, ein Zei- 
chen daftir, da’ dieses Maximum nur dann entstehen 
kann, wenn der Boden frei mitschwingt. Deswegen 
sollte er im fraglichen Frequenzgebiet eine oder 
besser mehrere Resonanzen besitzen, damit das ge- 
wiinschte Maximum der Abstrahlung zustande 
kommt. 


3. MeBanordnung 


Ein Schwebungssummer wurde kurzzeitig an ein 
Anregungssystem R12b (Neumann) gelegt, wo- 
durch impulsartig Wechselstrome mit konstanter 
Amplitude aber verschiedener Frequenz die Mitte 
des Geigensteges in Richtung der Saitenschwingun- 
gen erregten. Mit dem gleichen Schalter wurde ein 
Kathodenstrahloszillograph gesteuert, so daB dessen 
Leuchtpunkt gleichzeitig von links nach rechts ab- 
gelenkt wurde, wahrend auf das andere Plattenpaar 
die verstarkten Mikrophonspannungen 
wurden. 


gegeben 
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Zur Schallaufnahme diente ein in 15 cm Abstand 
unter 45° zur Geigenebene aufgestelltes Mikrophon. 
Die Saiten waren durch eingeschobene Gummistiicke 
untereinander abgedampft. 

Die angewandte Impulsmethode liefert die An- 
regungsfunktionen ohne Beriicksichtigung der Un- 
stabilitaten des normalen Bogenanstrichs. 
wahrt auBerdem den Vorteil, 
nahezu ausgeschaltet sind. 

In bekannter Weise [11] wurde zunachst mit 
Pegelschreiber kontinuierlich die Resonanzkurve der 
zu untersuchenden Geige (Bild 1a—d) aufgenom- 
men und darauf an einigen ausgezeichneten Punk- 
ten: Flanken, Maxima und Minima der Resonanz- 
kurve die Einschwingvorgange registriert. Zur Prii- 
fung gelangten zwei mafsiige moderne Violinen und 
zwei klanglich hervorragende Geigen von Stradi- 
varius. 


Sie ge- 
dafS Raumeinfltisse 


4. MeBergebnisse 
Die aufgezeichneten Einschwingvorgange (Bild 2) 
lassen sich nach gewissen Typen ordnen: 
1. Wenn f= fp ist, nimmt die Amplitude kontinuier- 
lich, wie oben angegeben, bis zu ihrem statio- 
naren Wert zu. 
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Wenn / wenig kleiner oder gré8er als fp ist, wo- 
bei bei fy vorwiegend eine Resonanzspitze vor- 
handen ist, wachst die Amplitude tiber ein Maxi- 
mum und darauffolgendes Minimum bis zu einem 
stationaren Wert, der kleiner als der des Maxi- 
mums ist. Die Erscheinung folgt, wie oben an- 
gegeben, aus der Uberlagerung der abklingenden 
Systemschwingung mit der stationaren. Bei ge- 
ringer Dampfung des angeregten Systems ist die 
Modulation besonders deutlich. 


. Mehrere aufeinanderfolgende Maxima und Mi- 


nima im Einschwingvorgang deuten auf die An- 
regung mehrerer schwach gedampfter benachbar- 
(elektrisches Analogon; 
gliedrige gekoppelte Schwingungskreise) . 


er Resonanzen mehr- 
Befindet sich f in einem Minimum zwischen zwei 


(fo </<fo), so ist die 
Anfangsamplitude am groBten. Der Einschwing- 


starkeren Resonanzen 


vorgang verlauft in dem Fall invers, die Schwin- 
groBten Anfangswert 
schwebungsartig bis zu einer relativ kleinen sta- 


gung klingt von einem 


tionaren Amplitude ab. Die Energie pendelt da- 
bei zwischen der einen und anderen Resonanz 
hin und her. 
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Bild 1. Resonanzkurven der untersuchten Geigen; 
(a) maBige Geige, bezeichnet mit PTB-1, 
(b) maRige Geige, bezeichnet mit PTB-2, 
(c) Stradivarius-Geige von 1694, 
(d) Stradivarius-Geige von 1721. 
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. Falls f in einem Frequenzgebiet liegt, in dem sich 
mehrere Resonanzen zu einem breiten Resonanz- 
bereich tiberlagern, so wird der stationare Am- 
plitudenwert in kiirzester Zeit (in wenigen Perio- 
den) erreicht. é 
Auf die in diesem Zusammenhang wichtigen Be- 

rechnungen von K. W. Wacner [9] uber die spek- 

trale Darstellung der rechteckformig eingeschalteten 

Sinus- bzw. Kosinusschwingung soll hier nur _ hin- 

gewiesen werden. Es lat sich daraus entnehmen, 

dafs im kontinuierlichen Spektrum eines solchen Vor- 
gangs die Amplituden A(w) der dicht benachbarten 

Teiltone mit wachsender Frequenz w bis zu dem 

Maximum der tatsachlich geschalteten Frequenz zu- 

nehmen 


Q 
ON Grae 


(2 geschaltete Frequenz). Dariiber werden sie wie- 
der kleiner. Ein solches Spektrum ist nach Betrag 
und Phase mit der Ubertragungsfunktion des Re- 
sonanzsystems zu bewerten, woraus das resultierende 
Spektrum und der zeitliche Verlauf des Ausgleichs- 
vorganges durch Laplace-Transformation abgeleitet 
werden kann. 

Zunachst wurden aus der Resonanzkurve der 
Geige einige Frequenzen f) der Resonanzmaxima er- 
mittelt. Sie sind in Tabelle I zusammengestellt. Aus 
den zugehorigen Einschwingvorgangen wurden die 
Dampfungen in Form des Dekrements 
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und der Dampfungskonstante 6=/d ausgewertet. 
(Der Verlustfaktor des E-Moduls, der in der Arbeit 
von Skuprzyk angefihrt ist, betragt 1 = d/zx.) 

Die gewonnenen Werte erfassen die mittlere innere 
Dampfung sowie die Strahlungsdampfung. Sie be- 
ziehen sich auf die Transversal- und Biegeschwingun- 
gen, welche die Platten der Geige ausfithren. Wah- 
rend die von Sxuprzyk [4] ermittelten Werte fir 
Klangholz zwischen d=0,015 bei tiefen und 0,080 
bei hohen Frequenzen liegen, ist das hier gemessene 
Verhalten weniger eindeutig. Ahnlich ist nur der 
Befund, da die Dampfungen der alten Geigen unter- 
halb derjenigen der neuen liegen. Die Tendenz des 
Dampfungsanstiegs mit wachsender Frequenz wird 
uber ungefahr 700 Hz sichtbar. Darunter tritt ein 
gegenlaufiges Verhalten auf; die Dampfung wird 
mit abnehmender Frequenz grofer. Es wird ver- 
mutet, daf dafiir die grofere Strahlungsdampfung 
infolge der besseren Abstrahlung bei weniger unter- 
teilten Schwingungszustanden verantwortlich zu 
machen ist. 

Es sei noch bemerkt, dafS sich die GroBe der 
Dampfung wahrend des Einschwingens andert. Bei 
kleiner Amplitude ist das Dekrement kleiner als bei 
groBerer Amplitude, die zeitliche Amplitudenzunahme 
erfolgt also nicht streng nach einer e-Funktion. 

Um mit den erwahnten Lautsprechermessungen 
von Hrtmpotp [6] vergleichbare Ergebnisse zu er- 
halten, wurden die Ausgleichsvorgange auf die oben 
beschriebene Weise ausgewertet. Die Ergebnisse sind 
in Bild 3 zusammengestellt. Die Einschwingzeit 
wurde dabei in Zahl] der Perioden angegeben. Sie ist 


Tabelle I. 
Dampfungswerte (Dekrement d und Diampfungsfaktor 0) der untersuchten Geigen. 


Frequenz l 
Hz - = AAA eee 


Dampfungen 
Strad. 1694 | Strad. 1721 


d 6 | d ) 


0,159 
0,093 
0,073 


0,288 


60,5 
102 
0,128 36 
0,063 36 
0,082 53 
0,051 


0,148 2 
0,053 3 
0,041 3 
0,049 wtf 
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mit max. 50 Perioden wesentlich linger als die von 
Lautsprechern, bei denen nur bis drei gemessen 
wurden. Zur frequenz- und phasentreuen Ubertra- 
gung eines groBen Frequenzbandes ist natirlich die 
Geige nicht geeignet. 


Penioden 


© | 


a =| 
GUO reas eee 4 6 10H? 


Frequenz ——> 


Bild 3. Dauer des Einschwingens in Perioden in Ab- 


hangigkeit von der Frequenz; 

(a) bei den Geigen PTB-1 und 2, 

(b) bei den Stradivari-Geigen von 1694 und 
1720 


Die Auswertung zeigt ferner, dais die Kinschwing- 
zeiten der Stradivari-Geigen deutlich kiirzer als die 
der neuen Instrumente sind, trotzdem man ent- 
sprechend der kleineren Dampfung der alten Geigen 
langere Einschwingzeiten bekommen miiBbte. Die 
Erscheinung kann dadurch erklart werden, dafi man 
annimmt, das schnellere Einschwingen komme durch 
das Zusammenwirken mehrerer bandfilterartig ge- 
koppelter Resonanzsysteme zustande. Diese sind 
wohl einzeln weniger stark gedampft, indessen fiihrt 
die Uberlagerung der beteiligten Resonanzen zu 
kiirzeren Einschwingzeiten, woraus sich allerdings 
keine e-Funktion ftir den Einschwingvorgang ergibt, 
wie es an sich bei Auswertung der Dampfung vor- 
ausgesetzt wurde. 

An den Stellen, an denen die Resonanzkurve ein 
Maximum zeigt, ist die Einschwingzeit am langsten. 
Man sieht, dafi derartige Maxima bei den neuen 
Geigen hoch sind und relativ weit auseinanderliegen. 
Bei den Stradivari-Geigen liegen solche Maxima 
dichter und sind weniger hoch. 

Diese Unterschiede sind typisch und diirften 
wesentlich zur groéBeren Klanggiite der Stradivari- 
Geigen beitragen. 

Bei 750 Hz zeigt die PTB-1-Geige das erwartete 
Verhalten: Hier wurde eine aufergewohnliche ge- 
ringe Dampfung gemessen, der entsprechend dem 
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normalen Zusammenhang eine besonders lange Ein- 
schwingzeit entspricht. Die Resonanzkurve zeigt an 
dieser Stelle eine besonders scharfe Resonanz mit 


.steilem Abfall nach tieferen Frequenzen. 


So ist es erklarlich, da} es fiir die klangliche Gite 
einer Geige giinstig ist, wenn mehrere schwach ge- 
dampfte Resonatoren benachbarter Frequenzlage vor- 
handen sind. An solchen Stellen ist die Ansprache 
rasch beendet und die Abstrahlung, d.h. der Wir- 
kungsgrad relativ gut. Bei einzelnen schmalen Re- 
sonanzen kann wohl der Wirkungsgrad auf dem 
Maximum noch ‘besser sein, indessen ist die An- 
sprache benachteiligt. Der Wirkungsgrad ist beson- 
ders klein, wenn die gespielte Frequenz auf ein 
Minimum der Resonanzkurve fallt. Sind die benach- 
barten Resonanzspitzen so weit voneinander entfernt, 
dai die auftretende Schwebungsfrequenz tiber etwa. 
20 Hz liegt, so treten schnelle Modulationen auf und 
der Klangeinsatz wirkt rauh. Das Vorhandensein 
mehrerer tief und breit eingeschnittener Minima 
kann daher allgemein als ungiinstig angesehen wer- 
den. Auch aus diesem Grund diirfte es sich empfeh- 
len, die Resonanzkurven in logarithmischem Ampli- 
tudenma8stab aufzuzeichnen, weil bei linearem tiber 
die Minimastellen nicht viel ausgesagt werden kann. 

Bei der Diskussion von Resonanzkurven ist daher 
auf folgende Punkte zu achten: 

a) Breite und Frequenzlage der Resonanzmaxima, 

b) Vorkommen von bandfilterartigen Resonanz- 
gebieten, 

c) Breite und Frequenzlage der Minima. 


Es kommt nattirlich selten vor, da die harmo- 
nisch liegenden Teiltone der Saitenschwingung alle 
auf Resonanzmaxima fallen. Meist werden nur 
einige Teilténe eines Klanges durch Resonanz ver- 
starkt. Solche erscheinen dann im _ Einschwing- 
vorgang relativ spat; andere Obertone, die auf Flan- 
ken, auf Resonanzgebiete oder Minima fallen, sind 
friiher aufgebaut. Durch diese verschiedene zeitliche 
Entwicklung der Klangkomponenten wechselt daher 
auch der subjektive Eindruck des Einschwingvorgan- 
ges von Ton zu Ton, wie auch im stationéren Vorgang 
die Klangfarbe von Ton zu Ton wechselt. Sind nun 
aber, wie oben bei b) angegeben, Resonanz ge biete 
vorhanden, so sind die Anderungen hinsichtlich des 
Klangaufbaues und des stationdren Vorganges ge- 
ringer, die Klangfarbe ist bei Grundtonanderung 
weniger sprunghaften Anderungen unterworfen. 
Solche Erscheinungen konnten bei den Stradivari- 
Geigen beobachtet werden. Wenig giinstig ist es, 
wenn mehrere Teilténe auf Minima der Resonanz- 
kurve fallen. Solche Klange wirken ausgesprochen 
matt. 


Auf einige Einzelheiten der Resonanzkurven soll 
noch hingewiesen werden: Die Stradivari-Geige von 
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1721 besitzt ein ausgesprochen starkes und breites 
Resonanzgebiet bei etwa 280 Hz. Mift man die 
Frequenzbreite dieser Resonanz (Hohlraum) bei den 
verschiedenen Geigen bei dem willkiirlich gewahlten 
Pegel von 30 dB, so ergeben sich folgende Grenzen 
bei PTB-1: 250 bis 270 Hz; bei PTB-2: 270 bis 
290 Hz; bei Stradivari 1721: 245 bis 310 Hz und 
bei Stradivari 1694: 265 bis 330 Hz, d.h., die so 
ermittelte Resonanzbreite ist also bei Stradivari mit 
65 Hz mehr als dreimal so gro als bei den moder- 
nen Instrumenten. Aber auch die Hohe der Spitzen- 
werte ist deutlich verschieden. Bei gleicher An- 
regungsenergie liegt diese bei den altitalienischen 
_ Geigen bei 46 und 41 dB gegeniiber 39 und 37 dB 
bei den Vergleichsgeigen. Die Grundtone auf der 
e°-Saite von Stradivari 1721 sind daher aufer- 
gewohnlich intensiv und sprechen in relativ kurzer 
Zeit an. Von Kinstlern wird diese Geige wegen ihrer 
Fille und Tragfahigkeit auBerordentlich geschatzt. 
Wie auch friher von H. Metner [2] angegeben 
wurde, kommt den Grundtonen einer Geige eine be- 
sondere Bedeutung zu. Auf die interessante tiefere 
Resonanz bei 150 Hz, welche durch das System 
Saitenhalter + 4 Saiten zustande kommt, sei beson- 
ders hingewiesen. 

Es fallt weiterhin auf, daB bei beiden Geigen der 
Frequenzbereich zwischen der Hohlraumresonanz 
und den Hauptmaxima auBergewohnlich gleichmabig 
und von mittlerer Intensitat ist. Minima treten auf- 
fallend wenig auf und sind weder breit noch tief 
eingeschnitten. So sprechen die Teiltone 1. und 
2. Ordnung in diesem Bereich schnell an und haben 
eine ausreichende Intensitat. Unterschiedlich ist der 
Frequenzbereich zwischen 500 und 1000 Hz aus- 
gebildet. Wahrend die Stradivari 1721 ein einziges 
starkes Maximum bei 580 Hz besitzt, hat die Stradi- 
vari 1694 deren zwei bei 580 und 660 Hz. Beson- 
ders bemerkenswert ist dabei, dafB bei derselben 
Geige vor zwei Jahren nur ein Hauptmaximum bei 
der tieferen Frequenz gemessen wurde. Leider lassen 
sich die Griinde fiir dies Verhalten nicht mehr re- 
konstruieren; es diirfte aber auf Grund der laufen- 
den Messung anderer Geigen feststehen, dai sich die 
Hohe der Resonanzen in Abhangigkeit vom Wetter 
verandert, Hinsichtlich des Einschwingens diirfte die 
Aufspaltung der Hauptresonanz nicht ungiinstig sein. 
Im Klange wirkt sich diese Eigenschaft in der Rich- 
tung aus, da die Geige etwas heller als die von 
1721 zu klingen scheint. In der gleichen Richtung 
wirkt das Resonanzgebiet von etwa 820 bis 1000 Hz 
der Stradivari von 1694, das bei der Geige von 
1721 zwar unten stirker aber oben etwas schwacher 
ausgebildet ist. Abgesehen davon ist aber bei beiden 
Geigen die Betonung des ,,o“- und _,,a“-Formanten 
sichtbar. Im Bereich zwischen 1000 und 2000 Hz 
fallen bei der Stradivari von 1721 zwei breite und 
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tiefe Minima bei 1130 und 1600 Hz auf. Trifft der 
2.Teilton auf das untere Minimum, so liegt der 
Grundton fast auf dem Hauptmaximum. Das be- 
deutet, daB die Lticke im Teiltonaufbau praktisch 
kompensiert ist. Der 3.Teilton fallt ebenfalls auf 
ein Maximum, der Klang wird also voll und kraftig 
wirtken. FaBt man 1130 Hz als 3.Teilton auf, so 
liegt der 2. auf der unteren Flanke eines Maximums, 
spricht also rasch an, der erste liegt auf der unteren 
Kante eines Resonanzgebietes, ist also ebenfalls 
schnell aufgebaut. Der 4.Teilton liegt auf der oberen 
Flanke eines Resonanzgebietes; fiir ihn gilt das- 
selbe. Man findet weiterhin bei dem 1600-Hz-Mini- 
mum, dal} es beim Auftreffen des 2. Teiltones wegen 
des starken Grundtons und seines raschen Aufbaus 
unwirksam ist. Wenn der 3. Teilton auf das Minimum 
fallt, ist der zweite Teilton relativ stark und der 
erste liegt auf der unteren Flanke der Hauptresonanz. 
Der 4. Teilton ist bei 2130 Hz ebenfalls stark. Auch 
wenn der 4.., 5., 6. usw. Teilton auf dieses Minimum 
fallt, sind soviel kompensierende Maxima bei den 
anderen harmonischen vorhanden, daf sich die Liicke 
nicht schadlich auswirkt. AuBerdem ist bekannt, da$ 
gute Geigen bei dieser Frequenz eine Minimumstelle 
besitzen sollen, damit die Gefahr des Niselns ver- 
mieden wird. Bei der Stradivari von 1694 sind die 
analogen Minima etwas verschoben und weniger 
breit (1150 und 1700 Hz). Wegen der starken und 
dichten Resonanzstellen bei tieferen Frequenzen ist 
aber vermieden, daf sich diese Liicken schadlich aus- 
wirken. 

Zahlt man vergleichsweise alle Resonanzmaxima 
der vier untersuchten Geigen uber 30 dB zusammen, 
so bekommt man bei PTB-1 und 2 je 11 Resonanzen, 
bei Stradivari 1721: 21 und Stradivari 1694: 17, 
wobei die obere Frequenzgrenze bei etwa 3500 Hz 
liegt. Der Ubertragungsbereich ist also entsprechend 
Bild 3 bei wertvollen Geigen mit starkeren und dich- 
ter liegenden Resonanzen besetzt. 


Die Partien tiber etwa 2000 Hz nach oben wirken 
sich besonders auf den sogenannten Glanz der Geige 
entsprechend. der Helligkeit des ,,e“- und_,,i“-For- 
manten aus. Bei der Stradivari 1694 sind auffallend 
starke, relativ breite Resonanzen zwischen 1900 und 
4000 Hz vorhanden, dariiber hinaus aber praktisch 
nur noch geringfiigige. Diese Violine wirkt daher 
besonders hell und strahlend, doch nie scharf. Als 
noch besser wird allerdings die Stradivari von 1721 
bewertet, deren Maxima im fraglichen Bereich weni- 
ger dicht liegen und weniger hoch sind. Mit grofen 
Liicken ist besonders der Bereich zwischen 2000 und 
3000 Hz 3000 und 
4000 Hz starkere Spitzen zu erkennen sind. Diese 


versehen, wahrend zwischen 


Eigenschaft bewirkt offenbar, daB der Stradivari 
1721 von Kennern der besonders ,,milde Glanz“ 
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nachgertihmt wird, welcher den Spitzeninstrumenten 
neben ihrer Sonoritat besonders eigen ist. 
Entsprechende Auswertungen lassen sich an mo- 


dernen Vergleichsinstrumenten anstellen. Es sei hier - 


lediglich bei PTB-1, welche als wenig gut anzusehen 
ist, auf das starke Naseln aufmerksam gemacht, das 
durch das ausgedehnte Resonanzgebiet um 1500 Hz 
und eine entsprechend ungiinstige Intensitatsvertei- 
lung der anderen Teiltone hervorgerufen wird. Wenn 
namlich der 3. Teilton auf das Maximum bei 1500 Hz 
fallt, liegt der 2. Teilton bei 1000 Hz, ist nach der 
Resonanzkurve demnach schwach, der 1. ist, weil er 
nahezu auf dem Hauptmaixmum liegt, stark, der 
4.. ist schwach. Bei solcher Zusammensetzung (ohne 
starkere Oktave) ist die nasale Wirkung besonders 
gro}. Wenn der 2.Teilton auf 1500 Hz fallt, ist 
keine ungiinstige Wirkung zu befiirchten, ebenso 
nicht bei dem 4., 5. und 6., bei letzterem sicher nicht, 
weil der 2. und der 3. Teilton auf Resonanzen der 
Kurve fallen. 

Bei der Anregung unterhalb der Hohlraumreso- 
nanz (Bild 2) fallt auf, ein wie starker Klirrfaktor 
[12] in Erscheinung tritt. Er ist dadurch erklarlich, 
daf$ im Impulsvorgang die hoheren Harmonischen 
angestoen werden und sich dadurch, dafi in diesen 
Frequenzlagen eine erheblich grofere Resonanz vor- 
handen ist als beim Grundton, stark bemerkbar 
machen. 


5. SchluB 


Durch die dargestellten Untersuchungen lassen 
sich die im Eingang erwahnten Beobachtungen von 
H. von Hermuortz erklaren. Die Unstabilitaten des 
Bogenanstrichs verursachen sowohl Amplituden- wie 
Frequenzschwankungen der Saitenschwingung. Be- 
sitzt der Resonanzkorper nur wenig Resonanzen re- 
lativ geringer Dampfung und dazwischen mehrere 
Minima, so werden auch geringe Frequenzschwan- 
kungen in erhebliche Amplitudenschwankungen um- 
gesetzt. Das System der Saite arbeitet zeitweise auf 
Resonanz, zeitweise auf ein Minimum, woraus folgt, 
da} die Phasenlage der Schwingung immer wieder 
wechselt, es demnach schwierig ist. eine stabile 
Schwingung zu erzeugen. Wenn dagegen der Re- 
sonanzkorper eine grofere Zahl von dicht benach- 
barten Resonanzen ohne starkere Minima besitzt, so 
wirken sich die Unstabilitaten des Bogenanstrichs 
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wesentlich geringer aus. Bei geringen Veranderungen 
der Grundfrequenz andern. sich die Resonanzbedin- 
gungen nur wenig, es ist also modglich, auch bei 
Frequenzanderung stabile Schwingungen aufrecht- 
zuerhalten. 

Mit Hilfe der dargestellten Methode 1a8t sich also 
die Gite einer Geige von dem Gesichtspunkt der 
Ansprechfahigkeit, also derjenigen Eigenschaft, 
welche fiir den Spieler von besonderer Wichtigkeit 
ist, begrunden. 


Der Verfasser dankt der Deutschen Forschungs- 
gemeinschaft fiir die Bereitstellung von Mitteln fiir 
Geigenuntersuchungen, Herrn Prof. Dr.-Ing. Lurz, 
Technische Hochschule, Braunschweig, fiir die Ver- 
mittlung der beiden Stradivari-Geigen. 

Den Herren Katkxe und Simon danke ich fiir 
ihre Hilfe bei den Messungen. 

(Eingegangen am 6. Juni 1957.) 
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OBLIQUE AIR-TO-GROUND SOUND PROPAGATION 


OVER BUILDINGS 
by P. H. Parkin and W. E. Scuotzs 


Department of Scientific and Industrial Research, Building Research Station, Garston, 
Watford, Herts 


Summary 


The propagation sideways to the ground of the sound emitted by helicopters has been 
investigated in a built-up area. Over the horizontal distance covered (about 920m 
(3.000 ft) it was found (a) that the attenuation (i.e. all losses in excess of the inverse 
square reduction) at frequencies up to 1000 c/s was small and independent of the weather 
for helicopter heights of greater than about 100m (300 ft); and (b) that at lower heights, 
the attenuation increased considerably. 


Sommaire 

On a étudié dans une zone urbaine, comment le son émis par un hélicoptére se propa- 
geait obliquement vers le sol. La portée horizontale était de l’ordre de 920m; lorsque 
Valtitude de Vhélicoptére est supérieure 4 100m environ, l’affaiblissement qui vient 
s’ajouter a la variation en raison de l’inverse du carré de la distance est faible pour des 
fréquences <_ 1000 Hz, et ne dépend pas des conditions atmosphériques. Pour de plus 
faibles altitudes, Paffaiblissement augmente trés notablement. 


Zusammenfassung 

Es wurde in bebautem Gelande die Ausbreitung der Gerdusche von Hubschraubern iiber 
dem Erdboden untersucht. In dem beobachteten Bereich (Horizontalentfernung bis ca. 
920 m) ergab sich: a) Die Dampfung (d.h. alle Verluste, die iiber das quadratische Ab- 
nahmegesetz hinausgehen) ist klein fiir Frequenzen bis 1 kHz und unabhangig vom Wetter 
fiir Flughéhen grofer als etwa 100m, b) bei geringeren Flughdhen steigt die Dampfung 


betrachtlich an. 


1. Introduction 


A previous investigation [1] had shown that 
sound propagated vertically downwards from an air- 
craft in flight was attenuated only by molecular 
absorption (apart of course from the inverse square 
reduction). The investigation reported now was 
concerned with the oblique propagation of sound, 
i.e. the sound propagated to the ground sideways 
from an aircraft in flight. It might be expected that 
such propagation would be affected by wind and 
temperature gradients, by the presence of the ground 
and by the molecular absorption. In what follows 
the effects, if any, of all these factors are combined 
under the term ,,attenuation™. 

The aircraft used were Westland’s type S.55 
helicopters, fitted with exhaust silencers. They flew 
along the line of the River Thames in London, and 
the sound pressure levels were measured by ten 
microphones at roof top level (21.5m (70 ft) to 
36.5 m (120 ft) above ground) situated at roughly 
equal intervals along a line at right angles to the 
line of flight. Fig. 1 shows the plan of the area and 
Fig. 2 the section, with the microphone positions and 
line of flight of the helicopters indicated. 

The investigation was in two parts. The first part 
consisted of measurements on one calm day with.a 


helicopter flying at various heights of between 46 m 
(150 ft) and 305 m (1000 ft) relative to the mean 
height of the microphones. The objects were to 
determine the directivity of the helicopter as a noise 
source and to investigate the oblique attenuation as 
a function of height. The second part was spread 
over 28 days throughout the period December of 
one year to May of the following year, using iden- 
tical helicopters which were engaged in regular 
passenger service between the centre of London and 
London Airport. These helicopters flew at heights 
of between 105m (340 ft) and 230m (750 ft) all 
relative to the microphones, with an average value 
of 160 m (520 ft). The object was to investigate the 
effect of the weather on the attenuation. 


2. Technique 


The ten microphones on stands 1.5 m (5 ft) high 
were placed on the highest accessible positions on 
the buildings. (All the microphones except the first 
were on flat roofs; the first had behind it a pitched 
roof, see Fig. 2). They were protected by plastic 
covers and windshields, which had no effect on their 
sensitivities at frequencies up to 1 000 c/s, and little 
effect between 1 000 and 2000 c/s. After amplifica- 
tion the signals were fed via telephone land-lines to 
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Fig. 1. Plan of microphone layout. 
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Fig. 2. Section of microphone layout. 


a central control room where they were recorded on 
ten tape-recorder channels. An overall electrical cali- 
bration was made every day. 

The heliocopters flew nominally along the centre 
of the river, and their position in space and in time 
relative to the recordings was determined by a 
calibrated camera. In order to avoid interference 
with the recordings at the critical moments when a 
helicopter was opposite the line of microphones, the 
indicating marks on the tapes (a short circuit of 
microphone signals) were delayed automatically 
until 10 seconds after the photograph. 

The recordings were replayed through octave 
band-pass filters onto a level recorder set at a writ- 
ing speed of 50 dB per second. From the short 
circuiting marks and by allowing for the time taken 
for the sound to travel, a point on the trace could be 
determined corresponding to the time of arrival at 
the microphone of the sound emitted by the heli- 
copter when it was at right angles to the line of 
microphones. The level recorder traces fluctuated 
rapidly by about +4 dB, and the sound pressure 
level was taken from the value of a smooth line 
drawn through these fluctuations averaged over an 
interval of about +1 second. 

In the lowest octave (37 — 75 c/s) the back-ground 
noise at most of the microphone positions was as 
great or greater than the helicopter sound, so no 
measurements were possible in this band. The highest 
octave measured was 1 200 —2 400 c/s; the plastic 
covers prevented measurements at higher frequencies, 
and in any case it was expected that the molecular 
absorption would predominate. 


3. Analysis of results 


3.1. Directivity of source 


The directivity pattern of the helicopter noise 
source in the plane at right angles to the line of 


flight will obviously affect the measured attenuation, 
and although the two effects — directivity and true 
attenuation — can not be completely separated it is 
thought that the technique here described was suf- 
ficiently accurate. 


On a day when the weather was calm the heli- 
copter flew at heights (relative to the microphones) 
of between 46m (150 ft) and 305m (1000 ft). 
The relative sound pressure levels, after correction 
for the inverse square law, at the first three micro- 
phone positions were plotted as a function of ele- 
vation of the helicopter with respect to each micro- 
phone, the two directions of flight being plotted 
separately. For example, Fig.3a—c shows the 
results for the three microphones for the 300 to 
600 c/s octave. There is of course no reason why 
the directivity should be a linear function of angle, 
but the results do not justify any more accurate © 
evaluation, so the directivities are summarised in 
Table I in terms of decibels per degree. 


It is seen that at all frequencies the results for 
microphones 2 and 3 were similar but were from 
two to four times as great as the results for micro- 
phone 1. Now microphone 1 was on a ledge with a 
large roof sloping at about 60° rising above and 
behind it, and it is suspected that the reflections 
from this roof were responsible for the differences. 
Because of this the results from microphone 1 have 
been omitted from the directivity calculations, and 
the results from microphones 2 and 3 have been 
averaged. 


For microphones 2 and 3 the differences between 
the directivities on the port and starboard sides were 
probably due to the exhaust being on one side. 
However, the maximum difference between the two 
sides is 0.05 dB per degree (in the 1 200 — 2 400 e/s 
octave). It is more convenient when analysing the 
other results not to have to distinguish between the 
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Table I. Directivity of source; dB per degree. 


Microphone 1 


Microphone 2 
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| 0.02 0.08 0.08 0.08 | 0.06 
| oo8 | o16 | 008 | o14 | O12 


not necessarily be a linear function of horizontal 
distance, the results obtained in these experiments 
justified the assumption made in the analyses that 
the attenuation could be treated as a linear function, 
at least over the total distance of about 920 m 
(3 000 ft) covered by the line of microphones. 
Accordingly, the attenuation (in dB per km (dB per 
1000 ft)) has been taken from regression lines 
calculated for each day. 

For example, for Fig. 4 the slope of the regression 
line was 5 dB per km (1.5 dB per 1 000 ft). 

The average attenuations for all the 28 days, and 
the standard deviations and the ranges for the five 
octave bands are given in Table II, and graphically 
in Fig. 5. 


Octave band | 
e/s | 
| 
75— 150 | 0.06 
150— 300) 0.04 
300— 600 0.06 
600—1200, 0.02 
1200—2400| 0.06 
é 3 
3 Starboard 
g 
Fig. 3. 
Relative sound pres- 
sure levels as a func- 
L tion of elevation of 
Starboard helicopter (300 to 
600 ¢/s) ; (a) Micro- 
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Fig. 4. Typical plot of excess attenuation. 


two directions of flight. The error in the subsequent 
calculations of attenuation introduced by the dif- 
ference in the directivity between port and starboard 
sides will be greatest when the helicopter was at its 
geatest height, but even_then the error amounts to 
less than +1.5 dB per km (+0.5 dB per 1 000 ft) 
for this 1200—2400c/s octave, and for other 
octaves and lower heights the error is smaller still. 
Thus, in the further analyses no distinction is made 
between the two directions of flight. 


3.2. Variation of attenuation with weather 

The average number of measured flights per day 
for the 28 days throughout the six month period was 
5.5. The results for each day have been averaged 
because the weather conditions, on these days, did 
not change appreciably between flights. 


Fig. 4 shows a typical plot of the attenuations at 
the ten microphone positions for the 600 to 1 200 c/s 
octave on one of the days. The values have been 
corrected for the inverse square law and for the 
directivity of the source. While any attenuation due 
to the ground surface, or to any other cause, will 
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Fig. 5. Average attenuation over 28 days. 
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Table II. 
Attenuations in dB per km (dB per 1000 ft). 
Octave band | Standard 
c/s | verre deviation Penge 
75— 150 | — 0.8 (— 0.25) a3 Gen): 4 — 6.5 (2 
| | tos S32 tones 
150— 300 0 (0) | g(elye 38 (sc 
| to + 6.5 to + 2) 
| 
300— 600 + 0.7 (+ 0.2) | Sry athe ae S35 oar 
to + 8.5 to + 2.5) 
600—1200 gee ee ahh) | +3(+1) | an; (1.5 
| to + 8.5 to + 2.5) 
1200—2400 + 11.5 (+ 3.5) 43 (Ae eae 6.5 (+ 2 to 
to 13" 


3.3. Variation of attenuation with height of heli- 
copter 


On a calm day the helicopter flew at heights of 
between 46m (150 ft) and 305m (1000 ft) (with 
respect to the mean height of the microphones), four 
flights at each height. The results for each height 
were averaged, and analysed as previously described 
(except that two of the ten microphones were not 
used). The attenuations as a function of helicopter 
height are plotted in Fig. 6. 


4. Discussion of results 


Considering first the results for the 28 days, 
although the weather was not necessarily an exact 
sample of all possible weather conditions, it was 
varied and reasonably typical of this locality. For 
example, the turbulence varied between ,,light“ and 
»severe’ the temperature (at ground level) between 
= 45° Ci (24°F) and-20:57€ “(71> FR), «andthe 
relative humidity (at ground level) between 30 per 
cent and 85 per cent. The maximum wind speed 
within +20° of the direction of the sound pro- 
pagation was 16km per hour (10 mph) at the 
surface and 32km per hour (20 mph) at a height 
of 305m (1000 ft). The maximum wind speed in 
the opposite direction was 32 km per hour (20 mph) 
at the surface and 64km per hour (40 mph) at a 
height of 305m (1000 ft). Nevertheless, in spite 
of this wide variation of weather conditions, there 
were only slight changes in the attenuation, and 
even these it was not possible to correlate with the 
weather. It appears therefore that up to about 
1000 c/s and at these helicopter heights no signifi- 
cant attenuation is caused by wind and temperature 
gradients or by the presence of the ground. 

In the highest octave measured (1200 to 
2400 c/s) there was some appreciable attenuation. 
This was probably due to the molecular absorption 
of the air, although it could not be correlated with 
temperature and humidity. However, an octave band 
is rather too wide for such a correlation. 


It might be added that even if the assumption made 
in the directivity and attenuation calculations were 
not valid, it would still be true that the weather did 
not vary the attenuation. This is shown by the fact 
that the absolute sound pressure levels at the micro- 
phones did not vary appreciably from day to day. 
For example, the average value for all the 28 days 
in the 300 to 600 c/s octave at microphone 8 (610 m 
(2 000 ft) distance from the source) was 61 dB, with 
a standard deviation of +2 dB. 

All measurements were made in day-light hours, 
but the only factor likely to be consistently dif- 
ferent during darkness is the temperature gradient. 
But in any case temperature inversions which might 
occur at night should in theory tend to reduce the 
attenuation. 

It is seen from Fig. 6 that the attenuation in all 
octaves increased rapidly when the helicopter height 
was less than about 100 m (300 ft). These particular 
measurements were made on one day when the 
weather was calm. It may be that at these lower 
heights the weather will affect the propagation. 
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ETALONNAGES PAR RECIPROCITE EN CHAMP LIBRE 


DE TRANSDUCTEURS SOUS-MARINS* 
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Sommaire 


La méthode classique de réciprocité en champ libre permet 1’étalonnage de transducteurs 
linéaires sous marins pour les fréquences supérieures 4 5 kHz. Dans le présent article, 
Yauteur rappelle briévement le principe de cette méthode et insiste au passage sur les 
causes d’erreurs. Une bréve description est donnée de J’installation et des mesures faites au 
laboratoire de la Marine Nationale au Brusc. 


Zusammenfassung 


Lineare Wandler fiir Wasserschall kénnen bei Frequenzen oberhalb 5kHz nach der 
klassischen Reziprozitétsmethode im Freifeld geeicht werden. Der Verfasser gibt zundchst 
eine kurze Darstellung dieser Eichmethode und behandelt dann eingehend die einzelnen 
Fehlerquellen. AbschlieSend wird die MeBanordnung im Laboratorium der Marine Natio- 
nale in Brusc beschrieben und von den Messungen berichtet. 


Summary 


The reciprocity method used in a free field can be used to calibrate linear underwater 
transducers at frequencies above 5 ke/s. The author surveys briefly the principle of this 
method, pointing out sources of error. He then describes the installation and measurements 


carried out at the laboratory Brusc of the Marine Nationale. 


1. Objet 


Pour contr6ler ses installations d’acoustique sous- 
marine et pour mettre en ceuvre sur des bases so- 
lides un programme étendu de mesures de bruit, la 
Marine Nationale Francaise a eu besoin de trans- 
ducteurs électroacoustiques sous-marins étalonnés 
avec précision. Le Laboratoire de Détection sous- 
Marine du Brusc a donc créé et développé une in- 
stallation d’étalonnage répondant le mieux possible 
a ces besoins. 


2. Généralités 


La méthode d’étalonnage utilisée pour les fré- 
quences supérieures a 5 kHz est celle de réciprocité 
en champ libre. Bien qu’elle soit absolument clas- 
sique, jen rappellerai briévement le principe pour 
insister, au passage, sur les causes d’erreur. 

La méthode permet, on le sait, d’étalonner simul- 
tanément deux transducteurs que nous appellerons 
ae cekel se 

Comme cas particulier, il est possible d’étalonner 
simultanément un transducteur et un récepteur non 
reversible. 

L’application de la méthode suppose: 

— que les deux transducteurs sont linéaires, 
— que les conditions de champ libre sont remplies. 


* Communication présentée a Ja réunion du G.A.L.F. 
tenue 4 Milan du 9 au 12 octobre 1956. 


3. Linéarité 


La condition de linéarité postule simplement 1’ab- 
sence de phénoménes de saturation dans le trans- 
ducteur qui émet; ceci suppose qu’on ne lui ap- 
plique que des puissances faibles. Le critére de la 
linéarité est la proportionnalité entre la pression 
acoustique en un point quelconque du champ pro- 
duit par l’émetteur, et le courant qui l’alimente. 

Dans un transducteur qui regoit, la puissance 
transférée est extrémement faible. La tension recue 
est toujours proportionnelle a la pression acous- 
tique. Contréler le domaine de linéarité d’un émet- 
teur consiste donc simplement a rechercher dans 
quelle limite la tension aux bornes d’un récepteur 
placé de facon quelconque dans le champ, reste pro- 
portionnelle au courant dans |’émetteur. 

Il est entendu que les étalonnages dont nous par- 
lons ne sont valables que dans le domaine de linéa- 
rité commun a T, et T,. 


4. Champ libre 
Les conditions de champ libre justifient une petite 
discussion. 
4.1.0. Définissons la sensibilité d’un récepteur élec- 
tro-acoustique comme un quotient dont le 
numérateur est l’amplitude de la tension recue a 
circuit ouvert, et le dénominateur l’amplitude de la 
pression acoustique a l’endroit ow est placé le récep- 
teur. 


104. 


Une telle définition, et en particulier l’expression 
«a l’endroit ou» n’est claire et suffisante que pour 
une sonde, c’est-a-dire pour un récepteur trés petit 
devant la longueur d’onde et méme théoriquement 
ponctuel; il y a deux raisons a cela: 


— une sonde étant petite devant la longueur d’onde, 
Vamplitude de la pression acoustique est pratique- 
ment constante dans toute région du champ de 
dimensions analogues; ceci est vrai quelle que 
soit la forme du champ sonore et des systémes 
d’ondes stationnaires qui peuvent s’y développer. 


— l’introduction de la sonde dans le champ ne pro- 
voquant aucune diffraction, la pression effective- 
ment appliquée a la surface sensible de la sonde 
est la méme partout, en phase et module, et elle 
est égale a sa valeur avant introduction de la 
sonde. 


Il en résulte que le mode de vibration de la sur- 
face sensible n’est influencé ni par la forme du 
champ sonore, ni par l’orientation de la sonde; 
d’autre part, il n’est pas nécessaire, dans la défini- 
tion précédente, de préciser si le dénominateur de 
la sensibilité est la pression avant ou aprés l’intro- 
duction de la sonde. 


4.1.1. Cette simplicité disparait dans le cas d’un 

récepteur de dimensions non négligeables, 
voire méme grand devant la longueur d’onde. II 
faut préciser que le dénominateur de la définition 
est la pression en un point P du champ, avant lin- 
troduction du récepteur, puisque c’est bien cette 
pression préexistante qui nous intéresse quand nous 
utilisons la sensibilité. D’autre part, il faut préciser 
quel est le point lié au récepteur — ou point de 
référence — 
point P en plagant le récepteur dans le champ. Ce 
choix du point de référence influe comme I orienta- 
tion, sur la diffraction du champ préexistant, donc 
sur la valeur des pressions sur la surface sensible, 


que nous ferons coincider avec le 


donc sur la valeur de la tension aux bornes 4 circuit 
ouvert. 


4.1.2. Une sensibilité ne peut donc étre définie que 

dans un champ préexistant donné, pour un 
point de référence donné et une orientation donnée. 
C’est done une notion complexe et peu utilisable 
dans sa généralité. Lorsqu’on parle de sensibilité 
sans préciser, on entend la sensibilité en ondes pla- 
nes ce qui définit le champ préexistant. Le champ 
d’ondes planes a l’avantage en outre d’éliminer le 
choix du point de référence; dans un champ d’ondes 
planes idéal, l’amplitude de la pression est la méme 
en tout point et le dénominateur de la sensibilité se 
trouve alors trés bien défini quelle que soit la facon 
dont le récepteur est mis en place. La diffraction 
autour du récepteur est invariante par translation. 
Elle ne dépend que de Vorientation du récepteur 
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dans le champ c’est a dire de la directivité en récep- 
tion, Cette expression implique automatiquement 
quwil s’agit d’ondes planes. 


4.1.3. En pratique un champ d’ondes planes est une 


petite portion de champ sphérique. Envisa- 
geons, par conséquent, de placer notre récepteur 
dans un champ sphérique ot l’amplitude de la pres- 
sion est inversement proportionnelle au rayon du 
champ. A la région destinée a étre occupée par le 
récepteur correspondent deux valeurs extrémes de 
rayon. L’écart relatif entre ces deux valeurs est pré- 
cisément l’erreur relative maxima que l’on peut 
commettre sur le dénominateur de la sensibilité, en 
disant qu’il est égal a la pression «a l’endroit ot» 
l’on place le récepteur. Si cette erreur est faible de- 
vant la précision désirée, c’est a dire si les dimen- 
sions radiales du récepteur sont suffisamment petites. 
devant le rayon du champ, il n’y a pas lieu de choisir 
un point de référence sur le récepteur. Cette con- 
dition donne, en quelque sorte, le rayon minimum 
de champ sphérique qui dispense du choix d’un 
point de référence sans erreur appréciable sur le dé- 
nominateur de la sensibilité. 


4.1.4. Par ailleurs si on considére la dimension du 

récepteur normale au rayon du champ, son 
quotient par le rayon est une quantité petite; la dif- 
fraction en ondes sphériques et la diffraction en 
ondes planes ne différent que du second ordre par 
rapport a cette quantité. L’écart sur le numérateur 
de la sensibilité, c’est-a-dire sur la tension recue a 
circuit ouvert, est donc aussi du second ordre. Cette 
condition fournit donc le rayon minimum de champ 
sphérique qui donne une erreur négligeable sur le 
numerateur de la sensibilité. 


4.1.5. La plus sévére des deux conditions précéden- 

tes permet de définir le rayon minimum du 
champ d’ondes sphériques presque planes dans le- 
quel la sensibilité peut étre considérée a la précision 
voulue, comme égale a la sensibilité en ondes planes. 
Nous appellerons M cette sensibilité qui sera celle 
des étalonnages. 


4.1.6. En résumé, soit un champ sphérique qu’on 

peut considérer comme produit par un point 
source 2 (Fig. 1); soit une région R_ susceptible 
de contenir le transducteur T et remplissant dans 
ce champ les conditions de la «presque planéité» 
que nous venons de définir en fonction de la pré- 
cision désirée. 

Soit 7 amplitude de la pression dans la région R 
en absence de T, et soit e la tension recue a circuit 
ouvert par T lorsqu’il est mis en place. La valeur 
de M est par définition e/z. 

M=e/n. (1) 


La pression 2 qui est inversement proportionnelle 


a 
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Fig. 1. Champ d’un point source. 


au rayon peut s’exprimer en fonction du débit g 
du point source 2, c’est-a-dire de l’expression 
§§ wds ow w est la vitesse du fluide et ou l’intégrale 
est prise sur une sphere infiniment petite entourant 
le point P. Cette relation résulte immédiatement des 
équations du champ sphérique et s’écrit: 


pag tt (2) 


© désigne la densité du fluide, f/ la fréquence. 


4.2.0. Précisons a présent la définition de la ré- 
ponse en émission. 

C’est un quotient dont le numérateur est l’am- 
plitude de la pression acoustique en un point A 
du champ produit par le transducteur; le dénomina- 
teur est l’amplitude du courant 7 qui l|’alimente. 
Nous admettons que nous sommes dans le domaine 
de linéarité et que par conséquent cette quantité ne 
dépend pas de 7. 


4.2.1. Pour un émetteur trés petit devant la lon- 

gueur d’onde, omnidirectionnel par consé- 
quent, et se rapprochant du concept de point source, 
la définition précédente est telle que dans un milieu 
illimité la réponse en un point donné fournit la ré- 
ponse en tout autre point. Pour un émetteur directif 
au contraire, et méme en milieu illimité, la réponse 
dépend du choix du point A, et plus précisément des 
trois coordonnées de A dans un triédre lié a |’émet- 
teur. Prés de |’émetteur, dans ce qu’on peut appeler 
sa zone d’interférence, la réponse est une fonction 
en général trés complexe de la position de A. Si le 
point A s’éloigne dans une direction définie par 
deux angles groupés sous la notation @, l’amplitude 
de la pression tend a devenirinversement propor- 
tionnelle a la distance. On peut définir la distance 
minima pour laquelle la loi est vérifiée a tant pour 
cent pres. On peut donc, autour de |’émetteur, et a 
Vextérieur d’une zone d’interférence d’autant plus 
étendue que la précision exigée est plus grande, deé- 
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finir une zone dite de champ libre dans laquelle le 
champ est en 1/r a tant pour cent prés dans toutes 
les directions. 
4.2.2. En principe, cette zone de champ libre dé- 
pend du choix de l’origine du triédre lié a 
l’émetteur, c’est-a-dire du choix d’un point de ré- 
férence. Si la condition précédente suffit a rendre 
négligeable les écarts relatifs de distance entre le 
point A et un point quelconque de |’émetteur, il n’y 
a pas lieu de préciser le choix du point de référence. 
S’il n’en est pas ainsi, il faut dilater la zone d’inter- 
férence jusqu’a ce que les écarts relatifs soient de- 
venus négligeables. 

De toute fagon, la plus sévére des deux conditions 
définit le champ libre a tant pour cent prés dans 
lequel la distance r d’un point A peut étre comptée 
sans erreur appréciable a partir d’un point quel- 
conque de |’émetteur. 


4.2.3. Donc, en champ libre, la pression acoustique 

est le quotient, par la distance r d’une fonc- 
tion de la direction c’est-a-dire de 9. II est naturel 
de rapporter la valeur de la pression a sa valeur 
pour r=1 méme si r=1 ne correspond pas phy- 
siquement a un point en champ libre. Nous appelons 
réponse en champ libre S, le quotient par le courant 
i de la valeur de la pression pour r=1. Celle-ci est 
donc égale par définition au produit Sz et la pres- 
sion en champ libre a la distance r est égale par 
définition a: 


p=Sir. (3) 


La réponse S est évidemment fonction de @ et 
cette fonction est la directivité en émission du trans- 
ducteur. Cette expression contient implicitement la 
notion de champ libre. 

Une mesure de directivité en émission n’est cor- 
recte que si le champ est en 1/r. Il en est de méme 
pour les étalonnages destinés 4 mesurer S. L’écart 
par rapport a cette loi représente |’erreur de champ 
libre sur S et il est important de Ja connaitre pour 
apprécier la précision de la mesure de S, 

On peut dire en somme que, dans la zone de 
champ libre, l’émetteur se comporte comme un point 
source situé dans la région qu’il occupe, en un point 
qu’on peut négliger de préciser et dont la puissance 
varierait avec la direction suivant la méme loi que S. 


5. Conséquences de la réciprocité 


5.1. Entre la sensibilité M en ondes planes ou 

presque planes et la réponse S en champ libre, 
telles qu’elles ont été définies et dans les limites 
d’application ot leur définition est correcte — il 
existe une relation simple qui résulte du principe 
de réciprocité. On sait déja dans la pratique que la 
directivité en réception et en émission correctement 
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mesurées sont identiques; ceci démontre que le rap- 
port M/S ne dépend pas de la direction. On peut 
de plus montrer qu’il ne dépend pas du transducteur 
considéré. 

Dans un quadripéle électrique la réciprocité se 
traduit de la fagon suivante: 

Injectons un courant d’un cété du quadripdle; 
nous recueillons une certaine tension a circuit ouvert 
de l’autre cé6té; le rapport de cette tension au cou- 
rant injecté est une impédance de transfert. Nous 
pouvons de la méme facon mesurer l’impédance de 
transfert dans l’autre sens. La réciprocité est l’éga- 
lité de ces deux impédances de transfert. Elle se tra- 
duit par une relation entre les 4 coefficients des 
équations de quadripole. 

5.2. Prenons un transducteur T immergé dans 

un milieu illimité et considérons un point P 
de ce milieu, en champ libre. Aux bornes de T peu- 
vent apparaitre les variables électriques e tension, 
et ¢ courant. Au point P peuvent apparaitre les 
variables acoustiques homologues, pression p, et dé- 
bit g au cas ow il existerait une source en P. Entre 
ces 4 variables, et dans les limites de linéarité du 
transducteur, on peut écrire des relations linéaires 
et homogénes analogues 4a celles du quadripéle. Les 
bornes de T et le point P forment un quadripdle 
électro-acoustique auquel nous appliquons le prin- 
cipe de réciprocité (Fig. 2). 
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Fig. 2. Relations de réciprocité. 


Injectons un courant z dans le transducteur. Nous 
créons en P une pression «a circuit acoustique ou- 
vert» puisqu’il n’y a pas de débit. Le point P étant 
en champ libre, cette pression vaut p=Si/r (3) par 
définition de S. L’homologue de l’impédance de 
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transfert est 
pji=S/r. (4) 


Supposons a présent qu’un point source en P, de 
débit g, crée un champ acoustique et que le trans- 
ducteur T recoive a circuit ouvert une tension e. 
L’impédance de transfert dans ce sens est e/q. 

Or e est égal au produit Ma ou M est bien la 
sensibilité en ondes presques planes de T, et ou 2 
est la pression créée par le point source a l’emplace- 
ment de T et en l’absence de celui-ci. 

On a vu que a=q(of/2r) (2). Il résulte que 
Vimpédance de transfert 


~ elg=M(of/2r). (5) 


En égalant les deux impédances de transfert on 
trouve que: 


M/S = 2/o f (6) 


cest-a-dire que le rapport M/S ne dépend ni du 
transducteur, ni de la direction, mais seulement de 
la densité du milieu et de la fréquence. On pose en 
général ce rapport égal a J 


2/of=J. (7) 


Si l’on exprime M en volt par barye et S en barye 
par amp. a 1 métre, le coefficient J vaut, dans eau 
USO ol: 

f=210 37 


(f en Hz). (8) 


5.3. Il est évidemment important de se demander 

dans quelle mesure on a le droit de considérer 
une transformation électro-acoustique comme réci- 
proque. 

En vertu du principe de Hetmotz, le probléme 
ne se pose en réalité que pour les transformations 
électromécaniques. Cette question sort de notre sujet. 
Il ne semble pas que les conditions suffisantes mini- 
mum de réciprocité soient encore clairement établies. 
Par ailleurs des auteurs comme MacMr1tan ont pu 
montrer que des systémes mixtes couplant mécanique- 
ment deux transducteurs de type différent peuvent 
mettre le principe en échec!. Mais sa validité n’est 
pas contestée pour des transducteurs simples par 
exemple purement magnétostrictifs ou purement 
piézoélectriques dans la mesure ou précisément ils 
n’utilisent qu’un procédé de transfert de |’énergie. 

Nous nous en tiendrons a ces cas simples qui seuls 
nous intéressent pratiquement. 

Il est facile d’ailleurs de vérifier expérimentale- 
ment deux conséquences nécessaires de la récipro- 
cité. La premiére, bien connue, est lidentité des 
directivités en émission et en réception. La seconde 
est la suivante; si |’on place dans l’eau deux trans- 
ducteurs réciproques T, et T,, le quadripéle élec- 
trique formé par leurs quatre bornes est également 
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réciproque. Ceci résulte de ce que chaque transduc- 
teur associé au méme point P du milieu constitue, 
comme on l’a vu, un quadripdle électro-acoustique 
réciproque et que le quadripdle électrique T, Ts 
peut étre considéré comme le résultat de leur mise 
en cascade. Or deux quadripéles réciproques mis en 
cascade donnent toujours un quadripéle réciproque. 

L’inverse n’est malheureusement pas exact, car il 
pourrait se faire que le défaut de réciprocité de deux 
quadripéles se compensent exactement dans la mise 
en cascade. Ainsi, bien que la réciprocité du quadri- 
pole électrique T, T, se vérifie trés bien expérimen- 
talement, a la précision des mesures, elle ne consti- 
tue pas une condition suffisante de réciprocité. On 
peut se demander cependant si la réciprocité élec- 
trique de tous les couples T, T,, T, T,, T, Ts; etc. 
de plusieurs transducteurs de types différents ne 
constitue pas une preuve suffisante de réciprocité 
acoustique pour chacun d’eux: 


6. Champ libre entre deux transducteurs 


Sur le transducteur T,, choisissons un point O, 
et une 4 droite O, 6, qui représente la direction 
dans laquelle on veut |’étalonner; soient M, et S, 
sa sensibilité et sa réponse dans cette direction. Sur 
le transducteur T, nous aurons de méme le point O, 
et la $ droite O, 6, A laquelle correspondent M, 
et Sy. 

Placons T, et T, dans l’eau de facgon que les deux 
directions s’opposent sur la méme droite (Fig. 3). 
Les conditions de champ libre seront remplies entre 
T, et T, pour le couple de directions O, 6, , O» 5 
si chacun d’eux se comporte par rapport a |’autre 


aE 


s 
8dB 
20 2,24 251m 


1,26 141 158 178 
Distances ——= 


p Ww 


Fig. 3. Contréle de la condition de champ libre. 
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comme un point source qu’on peut confondre sans 
erreur appréciable avec son point de référence. 

Le critérium pratique est que T,; par exemple 
émettant a courant constant, la tension recue par 
Ty soit inversement proportionnelle a la distance 
O, O,. La réciprocité du quadripéle électrique T, Ts 
donne alors la méme loi, pour une tension recue 
sur T, avec T, comme émetteur. 

La vérification de la condition de champ libre 
qui est trés importante est exécutée avant tout éta- 
lonnage sur des graphiques logarithmiques comme 
ceux de la Fig. 3. On porte en abscisse la distance 
O, O, en dB par rapport a 1 métre. En ordonnée 
on porte la tension recue aux bornes de T,, ou une 
quantité proportionnelle si on utilise un amplifica- 
teur, exprimée en dB par rapport a une valeur quel- 
conque. 

Sur un tel graphique le champ libre idéal se tra- 
duit par une droite de pente négative inclinée a 45°. 
La courbe réelle, dans le cas moyen, présente l’aspect 
de la Fig. 3; elle oscille autour d’une droite idéale; 
les oscillations sont en général assez fortes aux 
faibles distances a cause des ondes stationnaires qui 
s’établissent entre les deux bases, et aux fortes dis- 
tances parce qu’alors l’onde directe est perturbée 
par des réflexions parasites sur la surface, sur le 
fond, ou sur des obstacles qu’on n’a pas pu éliminer. 
I] existe en général une région moyenne ou le champ 
est a peu prés régulier, et qui est la plus favorable 
a la mesure. 

On n’a pas toujours intérét a éloigner le plus 
possible les deux transducteurs lun de l’autre. Une 
telle courbe dépend en fait des directivités combi- 
nées des deux transducteurs, done de la fréquence, 
des dimensions du bassin de mesure et du pouvoir 
absorbant de ses paroies. En particulier, l’emplace- 
ment de mesure limite trés nettement la fréquence 
inférieure a laquelle on peut faire des mesures entre 
qui est représentée, en decibels par l’écart entre la 
courbe deviennent si fortes qu’elle est pratiquement 
inutilisable. 

Quand elle l’est, en tout cas, elle tient compte de 
toutes les causes d’erreur possible sur les conditions 
de champ libre et elle permet de mesurer cette erreur 
qui est représentée, en decibels par |’écart entre la 
courbe et la droite moyenne idéale. 

Ainsi, si on adopte une distance r (Fig. 3) ot 


Pécart vaut ab — il est commode pour la stabilité 
des lectures de se placer au sommet d’une oscilla- 
tion par exemple — on corrigera de la valeur ab, 


la lecture de la tension recue. 


7. Impédance de transfert 


Entre les quatre quantités M,, S;; M,, Ss qu’il 
s’agit de mesurer, la réciprocité nous donne deux 


108 


relations tirées de la formule (6) 


M,/S,=M,/S.=J. (6a) 


La mesure de l’impédance de transfert du quadri- , 


pole électrique T, T, nous en donne une troisiéme, 
pourvu que les conditions de champ libre soient cor- 
rectement réalisées. Soit alors 7, le courant alimen- 
tant T, et e, la tension recue a circuit ouvert sur Ts. 
Cette tension ey a été au besoin corrigée comme on 
l’a vu précédemment de fagon a correspondre au 
champ libre. Dans ces conditions la pression créée 
par T, a l’emplacement de T, en l’absence de celui- 
ci est, comme on I’a vu: 

m= S,i,/r 2 


et la tension e, est égale a M, a d’ou (Fig. 4) 
e, = M, S, i,/r. (9) 


Fig. 4. Mesure de l’impédance de transfert. 


L’impédance de transfert du quadripdéle électrique 
est 


€o/t; = M, S,/r. (10) 


€y et i, se mesurent séparément. Comme on connait 
r on peut calculer la quantité re /i,; que nous ap- 
pelons Z et qui est l’impédance de transfert ramenée 
a la distance unité, dans la pratique, a 1 métre. On 
a donc: 


TOM. Sic (11) 


La réciprocité du quadripéle électrique entraine 
évidemment 


mais la derniére égalité est déja contenue dans la 
formule (6a). Comme on |’a dit, il est facile de véri- 
fier ’égalité des impédances de transfert dans les 
deux sens. 

En définitive, le groupe de relations (6a) et (11a) 
ou J et Z sont des quantités connues, constituent 
trois relations indépendantes. Une quatriéme est 
nécessaire pour calculer M, S, My S,. 


8. Rapport des sensibilités 


La relation supplémentaire nécessaire peut s’ob- 
tenir en mesurant soit le rapport S,/Sy soit le rap- 
port M,/M,. Pratiquement, c’est M,/M, que l’on 
mesure. I] faut disposer pour cela d’un émetteur E 
rapporté lui aussi a un point de référence O et une 
demi-droite O 0. Il est alimenté par un courant qui 
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reste constant pendant toutes les mesures et donne 
ainsi un champ acoustique constant. Placons dans 
ce champ le transducteur T, de facon que les deux 
directions O,6, et OO s’opposent sur la méme 
droite. Vérifions les conditions de champ libre entre 
T, et E. Avec ou sans correction nous mesurons une 
tension e, a circuit ouvert aux bornes de T,, pour 
une distance OO,, égale a 7,. Faisons ensuite la 
mesure homologue avec T,. Nous avons la tension 
ey pour la distance ry. Il est clair que le rapport 
des sensibilités de T, et T, vaut alors 

M,/My= ey 1,/es To. (12) 
On peut donc le calculer et appelons R sa valeur: 


R=M,/M,. 


9. Formules de résolution 


Nous disposons a présent d’assez de relations 
pour létalonnage simultané de T, et T,. Les for- 
mules sont: 


M,/S,=M,/S,=J, J=2/ef; 

MAS Moya M,/M,=R, 
M,=VZJ/R, S,=VZRI/J, 
M,=VZJ/R, So=VZ/I dee 


Remarquons que l’erreur relative sur les M et 
sur les S est égale a la demi-somme des erreurs 
relatives sur Z et sur R. 


La mesure de Z comporte: 
1 vérification de champ libre, 
2 mesures électriques — tension et courant. 


La mesure de R comporte: 

2 vérifications de champ libre, 

2 mesures électriques: deux tensions. 

Un étalonnage par comparaison avec un étalon 
maitre comporterait les mémes opérations que la 
mesure de Z et l’erreur relative serait par consé- 
quent égale, dans des conditions identiques, a l’erreur 
relative sur Z. 

Si on fait le bilan des deux méthodes on s’aper- 
goit que la mesure par comparaison ne l’emporte en 
précision, que si l’erreur relative de |’étalon maitre 
est inférieure a la moitié de l’erreur qu’on peut com- 
mettre sur une vérification de champ. Or, dans la 
gamme de fréquence ot l’emplacement de mesure 
nous a permis d’opérer jusqu’a présent, il est assez 
facile — les graphiques de vérification de champ 
le montrent — d’obtenir une précision de +0,3 dB 
(+3,5%). Il faudrait donc un étalon maitre stable 
et précis 4 mieux que 3%. Nous ne pensons pas 
qu’un tel étalon, s’il existe, puisse étre utilisé dans 
les conditions proches de la pratique oi nous opé- 
rons. C’est pourquoi |’étalonnage direct des trans- 
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ducteurs d’usage courant nous parait préférable. 
Comme on en étalonne deux a la fois, le bilan des 
temps de travail n’est pas lui non plus défavorable ?. 

Enfin, la méthode posséde un grand avantage 
psychologique, parce qu’elle dispense précisément 
d’avoir confiance dans un étalon maitre, dont on 
peut toujours redouter les variations ou les détério- 
rations, 


10. Précision calculée 


10.1. Un étalonnage comporte trois vérifications de 
champ et quatre mesures électriques. La pré- 
cision des verifications de champ dépend des dimen- 
sions du bassin de mesure et en particulier de sa pro- 
fondeur. Actuellement, pour les fréquences supérieu- 
res a 8kHz, nous comptons sur une précision de 
+0,3 dB en général. Nous pouvons encore travailler 
entre 6 et 8kHz a +0,5 dB. En dessous de 6 kHz 


les mesures sont trés approchées. 


10.2. En ce qui concerne les mesures électriques, 

nous n’avons pas encore cherché a pousser la 
précision au maximum. Avec un ampéremétre a 
thermocouple et des voltmétres fréquemment con- 
trélés, nous comptons également sur une précision 
de £0,3 dB (3,5%). 


La demi-somme des sept erreurs donne donc: 


+ 1,05 dB au-dessus de 8 kHz, 
+ 1,35 dB entre 6 et 8 kHz. 


Comme il y a bien peu de chance que sept erreurs 
s’ajoutent intégralement, nous ne pensons pas étre 
optimistes en comptant sur une précision de + 1 dB 
au dessus de 8 kHz et +1,2 dB entre 6 et 8 kHz. 
Nous verrons qu’une mesure plusieurs fois répétée 
donne une convergence meilleure. 

La précision de +1 dB qui était notre but initial 
est a peu prés suffisante pour nos besoins. Pour 
Yobtenir a coup sir en dessous de 8 kHz un bassin 
plus vaste, plus profond, suffirait. Mais l’intérét de 
ces fréquences n’a pas paru justifier jusqu’a présent 
la dépense a engager. Pour des fréquences beaucoup 
plus basses, par exemple inférieures a 3 kHz, la mé- 
thode en champ libre cesse d’étre praticable et il 
faut faire appel a d’autres méthodes. 

Les vérifications de champ qui sont a la base de 
la méthode sont, dans la pratique courante, trés 
simplifiées par l’entrainement des opérateurs. Comme 
on a trés souvent le méme type de transducteur a 
étalonner, on sait assez vite que pour un couple 
de type donné, a une fréquence donnée, la précision 
de champ est suffisante et i] n’ést pas nécessaire de 
refaire chaque fois la mesure. 


_ ® Une mesure supplémentaire de tension recue, per- 
met, au besoin, d’étalonner du méme coup le transduc- 
teur qui sert de source auxiliaire FE, dans la mesure 
de R. Le bilan deyient alors nettement favorable. 
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11. Installation 


La Marine Nationale a installé a proximité de son 
Laboratoire du Brusc un chaland-laboratoire ou se 
font les mesures d’étalonnage. Ce chaland est installé 
sur une fosse draguée dans un haut fond (Fig. 5 47), 
a 200 m environ du rivage. La fosse a une vingtaine 
de métres de dimension horizontale et entre 4 et 
5 metres de profondeur. Le fond du chaland est 
percé de deux cuves. Celles-ci sont bordées de rails 
qui servent a des chariots. Ceux-ci sont a leur tour 


Figs. 5, 6 et 7. Vues du chaland-laboratoire du Bruse. 
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parcourus par des chariots plus petits qui supportent 
les axes de suspension des transducteurs. 
On peut ainsi déplacer ceux-ci dans deux dimen- 


sions: les faire tourner autour d’un axe vertical et | 


régler leur immersion. On travaille en général a des 
immersions comprises entre 2m et 2,50 m. 
Le banc de mesure comporte (Fig. 8): 


a) du cété émission 

— un générateur de tension sinusoidale, 

— un amplificateur de puissance 10 watts dont la 
sortie est adaptable aux impédances des trans- 
ducteurs mesurés, 

— un ampéremétre a thermocouple (3%) est mis en 
circuit pour faire la mesure du courant puis aussi- 
tot remplacé par une résistance égale, 

— un voltmétre V, mesure la tension aux bornes 
de |’émetteur, 

— éventuellement, on peut controler le taux d’har- 
moniques de la tension sur |’émetteur; 


Atténuateur i he 
Amplificateur 
a filtre étroit 


Indicateur 


Générateur 


Arplificateur 
de puissance 


Fig. 8. Schéma de principe du banc de mesures d’éta- 
lonnage. 


b) du cété réception 

— un amplificateur 4 haute impédance d’entrée (su- 
périeure 4 100 kQ entre 5 et 50 kHz) qui permet 
de recevoir a circuit ouvert entre 5 et 50 kHz sur 
tous les hydrophones utilisés jusqu’ici. 

— Cet amplificateur est accordé sur la fréquence de 
mesure entre 5 et 50 kHz par un changement de 
fréquence et un filtre de 500 Hz de largeur de 
bande. Ce filtre étroit élimine les bruits parasites 
et permet de mesurer des niveaux sonores de 
— 950 dB par rapport a une barye, 

— un indicateur en bout de chaine permet de repérer 
les tensions a la sortie de l’amplificateur. Cet in- 
dicateur a son cadran amorti par de fortes con- 
stantes de temps de facon a pouvoir faire des 
lectures correctes m€me par houle moyenne, 

— une commutation permet de couper la chaine 
acoustique et d’injecter a l’entrée de l’amplifica- 
tion, une tension réglable par un atténuateur et 
mesurable avec le voltmétre V,. 
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La mesure de la tension recue a circuit ouvert sur 
un hydrophone se fait done en repérant, aprés accord 
de l’amplificateur, la déviation de l’indicateur, En- 
suite, on coupe la chaine acoustique et on injecte sur 
lampli une tension de référence mesurée avec pré- 
cision par V,. Ce procédé permet d’éliminer les 
erreurs dues a |’amplificateur. 

L’impédance de transfert est donc, en définitive, 
le rapport d’une tension de référence et du courant 
dans |’émetteur. 

Le rapport des sensibilités est le rapport de deux 
tensions de référence. Toutefois, pour des valeurs de 
ce rapport ne dépassant pas 5 ou 6 dB, on utilise 
directement le cadran de l’indicateur gradué en 
decibels et qui donne une précision suffisante. 


12. Réponse par volt 


Nous avons appelé S la réponse en émission, c’est- 
a-dire la pression sonore 4 1 métre par unité de 
courant ou encore les baryes par ampere a 1m. On 
peut également définir la réponse par volt a 1 m qui 
est égal au quotient de la précédente par le module 
de Vimpédance. Cette réponse est plus facile a uti- 
liser pour qui veut produire avec le transducteur un 
niveau connu, car on dispose plus fréquemment d’un 
voltmétre que d’un ampéremétre. C’est pourquoi, 
on a placé dans le banc de mesure le voltmétre V, 
qui donne la tension aux bornes de |’émetteur. On 
peut ainsi calculer le module de l’impédance et par 
conséquent la réponse par volt. 

Pour le transducteur placé en réception, une me- 
sure d’impédance supplémentaire est nécessaire. Les 
mesures des réponses par volt exigent donc une ou 
deux mesures électriques de plus et par conséquent 
comportent des erreurs supplémentaires. 


13. Convergence des résultats 


La courbe de la Fig.9 donne la valeur de la 
sensibilité d’un transducteur piézoélectrique entre 


oo 
is) 


1 
foe} 
[Sal 


Sensibilité @ la réception en dB par rapport 41/1 barye—>- 


pe Re 0. > Meo Tae 


Fréquence —> 


Fig. 9. Etalonnage par réciprocité d’une base piezo- 
électrique (sensibilité). 


— ‘ 
: 


ACUSTICA 
Vol. 8 (1958) 


6 et 30 kHz. Chaque point de mesure correspond 
a une mesure par réciprocité faite avec un trans- 
ducteur auxiliaire différent, une source E différente 
et des conditions de mesure variables. 

Si l’on excepte les deux fréquences extremes ot 
la pente de la courbe justifie une certaine dispersion, 
on voit qu’on a obtenu couramment sur la méme 


fréquence, plus de 5 points de mesure dont la dis- 


_ persion totale ne dépasse pas 1 dB. Ceci autorise a 


. 
| 


penser que la précision est bien proche de + 0,5 dB. 

La Fig. 10 concerne la réponse par volt du méme 
transducteur. La dispersion est un peu plus forte 
mais laisse espérer une précision de +0,7 dB. 

Sans nous engager sur ces chiffres qui ont besoin 
dune confirmation statistique plus complete, il 
semble que la précision «officielle» de +1 dB soit 
asez facilement tenue. 


“i 
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Réponse par volt en dB par rapport 1 barye/1V a 1m —= 
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wm 
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Fig. 10. Etalonnage par réciprocité d’une base piezo- 
électrique (réponse par volt). 


14. Utilisation éventuelle des impulsions, 
inconvénients 


On peut réaliser la condition du champ libre 
méme dans une cuve de petites dimensions, a con- 
dition d’opérer en impulsion et de ne tenir compte 
sur le récepteur que de l’onde directe séparée des 
échos parasites par le balayage d’un oscilloscope. 
Une telle méthode permet de se passer d’un emplace- 
ment de mesure de grandes dimensions. Mais |’ob- 
jection majeure a un tel procédé est que les mesures 
électriques deviennent des mesures a l’oscilloscope 
et perdent ainsi beaucoup de précision par rapport 
aux mesures en régime permanent. 

De plus, V’obligation d’utiliser des impulsions 
courtes (de l’ordre de la milli-seconde) rend im- 
possible ces mesures sur la fréquence de résonance 
d'un transducteur dont le coefficient de surtension 
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mécanique serait élevé. En effet, dans ces conditions 
ni le courant dans |’émetteur, ni la tension sur le 
récepteur n’ont le temps de s’établir et la définition 
de leur valeur devient difficile. 

Enfin, il n’est pas exclu a priori que, sur des im- 
pulsions trés courtes, les réponses des transducteurs 
ne soient pas exactement les mémes qu’en régime 
permanent en raison de phénomeénes transitoires 
éventuels. Les mémes objections peuvent étre adres- 
sées aux méthodes de self-réciprocité. C’est pour- 
quoi il a paru préférable de réserver les méthodes 
d’impulsions pour des étalonnages moins précis faits 
par comparaison, pour rechercher le maximum de 
précision en opérant en régime permanent et, dans 
la mesure du possible, loin de tout obstacle. 


15. Possibilités d’amélioration 


Rien n’exclut, avec la méme méthode, une pré- 
cision plus grande. On pourrait gagner sur la con- 
dition de régularité de champ, dans un volume d’eau 


' plus vaste et plus calme que la fosse marine utilisée. 


Ceci est vrai surtout aux basses fréquences, mais 
suppose l’aménagement d’un site — petit lac ou 
plan d’eau calme et assez profond éloigné des sour- 
ces de bruit — donc d’assez gros moyens matériels. 
Il faudrait également rapprocher les procédés de 
mesure de ceux des Laboratoires de Métrologie. On 
atteindrait certainement assez aisément la précision 
de +0,5 dB a chaque mesure. 

Aux fréquences élevées — de l’ordre de 100 kHz 
— ou, pour des directivités équivalentes, les émet- 
teurs sont plus petits et plus maniables, ot les ondes 
sonores sont plus facilement absorbées, les condi- 
tions de champ libre sont plus faciles a réaliser dans 
des cuves de dimensions restreintes. La méthode 
peut s’appliquer dans les conditions mémes de la 
recherche de Laboratoire. Elle parait étre ainsi mieux 
adaptée encore a ces fréquences qu’aux fréquences 
usuelles de l’acoustique sous-marine. 

(Regu le 11. Novembre 1956.) 
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SIGNAL/NOISE PERFORMANCE OF SUPER-DIRECTIVE ARRAYS 


by D. G. Tuckrr 
Electrical Engineering Department, University of Birmingham 


Summary 


A super-directive array is defined as one whose effective aperture exceeds its physical 
aperture, and it is shown that while such an array has a Directivity Index somewhat better 
than that of an ordinary array, nevertheless its Noise Factor is much worse. Consequently, 
in many receiving systems, the use of super-directivity is not be recommended. The con- 
ception of super-directivity is extended to bearing-determining arrays which give a null 
response on the axis, and it is shown that the same disadvantage of a bad Noise Factor 
applies to them also. CES 


Sommaire 

On dit qu’une baie de réception électro-acoustique est «superdirective» lorsque son 
ouverture efficace est supérieure 4 son ouverture géométrique. Un tel dispositif a un indice 
de directivité quelque peu supérieur a celui d’une baie normale, mais on montre que son 
facteur de bruit est beaucoup moins bon: de sorte que de recourir a la superdirectivité 
n’est pas 4 recommander, dans beaucoup de problemes de réception. 

On étend la notion de superdirectivité aux baies de goniométrie a réponse nulle sur 
Vaxe: on montre qu’elles aussi présentent l’inconvénient d’un mauyais facteur de bruit. 


Zusammenfassung 


Als Anordnung mit Uber-Richtwirkung wird eine solche bezeichnet, deren effektive Off- 
nung groBer ist als die physikalische. Es wird gezeigt, dal} eine solche Anordnung einen 
etwas besseren Richtfaktor hat als eine gewoéhnliche, jedoch ist die Rauschzahl ungiinstiger. 
Folglich ist fiir viele Empfangssysteme die Anwendung der Uber-Richtwirkung nicht zu 
empfehlen. Das Prinzip der Uber-Richtwirkung wird ausgedehnt auf Peilanordnungen mit 
Nullanzeige in Achsenrichtung. Auch fiir diese Systeme bestehen die Nachteile der un- 


giinstigen Rauschzahl. 


1. Introduction 


This paper is a sequel to a previous one [3] 
which developed a new approach to the calculation 
of signal/noise performance of arrays. It was shown 
that, while the conventional directivity factor (or 
directivity index) measures the signal/noise per- 
formance of a beam-forming array as far as noise 
received from the medium is concerned, it does not 
provide a true measure of the performance in respect 
of thermal agitation noise in the transducers or 
receiving amplifiers, except when the array is large 
so that all the effective angular range of response 
is confined to a comparatively small angle. For tak- 
ing account of thermal agitation noise, and to allow 
for arrays which do not form a conventional beam 
but instead produce a null on the axis, and for ar- 
rays in which the outputs from two halves are multi- 
plied together, the signal/noise performance factor 
(called Noise Factor) was proposed. This Noise 
Factor in the ratio by which the signal/noise per- 
formance of a receiving system is worse than that 
of the same array used with uniform and equi- 
phase sensitivity along its length. In addition to a 
rather wider range of application, Noise Factor has 
the advantage over directivity factor of being much 
more easily calculated. 


The present paper is concerned with arrays of the 
type frequently called “Super-Directive”. While no 
formal definition of this term appears in standard 
glossaries, it can probably be defined by usage as 
describing an array whose effective aperture exceeds 
its physical aperture; in other words the super-direc- 
tive array has a directivity index (or power gain in 
the terminology of electro-magnetic aerials) greater 
than that of the uniform and equi-phase array‘. 
This result can be obtained only at the expense of 
a very substantial response in the range of imagi- 
nary angles, where |sin@|>1, 0 being the angle of 
any direction relative to the normal to the face of 
the array. 

It might be thought that the super-directive array, 
having such a favourable directivity factor, would 
have considerable merits as a transmitting array. 
This is in principle true enough, although in prac- 


1 This use of the term “Super-directive” appears uni- 
versal in published work on the subject; see, for ex- 
ample, Prircuarp [5], [6]. It should be noted, how- 
ever, that Prircuarp gives the impression that “super- 
directive arrays” and “arrays with reversed-phase ex- 
citation” are synonymous; but, in fact, reversed-phase 
excitation (or sensitivity, in the case of a receiving ar- 
ray) can be used advantageously for purposes other 
than to obtain super-directivity —see, e. g., Section 10 
of reference [1]. 
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tice it may be difficult to realise with sufficient ac- 
curacy the required amplitude distribution over the 
array; and the power efficiency is low. For re- 
ception, however, the use of superdirectivity is, 
even in principle, of questionable value, since it 
will be shown in this paper that the Noise Factor 


' is worsened at a much faster rate than the direc- 


tivity factor is improved. It is in fact in relation 
to super-directive arrays that the use of the Noise 


‘Factor as a criterion of performance is most 


valuable. In the previous paper where the Noise 
Factor was discussed in relation to ordinary arrays, 
its merits were less emphatically established. 

For simplicity, discussion in this paper, as in the 
previous papers, is confined to strip or line arrays, 
i.e. to arrays of negligible width, although of finite 


- length. This means that the directional patterns can 
_ be considered in two dimensions only. This simpli- 


fication does not affect any of the principles in- 
volved. It is also assumed in the illustrative examples 
that the arrays are formed from a continuous trans- 
ducer divided into equal sections. Only “broadside” 
arrays are specifically considered —i. e. only arrays 
whose acoustic axis is normal to the line of the 
array — but “endfire” arrays give the same kind of 


_ result. 


Calculations of the distribution of sensitivity over 
the length of the array, and of the number of sec- 
tions required, are not given here, as their basis has 
been given in an earlier paper [1]. 


2. The design of a super-directive array 


As explained above, the super-directive array dif- 


fers from the ordinary one in having a substantial 


response in the range of imaginary angles, i. e. where 
|sin@|>1. Directional responses are normally cal- 
culated in terms of a variable #, which is given by 


x sin 0, 
where / is the length of the array and / is the wave- 
length. On the x-scale, therefore, the limit of real- 
angle response occurs where x= +/1x//. It is clear 
then that the design of a super-directive array de- 
pends on knowledge of the length (or size in the 
general case) of the array and cannot be worked out 
in more general terms. Reference has already been 
made to the region of imaginary angles, where 
|x|>In/2; the significance of this range will be 
discussed in more detail later. For the time being 
it need merely be said that this range plays an im- 
portant part in determining the signal/noise perfor- 
mance of the array. 

The design problem of any array is of course to 
find the distribution of sensitivity and phase, over 
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the length of the array, required to give the desired 
directional pattern. We are, in this paper, concerned 
only with receiving arrays, but the same directional 
response is obtained on transmitting if the distribu- 
tion of excitation corresponds to the distribution of 
sensitivity discussed here. The directional pattern 
and the distribution of sensitivity are related by the 
Fourier transforms shown below: 


+ occ 


D(x) = f eT (r) dr, 


+ co 
DG) =a cewe = Dard 


— co 


where T(r) is the so-called tapering function, i. e. 
the distribution of sensitivity over the length of the 
array, and D(x) is the directional pattern on the 
scale of x. Obviously D(x) is different at different 
frequencies, since x is a function of J and therefore 
these equations apply at any one particular fre- 
quency. It would in principle be possible to de- 
termine a suitable function, 7'(r), from these equa- 
tions to give the required pattern D(x). The pro- 
cess would in practice be difficult. A much simpler 
technique was described in [1], where the pattern 
of any desired degree of complexity can be synthe- 
sised by the super-position at suitable amplitudes of 
a number of elementary patterns. These elementary 
patterns are based on the pattern D(x) = (sin 2) /2 
which is obtained from the array when the sensiti- 


vity is uniform and the phase is constant over its 
length. By subjecting the successive sections of the 
array to successive phase-shifts (e. g. by connecting 
the sections to a delay line) the pattern can be de- 
flected to one side or other of the axis, so that it 
takes the form 


sin (x tn It) 
atnn 


The process of synthesis is thus very simple and its 
details need not be repeated here. When symmetrical 
patterns are being produced, it is unnecessary to use 
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Fig. 1. Example of synthesis of super-directive response 
with reduced secondary lobes and slightly-nar- 
rowed main beam. 
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real delay lines, since the overall effect is easily 
shown to be realised by giving different sensitivities 
to different sections of the array. 


The design of super-directive arrays by this pro-’ 


cess of synthesis is quite straightforward. It simply 
involves super-posing on any ordinary directive pat- 
tern some other patterns which have their peak 
responses in the range of imaginary angles”. Their 
secondary lobes are used to modify the ordinary 
response and in particular to cancel out its secon- 
dary lobes. The process will be made clear by some 
examples, 

Figs. 1, 2 and 3 show three simple cases of the 
synthesis of a super-directive response. In Fig. 1, 
the full line curve shows a normal (sin x) /x pattern 
of unity peak height. It will be assumed that the 
length of the array is twice the wavelength, so that 
the limit of real angle response is given by x= 2a. 
Within this range fall the first secondary lobes of 
the pattern with an amplitude of 0.22. These se- 
condary lobes can be reduced to very small values, 
with an accompanying narrowing of the main lobe, 
by adding in opposite sense the two deflected pat- 
terns shown by dotted lines. Each of these has a peak 
height of —1.0 occurring at r= +32, i.e. in the 
range of imaginary angles. It is seen that the resul- 
tant pattern shown in dashed lines is, in the range 
of real angles, very considerably improved. Similarly 
in Fig. 2, two patterns with peak amplitude of —1.5 


have been added to the ordinary pattern with still, 


further improvement in directionality. 


Fig. 3 shows an example of a slightly different 
kind. Here it is assumed that the array has a length 
of A, so that the limit of real angle response is given 
by «= +a. The aim has been to obtain a beam of 
one-third the width, in the 6-scale, of the ordinary 


2 The definition of super-directivity in terms of large 
response in the region of imaginary angles seems to 
have many advantages in spite of the difficulty of pro- 
viding a simple physical explanation of the significance 
of imaginary angles. For instance, it follows directly 
from this definition and from the delay-line method of 
synthesis that, since major peaks of response are to fall 
outside the range of physical angles, i. e. outside +x/2 
radians from the normal to the array, then the spacing 
between sections of the array cannot exceed half-a- 
wavelength. This is because the electrical delay net- 
works which couple the sections together cannot give 
phase differences between adjacent sections of more 
than x radians for steady-state signals in which one 
cycle is indistinguishable from another. It also follows 
from the definition and from the work of [3] that the 
thermal-agitation noise of the dissipation resistance of 
the array is proportional to the power gathered (notio- 
nally) by the whole directivity pattern including imagi- 
nary angles, whereas noise arising in the medium is 
received by the portion of the directivity pattern cor- 
responding to real angles only. This accounts for the 
poor Noise Factor discussed in Section 3. 
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beam corresponding to this length of array, or on 
the x-scale with a width of one-half that of the or- 
dinary pattern. This has been accomplished by ad- 
ding two deflected patterns with their peaks at ~ 
x= +22, and amplitudes calculated to give a zero © 
resultant at «= 2/2. This necessitated the ampli- 
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Fig.2. As in Fig. 1, but with secondary lobes prac- 
tically eliminated. ; 
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Fig. 3. Example of synthesis of super-directive response 
with main beam greatly narrowed. 


tude of the peaks in the range of imaginary angles 
being 7.5 times that of the ordinary pattern*. The 
resultant pattern shown in the lower graph in Fig. 3 
is seen to have secondary lobes, in the range of real — 
angles, of peak amplitude 0.44, i.e. rather larger 
than those of the ordinary pattern, but the main 
beam is narrowed, as desired. The directivity index 
of the new pattern is approximately 3 dB better than 
that of the ordinary (sin z)/x pattern. Evidently 
the secondary lobes can be reduced by the addi- 
tion of further deflected patterns. Even greater 
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3 It may be of interest to note that the variation of 
relative sensitivity along the length of the array is in 
this case given by 1+15 cos2ra, where the array ex- 
tends in length from r= —1 to r=+1. The minimum 
number of sections required in the array is five, and 
the maximum spacing between centres of the sections 
is therefore 2/5. 
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improvements can be obtained by the super-posi- 
tion of patterns whose zeros do not coincide. This 
_ matter is described by Woopwarp and Lawson [4], 
who show moreover the price that has to be paid 
_ for this improved directivity in terms of difficulty 
of construction. In extreme cases, extreme accuracy 
of distribution of sensitivity over the length of the 
array is required. In one of their examples a varia- 
_ tion of amplitude of the thickness of the line in the 
graph is sufficient to remove the beam altogether. 


_ It does not seem to have been pointed out pre- 
_ viously that in addition to the practical price paid 
i for super-directivity there is a fundamental price 

paid in the worsening of the Noise Factor *. This 
3 will be described in the next section, 


____ 3. Noise factor of super-directive arrays 


__ Even though the use of delay lines to produce the 
deflected beams used in the process of synthesis re- 
ferred to above is not necessary in practice, never- 
4 theless the retention of these delay lines as a stage 
; in the design process is very helpful in the calcula- 
e. tion of Noise Factor. Clearly, when a number of 
_ patterns are superposed the noise powers appearing 
at the terminals of the delay lines are added. It was 
f shown in [2] that for the values of phase-shift used 
_ in the delay lines for this process the outputs of 
_ noise from the various delay lines are all uncorre- 
_ lated with one another provided the noise voltages 
{ in the various sections of the array are themselves 
uncorrelated °, Thus it is clear that the improvement 
of directivity factor by the superposition of patterns 
with peaks in the range of imaginary angles adds 
merely to the noise output, while contributing no- 
thing to the peak signal output. The Noise Factor 
must therefore be considerably worsened. Moreover, 
while the improvements which have been achieved 
in the directivity factor are comparatively small, the 
worsening of the Noise Factor is comparatively 
high. For example, in Figs. 1, 2 and 3, the Noise 
Factors are worsened by approximately 5, 7 and 
20 dB respectively. In Fig. 3, then, an improvement 


4 Tt should be noted that Prircuarp [5] has given 
a hint of this in his statement, “However, it is possible 
that applications (of super-directive arrays) may exist 
‘in the field of acoustics, especially where the limiting 
noise in a receiving system is ambient noise in the 
‘medium rather than the limiting thermal noise in the 
transducers”’. 

5 These noise voltages are always uncorrelated as far 
as thermal-agitation noise in the dissipation resistance 
‘of the transducers is concerned; but this may not be 
so for noise arising in the medium, since if the sections 

f the array are very closely spaced, as in super-direc- 
tive arrays, there is bound to be some correlation, even 
‘if the noise is uniformly distributed in space. 
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of 3 dB in directivity is paid for by a worsening of 
20 dB in Noise Factor. For a transmitting array 
where noise factor is not relevant, and when ample 
power is available, no doubt super-directivity is 
worth attempting. On reception it will usually be 
undesirable. 


The extent to which super-directivity may be an 
advantage or a disadvantage for receiving arrays 
depends on the source of background noise against 
which the signal is to be detected. If the noise arises 
in the medium—for example, in an under-water 
array such noise would comprise thermal-agitation 
noise in the water and also what is generally called 
“sea-state’’ noise, being generated by the breaking 
of waves, by the rolling of shingle, etc., — then since 
the generation of this noise is physically confined 
to the range of real angles, clearly the directivity 
factor or index would be expected to be an appro- 
priate criterion of performance, and super-directivity 
would be expected to be an advantage. If on the 
other hand the dominant noise is that arising from 
thermal agitation in the dissipation resistance of the 
transducers forming the array together with the 
noise arising in the receiving amplifiiers (the signi- 
ficance of the latter was discussed in [3]), then the 
detection of signals against this background depends 
on the Noise Factor, and clearly super-directivity 
is a disadvantage. It should be remembered that a 
system where the basic limitation is noise from the 
medium can easily be changed to one where the limi- 
tation is resistance and amplifier noise by the use 
of super-directivity. 


4. Extension of the principle to 
bearing-determining systems 


In both references [1] and [3], some discussion 
was given of the design and performance of bearing- 
determining arrays whose directional pattern had 
a null at x=0, i.e. for signals arriving in a direc- 
tion normal to the line of the array. The principle 
of super-directivity, interpreted as referring to direc- 
tional patterns which have major peaks in the range 
of imaginary angles (i.e. in the range where 
|sin®|>1), can also be applied to this kind of 
array, and it will be shown by means of an example 
that the same disadvantages arise. 

In the previous papers, it was proposed that the 
signal/noise performance of bearing-determining ar- 
rays be measured as the ratio of the slope of the 
directional pattern through x=0 to the background 
noise. Assuming that multiplying systems are ex- 
cluded from consideration (these were fully dis- 
cussed in [3]), a good array arangement for bear- 
ing determination was shown to be that with “sine 
taper’, having the directional pattern 
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The slope of this pattern at x = 0 is easily seen to be 
8/x?. This gives a Noise Factor of —1.2 dB (the 
reference array being the standard “split-beam” 
array). 

A super-directive array, say of length 1, can be 
arranged to have the same slope at x=0. An ex- 
ample of a suitable basis of design is shown in 
Fig. 4. Two deflected patterns with peaks at +2 
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Fig.4. The principle of super-directivity applied to a 
bearing-determining array. 


and —2 2 on the z-scale, and with peak amplitudes 
of —8/nx and +8/x respectively, are added, and 
the resultant pattern 


8 |sin(e+2m)  sin(w—27n) 
Te x+27 x—-20 


Die 


is easily shown to have the same slope through x =0 
as the pattern D,. The response in the range of real 
angles is made up entirely of the secondary lobes 
of the component patterns. 

It can be seen by inspection of the resultant pat- 
tern, shown in Fig. 4, in which the real-angle res- 
ponse of D, is also shown, that as far as noise from 
real space is concerned, the super-directive pattern 
Dy, is superior to D,, the improvement being rather 
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less than 3 dB. But, when noise from the dissipation 
resistance and input circuits is considered, it is seen 
that the Noise Factor is worse by a ratio of 8/z, i. e. 
by 8.1 dB. Thus the use of super-directivity has the 
same objections in a bearing-determining system 
as in a beam-forming system. 


5. Conclusions 


In receiving arrays of super-directive type the 
directivity index can only be a significant measure 
of performance if the limiting factor is the noise 
arising from distant sources in the medium. When 
performance is dominated by the noise from the dis- 
sipation resistance of the transducers or from the 
receiving amplifiers, then the super-directive array 
is inferior to the ordinary array. The best Noise 
Factor that can be obtained is that of the plain 


array with uniform sensitivity. 
(Received March 21st, 1957.) 


List of Symbols 


Bae ae 

v=, sin 6, 

0 physical angle of any direction relative to the normal 
to the face of the array, 

l length of array, 

2 wavelength. 


References 


[1] Tucker, D.G., Some aspects of the design of strip 
arrays. Acustica 6 [1956], 403. 

[2] Tucker, D.G., The randomization of combinations 

of noise sources by means of delay lines. Acustica 

@ (195712399. 

Tucker, D.G., The signal/noise performance of 

electro-acoustic strip arrays. Acustica 8 [1958], 53. 

Woopwarp, P. M. and Lawson, J. D., The theoretical 

precision with which a radiation pattern may be 

obtained from a source of finite size. Proc. Inst. 

Elect. Engrs. (III) 95 [1948], 361. 

Prircnarp, R. L., Optimum directivity patterns for 

linear point arrays. J. acoust. Soc. Amer. 25 

[1953], 879. 

Pritrcuarp, R. L., Maximum directivity index of a 

linear point arrays. J. acoust. Soc. Amer. 26 

[1954], 1034. 


[3] 
[4] 


[5] 


[6] 


Hsia 


a 


ee 


SOUND ABSORPTION OF A STONE WALL 
by E. J. Evans * and P. H. Parkin ** 


Summary 


A Hollington stone wall at one time proposed for use in the new Coventry Cathedral 
was required to have an absorption coefficient of at least 0.2 between 250 and 1000 c/s. 
The absorption coefficient of the stone itself was less than 0.05 at all frequencies. The 
necessary absorption was obtained by leaving 2.5 cm gaps between the vertical faces of the 
stone blocks, these acting as resonators. They were filled with glass-wool and varied in 
depth to give the required broadband absorption. 


Sommaire 


On a demandé qu’un mur en pierre de Hollington pour la nouvelle cathédrale de 
Coventry ait un coefficient d’absorption d’au moins 0,2 entre 250 et 1000Hz. Or, le 
coefficient d’absorption de la pierre elle-méme était inférieur 4 0,05 pour toutes les fré- 
quences. La valeur demandée de l’absorption a alors été obtenue en laissant des intervalles 
de 2,5 cm entre les faces verticales des blocs de pierre, ces intervalles agissant comme des 
résonateurs. On les a remplis de laine de verre et leur a donné différentes profondeurs 
pour obtenir la largeur youlue de la bande d’absorption. 


Zusammenfassung 

Die seinerzeit fiir die neue Kathedrale in Coventry vorgeschlagenen Wande aus Holling- 
ton-Stein miissen mindestens einen Absorptionskoeffizienten von 0,2 im Bereich 250 bis 
1000 Hz haben. Der Absorptionskoeffizient des Steines selbst ist bei allen Frequenzen 
geringer als 0,05. Die erforderliche Absorption wird durch 2,5 cm breite Spalten zwischen 
den Steinen, welche als Resonatoren wirken, erzielt. Sie sind mit Glaswolle gefiillt und in 


unterschiedlicher Tiefe ausgefiihrt, um die notige Bandbreite zu erzielen. 


At an early stage of the design of Coventry Cathe- 
dral now being built in England, the internal wall 
surfaces were to be of stone. This has now been 
changed, but the sound absorption treatment that 
was devised for the stone may be of general interest. 


The absorption coefficient of the wall had to be 
at least 0.2 over the medium frequency range (about 
250 to 1000c/s) decreasing to 0.1 at lower and 
higher frequencies. The wall consisted of blocks of 
Hollington stone, (a fairly fine-grained, porous, 
triassic sandstone of moderate to low hardness), all 
the blocks being of the same depth (18.5 cm), but 
varying in width at random between 20cm and 
100 cm, and in height between 15cm and 40 cm; 
the stone was to form a facing only, with a solid 
backing. The solution adopted was to leave 2.5 cm 
wide vertical openings between each block. Half of 
these slots went back the full depth (18.5 cm) of the 
blocks; the other half were filled in solid from the 
back to half the depth. All slots were filled nearly 


to the surface with glass-wool. 


The absorption coefficients were measured in the 
reverberation chamber of the National Physical 
Laboratory (volume about 275’m*). The specimen 
wall 3.1 mx 3.1 m was laid on the floor. The “ver- 
* National Physical Laboratory, Teddington, Middle- 


sex. 
** Department of Scientific and Industrial Research, 


Building Research Station, Garston, Watford, Herts. 


tical” joints were left open as described; the hori- 
zontal joints were made with mortar in the usual 
manner. The ends of the “vertical” joints at the 
edges were sealed. The total number of vertical 
joints was 78, i. e. about one per 0.12 square metre. 

The absorption coefficients obtained are given 
in the Table, which also includes the results of two 
preliminary tests which provided the basis for the 
final design. 

The data in the Table refer to the following treat- 
ments; 


1. All slots 18.5 cm deep. Slots empty. 
2. All slots 18.5 em deep. Slots filled with glass- 


wool. 


3. One half slots 18.5 cm deep. Other half 9.3 em 

deep. Slots filled with glass wool. 

The absorption coefficients of the stone wall with 
all the joints closed were less than 0.05 at all fre- 
quencies. 

The results indicate that the slots act as resonator 
absorbers. Maximum absorption occurs when the 
depth is approximately an odd multiple of a quarter 
wave length, and the frequencies correspond to those 
of a stopped organ pipe of length approximately 
equal to the slot depth. For example, the empty slots 
show absorption maxima at 400c/s and 1200 c/s, 
with a minimum at about 800 c/s. These frequencies 
correspond to the resonance and anti-resonance fre- 
quencies of a stopped pipe 21.5cm long, i.e. the 
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Table I 
Reverberation absorption coefficient 
Mean (to nearest 0.05) 
eine Treatment 
(c/s) eee Ep af AA a er ae 
| 1 | 2 | 3 
125 Vet RUG) 0.15 | 0.10 
160 0.10 0.15 0.10 
200 0.15 0.15 0.15 
250 0.10 | 0.20 0.20 
320 | 0.30 0.35 0.25 
400 | 0.50 0.30 0.25 
500 0.20 0.20 0.25 
640 0.10 0.20 0.25 
800 0.05 0.20 | 0.25 
1000 0.05 0.25 0.25 
1250 0.15 0.20 0.25 
1600 0.10 0.20 0.20 
2.000 0.10 0.15 0.20 
4000 0.15 0.15 0.15 


depth of the slot plus an end correction of about 
3 cm. The introduction of the absorbent in treatment 
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2 gives slightly lower frequencies and broader 
maxima. It is clear that the resonance frequencies 
can be controlled in a simple manner by variation 
of the slot depth, and in treatment 3 a combination 
of two depths was adopted to provide uniform ab- 
sorption over a wide frequency range. 

The present measurements show that the slot 
resonator forms a useful addition to the range of 
available acoustical treatments. In particular it ap- 
pears to have application as a simple alternative to 
the Helmoltz resonator absorber. 
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4. Tonmeistertagung in Detmold 


Die vierte dieser schon Tradition gewordenen Tagun- 
gen fand im Einverstandnis mit der ICA vom 22. bis 
25. Oktober 1957 an der Nordwestdeutschen Musik-Aka- 
demie statt (Leitung E. Turennaus). Fast 200 musika- 
lische und technische Fachleute, u. a. aus acht aufer- 
deutschen Landern, waren zusammengekommen zu Re- 
feraten, Vorfiihrungen und gemeinsamen Diskussionen 
iiber kiinstlerische und technische Probleme der Musik- 
aufnahme und -wiedergabe, insbesondere iiber die Fort- 
schritte der stereophonischen Ubertragung. Hervorgeho- 
ben seien folgende Themen: 

Die geistige Welt des Musikinstrumens (L. June), 
Glockenspiele und Glockenspielmusik in den USA. (Prr- 
civat Price), Neuentwicklung in der raéumlichen Schall- 
wiedergabe (N. V. Franssen), Uber die Darstellung 
der Struktur ausgedehnter Schallquellen (F. Enxer), 
Versuch einer physikalischen Interpretation der kinstle- 


rischen Qualitaét musikalischer Klange (E. Turennavs), 


Musikalisch-akustische Kriterien von Konzertsdlen (F. 


Wincxet), Grofe Dirigenten — Randbemerkungen zum 
Problem der Begabung (H. P. Scumrrz) , Anwendung des 
MS-Stereophonieverfahrens zur Untersuchung raumaku- 
stischer Probleme (H. Lauripsen), Das stilisierte Real- 
gerausch in Horspiel und Film (W. Harr und J. Ascuen- 
BRENNER), Tonaufzeichnungsverfahren mit 1 bis 6 Kana- 
len (H. Fries), Kiinstlerische Probleme bei der Tonband- 
Musikaufnahme (J. UrsAnsxr), Methoden der Musik- 
aufnahme (P. Burxowi1rz), Stereophonie auf der Schall- 
platte (F. Scutecer), Uber den Neubau der Musik-Aka- 
demie in Warschau und ihre elektroakustischen Einrich- 
tungen (M. Rasewsxr). 

Ein Tagungsbericht wird herausgegeben. Die nachste 
Tonmeistertagung ist fiir den Herbst 1960 in Detmold 
geplant. H. Késrers 


Akustischer Kongre8 in Budapest 


Vom 4. bis 8. November 1957 veranstaltete der Opti- 
sche und Kinotechnische Verein in Budapest einen Kon- 
gre, auf dem Themen aus dem Bereich der reinen und 
angewandten Akustik sowie der Meftechnik behandelt 
wurden. Neben den Vortrégen wurden den Besuchern 
(ca. 100 aus Ungarn und 24 aus dem Ausland) meh- 
rere Besichtigungen und eine Ausstellung von MebB- 
geradten geboten. 

Der Kongre8 wurde von T. Tarnoczy erdffnet, der in 
seinem einleitenden Vortrag iiber die Bedeutung der 
akustischen Konferenzen sprach. 

I. P. Varx6, Budapest, machte mit der Theorie des 
Kardioid-Mikrophons bekannt. Mit Hilfe seiner Me- 
thode konnte er nicht nur die bisher bekannten Mikro- 
phontypen systematisieren, sondern er wies auch auf 
Méglichkeiten der bis jetzt noch unbeniitzten Mikro- 
phontypen hin. 


J. Mernavr, Prag, sprach iiber Intermodulationsver- 
zerrungen in Lautsprechern. Die objektive Qualitits- 
beurteilung der Lautsprecher wird durch seine Me8- 
methode und seine Nomogramme sehr erleichtert. 

F. TaxAcs, Budapest, diskutierte in seinem zusam- 
menfassenden Vortrag die Folgerungsverzerrung der 
Schallplattenwiedergabe und deren Zusammenhang mit 
der Spielzeit. 

Der Vortrag von G. Wacner, Prag, tiber analytische 
Methoden zur Lésung der mechanischen Schwingungs- 
systeme enthielt neue wertvolle Ergebnisse, die mit 
Hilfe der Kirchhoffschen Kombinationsregeln gewon- 
nen wurden. 

D. Huszty, Székesfehérvar, schlug als Qualifikations- 
angabe von Lautsprechern die in der Achse gemessene 
Schalldruckverteilung sowie die Messung der Gesamt- 
energie vor. 
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J. Kacprowsx1, Warschau, machte in seinem Vortrag 
Uber die Moglichkeiten der Abbildung der akustischen 
Impedanz des menschlichen Ohres mittels einer Ersatz- 
schaltung* mit der Planung und Verfertigung eines 
neuen und fiir Normalisierungszwecke geeigneten kiinst- 
lichen Ohres bekannt. 

P. J. Honort, Naerum, hielt einen zusammenfassen- 
den Vortrag tiber die Berechnung von Lautstarken. 

Verstandlichkeitsuntersuchungen in Filmstudios bei 
unterschiedlicher Mikrophonaufstellung und verschiede- 
nen Mikrophontypen wurden von F. Kemenzs, Budapest, 
behandelt. Die Messungen wurden auch bei Sprachauf- 
nahmen mit Begleitmusik und bei fertigen Filmen aus- 
gefihrt. Interessant ist die Bestimmung des Intelligenz- 
faktors der Logatomverstandlichkeit. 

W. Jassem, Warschau, sprach tiber das Durchschnitts- 


. spektrum der polnischen Sprache und verglich die Er- 


gebnisse mit den bei anderen Sprachen erhaltenen 
Durchschnittsspektren. Die Beachtung der statistischen 
Lautverteilung kann in seinem Untersuchungsverfahren 
als neu betrachtet werden. 

E. Maeyart, Budapest, berichtete iiber seine Unter- 
suchungen, die er im Zusammenhang mit der elektro- 
akustischen Regelung des Nachhalls ausfihrte. Er hat 
nachgewiesen, dafi die auf der elektrischen Seite ge- 
daimpften Lautsprecher einen Schluckgrad a=45 bis 
75% (abhangig von der Frequenz und dem Belastungs- 
widerstand) haben. 

M. Kwiex, Posen, bot eine schéne Zusammenfassung 
iiber die Zusammenhinge zwischen den rawmakustischen 
physikalischen Charakteristiken und den subjektiven 
musikalischen Empfindungseigenschaften. Bei einem 
Testkonzert wurden verschiedene akustische Eigenschaf- 
ten kiinstlich verdndert. Die bisher erhaltenen Ergeb- 
nisse, iiber die vorgetragen wurde, beweisen die Brauch- 
barkeit der Methode. 

T. Tarnoczy, Budapest, hielt einen Vortrag mit dem 
Titel ,,Planung und Untersuchung von Hallrdumen 
kleinerer Abmessungen“. Er wies nach, das man in 
einem geeignet ausgebildeten und nicht mit steifen 
Wanden begrenzten 10m* groBen Raum sowohl vom 
Gesichtspunkt der Nachhallcharakteristik und der Ener- 
gieverteilung, als auch der Frequenzirregularitat mit 
den groBen Salen gleichwertige Hallraume verfertigen 
kann. Ein weiterer Vorteil dieser Konstruktion besteht 
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darin, da} sich die Nachhallzeit fiir tiefe Téne nicht 
erhoht. Es ist interessant, da in diesem Raum nicht 
der von Scuréprer gefundene Zusammenhang zwischen 
der Frequenzirregularitat und der Nachhallzeit von 
F,,=1,4 T, sondern eine Beziehung F,,~ T zu messen ist. 

P. V. Briiex, Naerum, machte in seinem Vortrag iiber 
zeitgemdbe raumakustische Messungen mit den Auf- 
nahmen der Somervitreschen Kurvenscharen und den 
damit zusammenhangenden Auswertungsméglichkeiten 
bekannt. 

Frau M. LuxAcs, Budapest, untersuchte den inneren 
Gerduschpegel von Autobussen und parallel damit den 
Gehorverlust der Autobusfiihrer. Es ergab sich ein 
interessanter Vergleich zwischen dem von den verschie- 
denen Autobussen erzeugten Larm und dem Gehérver- 
lust der Autobusfiihrer. Neu scheint die Erkenntnis, 
da die Autobusfiihrer nicht auf ihrem rechten Ohr, 
auf der Motorseite, sondern auf ihrem linken Ohr, auf 
der Fensterseite, einen groéBeren Gehérverlust erleiden. 

A. Freunn, Wien, sprach iiber die elektroakustischen 
Einrichtungen der Bundestheater in Wien. 

P. Grueuss, Budapest, zeigte eine Hochfrequenz-Film- 
aufnahme tiber die in Fliissigkeiten erzeugten Stob- 
wellen. 

P. V. Briist, Naerum, behandelte in seinem Vortrag 
die Methoden der Schallanalyse bei Gerduschen in der 
Industrie und die aus den Messungen ableitharen Fol- 
gerungen. 

W. Ertur, Dresden, gab eine Zusammenfassung iiber 
die Theorie und die Mebergebnisse der Nebenweg- 
Schalliibertragung zwischen Raumen. Sein Vortrag bot 
eine gute Ubersicht iiber die neuesten Ergebnisse, die 
auf diesem Gebiet an verschiedenen deutschen For- 
schungsinstituten gewonnen wurden. 

G. Hecxenast, Budapest, machte mit einem neuen 
Studio-Magnettongerat fiir zwei Geschwindigkeiten be- 
kannt, welches vom ungarischen Rundfunk entwickelt 
wurde. 

E. Voet, Lahr (Baden), trug tiber neue elektro- 
akustische Me&gerate der ElektromeBtechnik vor und 
behandelte ein Magnetophon mit verdnderlicher Ge- 
schwindigkeit und die Kuntsche Nachhallplatte. 

Die Konferenz wurde mit einer Zusammenfassung 
von I. Barra geschlossen. 

T. Tarnoczy 


Société Internationale d’Audiologie 


IV° Congrés, Padova (Italia), 2—5 octobre 1958 


Les Secrétaires généraux (Dr. Circ et Dr. Trenque) 
de la Société Internationale d’Audiologie, nous font 
savoir que cette Société tiendra son [Ve Congrés a Pa- 
dova (Italia), les 2-3-4-5 octobre 1958. Cette trés im- 
portante réunion sera réservée aux audiologistes s'in- 
téressant particuliérement aux progrés de cette spécia- 
lité. Trois «Tables Rondes» sont prévues sur les sujets 
suivants: «Prothése auditive»; «Traumatismes sono- 
res»; «Interdépendance cochléaire». 


Toute correspondance pour ce congrés est a adresser au 

Secrétaire général du [Ve Congres de la Société Inter- 
nationale d’Audiologie, 

Prof. ArsLan, 

Clinique O.R.L. Universitaire, 

Padoya, Italia. 

Par ailleurs, c’est au printemps 1960, que se tiendra 
a Basel (Confédération helvétique), le Ve Congrés de 
cette Société. 


3. Internationaler Akustischer Kongre8 Stuttgart 1959 


Der 3. Internationale KongreB fiir Akustik findet in 
der Zeit vom 1.9 bis 8. 9. 1959 in Stuttgart (Deutsch- 
land) statt. Der Kongre8 wird auf Veranlassung der 
International Commission on Acoustics gema den Be- 
schliissen des 2. Kongresses 1956 (Cambridge) ver- 
anstaltet und vom Verband Deutscher Physikalischer 


Gesellschaften gemeinsam mit der Nachrichtentechni- 
schen Gesellschaft des Verbandes Deutscher Elektro- 
techniker und gemeinsam mit dem Verein Deutscher 
Ingenieure durchgefiihrt. 

Die Sitzungen werden in der Stuttgarter Liederhalle 
und in der Technischen Hochschule stattfinden. 
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Das Vortragsprogramm soll Themen aus dem gesam- 
ten Bereich der physikalischen und technischen Akustik 


umfassen, vorwiegend tiber: Larm- und Schwingungs- . 


abwehr, Elektroakustik und Raumakustik sowie phy- 
sikalische Akustik. 

Zusammenfassende Vortrage sind in Aussicht genom- 
men aus folgenden Gebieten: Akustische Methoden in 
der Festkérperphysik — Schallausbreitung in Gasen 
bei kleinen Drucken — Turbulenz und Schall — Kayvi- 
tation — Strukturforschung mit Ultraschall — Physio- 
logische Akustik — Psychologische Akustik — Stereo- 
phonie — Raumakustik — Bauakustik — Gerdausch- 
bekampfung — Informationstheorie — Sprache. 

Das Rahmenprogramm wird eine Fahrt nach Miin- 
chen einschlieSen mit der Moglichkeit, industrielle For- 
schungsstatten und kulturelle Sehenswiirdigkeiten zu 
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besichtigen, ferner Ausflugsfahrten in die weitere Um- 
gebung Stuttgarts (Schwarzwald, Bodensee). 


Eine Ausstellung zeigt moderne akustische und elek- ~ 


trische MeBgerate der Industrie. 

Auslandische Besucher, insbesondere  solche 
Ubersee, werden die Méglichkeit haben, vor oder nach 
dem KongreB deutsche Forschungsstatten zu besichtigen. 


Weitere Ankiindigungen tiber den Kongrefi werden ~ 


folgen, sobald die Einzelheiten festliegen. Interessen- 
ten werden schon jetzt gebeten, dem Sekretar der Ta- 
gung, Herrn Dr.-Ing. Eszrnarp Zwicker, Stuttgart-N, 
BreitscheidstraBe 3, ihre Adresse mitzuteilen, damit 
ihnen alle weiteren Nachrichten persénlich zugesandt 
werden kénnen. — Anmeldung von Einzelvortragen, die 
bis zum 15. Mai 1959 entgegengenommen werden, sind 
zu gegebener Zeit an die gleiche Adresse zu richten. 
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E. G. Ricnarpson, Technical aspects of sound, 
Vol. 2. Elsevier Publishing Co., Amsterdam 1957, 
410 pages, 206 figs, 23 cm X 16 cm, sh 70/—. 


The volume 2, edited by E. G. Ricuarpson, is the se- 
cond volume of three. 
The contributors to this volume are respectively 
1. H. H. Hussarp, Langley Aeronautical Laboratory 
(USA), 

2. L. W. Lassirer from the same Laboratory, 

3. E. Meyer, III. Physikalisches Institut, 
of Gottingen (Germany). 

4. B. E. Norrincx, Tube Investments Research Labora- 
tories, Cambridge (England), 

5. H. Oserst, Physikalisch-Technische Bundesanstalt, 

Braunschweig (Germany), 

6. A. A. Recmr, Langley Aeronautical 
(U.S.A.), 

7. E.G. Ricnarpson, King’s College, University of Dur- 
ham (England), 

8. K. Tamm, III. Physikalisches Institut, University of 

Gottingen (Germany). 

9. N. B. Terry, National Coal Board, Mining Research 

Establishment, Isleworth (England). 

These different contributors have been chosen each 
as a specialist in his own field. 

The book has been divided in four parts. The first 
chapter has been written by RicHarpson himself with 
the subject: Propagation of sound in the atmosphere 
and the sea. Thereafter follow three divisions: 
Division I: Ultrasonic transducers and appli- 
cations, written by B. E. Notrinex and N. B. Trrry. 
Division II: Underwater acoustics. This chap- 
ter is divided in four parts: 

1. Technique of measurements and apparatus, 
2. Air bubbles in water, written both by E. Meyer, 
3. Broad-band absorbers for water-borne sound, by 

K. Tamm, 

4. Resonant sound-absorbers, by H. Oxsersv. 

Division IIT: Aircraft noise and underwater 

propulsion noise, divided in four chapters: 

1. Propeller noise, by A.A. Recrer, H. H. Husparp 
and L. W. Lassirer, 

2. Jet noise, by the same authors, 

3. Other sources of aircraft noise and comparative noise 
levels, also by the same authors, 

4. Sounds of propulsion in water, by E. G. Ricuarpson. 


University 


Laboratory 


As we can see from the content of this book, it is 
an assembly of subjects on which practically no books 
have been written up to-date. Ricuarpson was fortunate 
that he was able to group a such excellent group of con- 
tributors. 

The book contains 407 pages and has been printed 
in a nice form. It is richly illustrated and it contains a 
large number of references on the different subjects. 

It is a book which gives a lot of information on the 
different recent problems in technical acoustics. It will 
certainly be in its right place in every laboratory on 
acoustics. A. van IrreRBEEK 


E.G. Ricuarpson, Relaxation spectrometry. 
North-Holland Publishing Co., Amsterdam 1957, 
140 pages, 88 figs, 15 cm x 23 em, sh 40/—. 


E. G. Ricuarpson is an active bookwriter. While in 


his earlier books, he treated essentially problems on : 


pure acoustics and ultrasonics in his recent book he 
discusses new problems. The subject of the book is con- _ 


general which appear during elastic deformations in 
solids, liquids, gases and. plastic substances. Acoustic 
relaxation problems are included. But at the end of 
the book in chapter VII the author discusses dielectric 
relaxation phenomena. 

Ricuarpson treats systematically the different phe- 


nomena and new fields of investigations which began to 


appear in the literature during the past years. At the 
end of each chapter an extended list of references of 
the newest literature is indicated. Ricuarpson’s book is 
devided into eight chapters. In the first chapter the author 
introduces the notion of relaxation on the basis of dif- 
ferent mechanical models. In the second chapter the 
relaxation phenomena are treated theoretically with 
the use of electromagnetic models. In chapter III are 


discussed the infrasonic properties of solids, together 


with the experimental methods. Chapter IV contains a 
description of the acoustic properties of solids (metals) 
in connection with internal friction. Chapter V is a dis- 
cussion of the ultrasonic relaxation phenomena in gases 
and liquids and finally in chapter VI is treated the ab- 
sorption in metals in connection with the experiments 
of Truett, Mackinnon, Mason etc. 

Finally it can be said the book has been edited in an 
agreeable form and it can certainly help young investi- 
gators who wish to start in this new field of research. 


A. van IrTERBEEK 


nected with the problem of relaxation phenomena in’ — 


aus” 
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